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level,  ASTM  slope,  and  density  are  presented  for  viscosity-pressure 
data_^These  correlations  are  tested  using  data  determined  by  a 
variety^of  investigators  over  a temperature  range  of  32°  to  275°F. 

The  correlations  developed  are  shown  to  require  less  physical 
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determined  using  a chemical  conversion  correlation  between  static  and 
dynamic  runs.  It  is  shown  to  be  higher  than  generally  used  flash 
temperatures.  Chemical  reactions  are  shown  to  play  an  important  role 
in  boundary  lubrication. 
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I.  fluid  and  lubricant  studies 


A.  SUMMARY . Studies  have  been  continued  to  understand  better  the  inter- 
action of  lubricants  and  wort  .ig  fluids  and  the  environment  in  which  these 
materials  must  function.  Film  thickness  of  the  lubricant  in  the  bearing  has 
particular  application  to  aircraft  uses  where  minimum  viscosity  levels  (high 
temperature  operation)  are  used  to  achieve  operability  similar  to  that  shown  at 
lower  temperatures.  Earlier  reports  include  the  development  of  correlations 
enabling  one  to  derive  a value  for  the  viscosity-pressure  coefficient  a.  The 
pressure  coefficient  determinations  carried  out  using  the  PRL  pressure  viscometer 
appear  to  be  pertinent  for  use  in  elastohydrodynamic  equations  for  the  measure- 
ment of  film  thickness.  Techniques  have  been  developed  which  allow  the  accurate 
determination  of  wear  debris  developed  using  the  four-ball  wear  tester.  Organo- 
metallic  materials  both  soluble  and  insoluble  in  the  test  oil  can  be  determined. 
Iron  and  iron  oxides  can  be  determined  in  the  presence  of  each  ether. 

Current  theory  indicates  that  pressure  viscosity  may  be  an  important 
parameter  for  hydrodynamic  and  elastohydrodynamic  lubrication.  Due  to  the 
difficulty  and  cost  required  to  obtain  experimental  data,  it  is  desirable  to 
develop  simple  pressure  viscosity  correlations  which  cover  a wide  range  of  fluid 
types . 


Gas  solubility  during  the  determination  represents  the  largest  error  in 
these  viscosity-pressure  data. 

The  objectives  of  this  study  are: 

1.  To  evaluate  the  quality  of  viscosity-pressure  data  obtained  from  the 
PRL  high  pressure  viscometer; 

2.  To  determine  the  effect  of  dissolved  gas  on  the  viscosity  pressure 
properties  measured  in  the  PRL  unit; 

3.  To  obtain  experimental  pressure  viscosity  values  on  a series  of  "dumb- 
bell" blends  of  low  and  high  molecular  weight  components  such  as  gas 
oil  bright  stock  and  paraffinic  resin-gas  oil  blends; 

4.  To  derive  correlations  for  viscosity-pressure  as  a function  of  fluid 
properties  from  the  PRL  pressure-viscosity  data; 

5.  To  test  the  correlations  with  the  available  viscosity-pressure  data 
from  other  sources  as  well  as  the  PRL  unit  data;  and 

6.  To  relate  the  correlating  parameters,  wherever  possible,  to  physical 
significance  and  or  theory. 


Gas  solobility  studies  have  been  used  to  determine  the  diffusion  effect  indicated 
by  viscosity  stability  as  a function  of  time  at  constant  pressure. 

The  pressure  viscosity  of  Pennsylvania  mineral  oils,  "dumbbell"  blends  of 
mineral  oil  fractions,  and  paraffinic  resin-solvent  blends  have  been  measured 
with  the  PRL  pressure  viscometer  at  pressures  up  to  3000  psig  and  temperatures 
of  323,  100°,  and  130°F.  The  results  show  that  the  pressure  coefficients  are 
independent  of  pressure  for  the  fluids  studied.  This  is  in  agreement  with 
previous  results  with  the  PRL  high  pressure  viscometer. 


Data  ter  pure  hydrocarbons,  mineral  oils,  nonhydrocarbons,  and  blends  of 
mineral  oils,  polymers  and  resins  at  100°F  were  used  to  establish  the  correlations. 
The  final  form  ot  the  correlations  for  viscosity-pressure  data  were  obtained  by 
regressional  analysis.  The  resultant  empirical  relationships  involve  viscosity 
level,  ASTM  Slope  and  density.  These  parameters  appear  to  be  related  to  such 
physical  characteristics  as  molecular  interlocking,  molecular  packing,  molecular 
rigidity,  and  size  of  the  flow  unit.  In  all  cases  the  parameters  used  in  the 
correlations  are  consistent  with  the  void  volume  theory  proposed  by  Cohen  and 
Turnbull  (1) . 

The  resultant  correlations  were  not  only  tested  with  viscosity-pressure  data 
at  100°F  from  a variety  of  investigators  but  were  also  tested  with  data  at 
temperatures  ranging  from  32°  to  275CF.  The  current  correlations  when  compared 
with  the  other  leading  correlations  were  shown  to  require  less  physical  property 
data,  apply  to  a wider  range  of  chemical  compositions,  cover  a wider  range  in 
fluidity  and  provide  better  accuracy. 

The  mechanism  of  boundary  lubrication  is  important  for  the  design  of  a 
lubr icant-bear ing  system  Wear  debris  analysis  has  been  used  to  measure  the 

producticn  of  metal  particles,  metal  oxides,  and  organometallic  reaction  products. 
The  separation  of  metal  particles  and  metal  oxides  can  be  done  by  preferentially 
dissolving  the  iron  oxides  using  an  acr id ine- inhib it ed  hydrochloric  acid  solution. 
The  iron  can  be  preferentially  dissolved  from  the  iron  oxides  by  the  use  of  nitric 
acid.  The  iron- iron  oxide  separation  can  be  measured  both  ways  with  relatively 
good  results. 

The  qualitative  results  in  the  four-ball  wear  tester  operated  in  the  normal 
way  (600  rpm,  40  kilograms,  167  F)  with  a hydrocarbon  lubricant  have  been 
analyzed.  The  results  show  that  in  all  cases  there  is  seme  iron  oxide  and  some 
organometallic  formed  in  the  wear  conjunction.  These  reactions  would  be  predicted 
for  high  temperatures  of  the  order  of  those  measured  in  elastohydrody namic  lubri- 
cation or  the  higher  temperatures  predicted  in  boundary  lubrication 

A series  of  four-ball  wear  tests  have  been  conducted  at  low  speed  to  limit 
the  heat  generated  at  the  bearing  conjunction  In  these  tests,  the  bulk 
temperature  has  been  used  to  control  the  bearing  temperatures.  These  tests  with 
hydrocarbon  lubricants  have  shewn  low  levels  of  chemical  activity  at  low  bulk 
temperatures.  On  the  other  hand,  there  is  a sharp  rise  in  chemical  activity  at 
higher  bulk  temperatures.  The  chemical  activity  to  predute  c t ganometall ic 
products  is  substantial  at  bulk  temperatures  of  500‘F  and  higher. 

There  is  clearly  a temperature  effect  n chemical  reu  iivity  that  follows 
the  bulk  temperature  in  a manner  similar  to  the  results  of  microcorrosion  studies. 
The  micrccor rosion  studies  are  conducted  in  a static  system  where  bulk  temperature 
is  t e only  source  of  energy  for  the  reaction.  One  of  the  ptcblems  in  comparing 
the  slow  speed  fout-ball  with  the  high  speed  four -ball  is  the  area  to  consider. 

The  high  speed  four-ball  test  at  low  bulk  temperature  focuses  the  chemical 
reaction  on  the  wear  scar  areas.  For  the  high  bulk  temperature,  low  speed  four- 
ball  tests  the  total  ball  area  may  be  more  important  than  the  wear  scars.  In 
order  to  determine  the  relative  importance  of  the  total  ball  area  and  the  wear 
scar,  low  speed  four-ball  tests  have  been  run  without  contact  between  the  balls. 
These  data  are  compared  with  the  high  speed  four-bali  and  the  micrccorrosion  tests. 
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The  information  from  microcorrosion  studies  indicates  that  the  reaction 
between  the  oxidized  hydrocarbon  and  the  iron  or  iron  oxide  surface  has  an 
initial  rate  that  shows  a typical  temperature  dependence.  In  the  case  of  the 
four -ball  wear  tester,  the  wear  scars  are  continuously  wiped  by  the  wear  action 
and  provide  a continuing  active  surface. 

Reduced  reaction  with  growing  wear  scar  and  decreasing  load  is  demonstrated 
by  sequential  tests  with  the  same  set  of  balls.  This  result  suggests  that  the 
temperature  of  the  reacting  environment  is  changing  with  time  and  actual  unit 


load  in  the  bearing  area.  Estimates  of  reaction  kinetics  are  based  on  the 
assumptions  of  no  change  in  the  reactivity  of  the  surface  and  an  average  wear  scar 
area  and  temperature.  It  is  apparent  that  the  wear  scar  area  and  surface 
temperatures  are  changing  continuously. 

Data  for  estimating  reaction  rates  and  temperatures  are  obtained  using 
standard  and  sequential  four-ball  wear  tests  with  a super -ref ined  mineral  oil. 

The  junction  temperature  rise  determined  using  a chemical  conversion  corre- 
lation between  static  and  dynamic  runs  is  shown  to  be  substantially  higher  than 
the  generally  used  Blok-Archard  flash  temperatures. 
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B.  VISCOSITY-PRESSURE  CORRELATION  OF  LIQUIDS.  Bearing  pressures  of  the 
order  of  50,000  to  500,000  psig  are  encountered  in  elastohydrodynamic  and  other 
thin  film,  fluid  lubricant  applications.  As  a result,  many  of  the  viscosity- 
pressure  studies  with  liquids  have  been  conducted  in  this  high  pressure  range. 
However,  the  recent  studies  by  Bell  (2)  and  Cheng  (3)  on  the  elastohydrodynamic 
theory,  suggest  that  the  inlet  film  thickness  is  largely  governed  by  the  pressure- 
viscosity  effect  in  the  low  pressure  region.  Westlake  (4)  further  pointed  out 
that  the  best  correlation  on  film  thickness  seems  to  be  obtained  with  viscosity 
data  obtained  at  10,000  psig.  It  is  the  current  theory  that  the  most  important 
viscosity-pressure  characteristic  of  liquids  is  the  pressure-viscosity  coefficient 

, where  U = viscosity,  p = pressure.! 

p=0  J 

Advanced  hydraulic  systems  typical  of  those  used  in  aircraft  and  aerospace 
applications  operate  at  bulk  pressures  of  up  to  10,000  psig.  Precise  viscosity 
properties  in  this  pressure  range  are  required  to  study  the  fluid  flow  in  hydraulic 
systems . 

Host  of  the  high  pressure-viscosity  data  were  obtained  with  rugged  equipment 
capable  of  achieving  pressures  of  100,000  psig  or  higher.  The  ruggedness  of  most 
of  these  viscometers  resulted  in  low  precision  data  at  low  pressures  (0  to  10,000 
psig).  The  development  of  a low  pressure  precision  capillary-type  pressure 
viscometer  by  this  laboratory  is  discussed  in  previous  reports.  The  use  of  this 
viscometer  to  measure  the  viscosity-pressure  characteristics  of  a wide  variety 
of  fluid  types  is  presented  in  Annual  Report  AFML-TR- 67-101 , Part  I.  The  use  of 
the  same  pressure  viscometer  to  measure  the  effect  of  polymers  and  resin  on  the 
viscosity-pressure  characteristics  of  the  solvents  is  included  in  Annual  Report 
AFML-TR-70- 304 , Part  II.  It  is  desirable  to  obtain  a general  pressure  coefficient 
correlation  to  include  both  polymeric  and  nonpolymer  fluids. 

1 . Review  of  Methods  for  Determining  and  Predicting  Viscosity-Pressure 
Properties.  Viscosity-pressure  data  for  liquids  are  relatively  limited  in  the 
literature  because  of  the  relatively  complex  apparatus  required  for  measurement. 
Pertinent  literature  concerning  the  experimental  determinations  and  the  corre- 
lations of  the  viscosity-pressure  characteristic  of  liquids  are  summarized  below. 

a.  Experimental  Determination.  There  are  basically  four  types  of  low 
shear  viscometers.  Bridgman  (5),  Kuss  (6),  and  the  ASME  Project  (7)  used  the 
falling  weight  viscometer.  A rolling  ball  viscometer  was  used  by  Webb  and  Dixon 
(8)  These  viscometers  were  designed  to  obtain  data  at  very  high  pressures. 
Pressures  up  to  170,000  psig  are  readily  obtainable  with  the  apparatus  and  with 
modification  pressures  exceeding  400,000  psig  may  be  reached.  However,  for 
pressures  below  10,000  psig,  this  apparatus  doesn't  produce  high  accuracy  pressure 
viscosity  data. 

By  determining  the  change  in  the  resonant  frequency  and  electrical 
resistance  of  a longitudinally  vibrating  probe,  a commercial  Instrument  called 
Utraviscoson  (9)  can  be  used  to  determine  viscosity-pressure  effects  at  low  shear 
rates.  The  maximum  pressure  for  this  apparatus  is  only  1000  psi. 

Data  obtained  using  a capillary  type  pressure  viscometer  are  discussed  in 
previous  Annual  Reports.  The  PRL  pressure  viscometer  was  specifically  designed 
to  obtain  highly  accurate  data  for  pressures  below  10,000  psig  and  under  low 
shear  rate. 
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There  are  also  three  major  types  of  high  shear  pressure  viscometers.  The 
problems  connected  with  the  high  shear  viscometers  are  that  it  is  difficult  to 
separate  the  shear,  pressure,  and  temperature  effects. 

Novak  and  Winer  (10)  constructed  a two-way  high  pressure  capillary  viscometer. 
Although  it  covers  a shear  stress  range  of  300  to  1.1  x 106  dynes  per  square 
centimeter,  the  uncertainty  of  the  pressure  level  increases  with  the  decreasing 
shear  stress.  Therefore,  data  obtained  by  this  apparatus  under  low  shear  stresses 
give  the  lowest  level  of  accuracy. 

Appeldoorn,  Okrent,  and  Philippoff  (11)  modified  the  Mason  (12)  vibrating 
crystal  viscometer  for  viscoelastic  behavior  studies  of  polymer  blended  solutions. 
Viscosity  and  elasticity  are  determined  by  the  change  in  the  resonant  frequency 
and  electrical  resistance  of  the  probe.  Limited  by  the  flatness  of  the  resonance 
peak  the  highest  atmospheric  viscosity  fluid  tested  was  69.28  centipoises  at  100°F. 

Optical  elastohydrodynamic  fluid  testing  was  developed  by  Cameron  (13,14) 
and  co-workers.  Optical  interference  was  used  to  measure  the  film  thickness 
between  a glass  or  quartz  plate  and  a specially  finished  rolling  steel  ball  which 
is  loaded  against  the  plate.  By  calibrating  the  apparatus  with  fluids  of  known 
pressure-viscosity  characteristics,  the  pressure-viscosity  of  unknown  fluids  can 
be  measured.  The  operational  range  of  the  rate  of  shear  is  between  0.1  and  40 
million  sec-1. 

Other  than  the  work  discussed  in  AFML-TR-70-304,  Part  IX  using  polymer 
solutions,  all  the  low  shear  viscosity-pressure  studies  were  with  pure  hydro- 
carbons, mineral  oils  and  synthetic  oils. 

From  these  studies,  the  following  generalizations  can  be  made: 

1.  Plots  of  the  log  of  viscosity  versus  pressure  up  to  10,000  psig  at 
constant  temperature  result  in  a straight  line  for  almost  all  fluids.  Some 
notable  exceptions  include  a slight  upward  curvature  for  phenylethers  and 
chlorinated  hydrocarbons; 

2.  Viscosity  is  a function  of  both  temperature  and  volume; 

3.  In  general,  temperature  effects  are  less  at  high  pressure; 

4.  For  the  lower  molecular  weight  range,  chain  length  has  a significant 
effect  on  the  pressure-viscosity  characteristics.  However,  the  effect  of  molecular 
structure  is  very  much  greater.  For  the  long  chain  molecule,  viscosity  level 
increases  with  chain  length  whereas  the  viscosity-pressure  characteristics 
remains  constant; 

5.  The  increased  internal  mobility  of  the  molecule  (e.g.,  building  ether 
linkage  into  a long  chain  molecule)  decreases  the  dependency  of  the  viscosity  on 
pressure; 

6.  The  viscosity-pressure  effect  increases  with  the  introduction  of  rigid 
or  bulky  side  groups; 

7.  Introduction  of  certain  functional  groups  (e.g..  Cl  or  OH)  can  increase 
the  viscosity-pressure  effect; 
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8.  The  largest  pressure  effects  are  for  those  substances  with  the  most 
complicated  molecular  structure; 

9.  For  a close  approximation,  the  effect  of  nitrogen  gas  saturation  on  the 
pressure  viscosity  of  fluids  is  to  reduce  the  pressure  viscosity  to  its  atmospheric 
level;  and 

10.  Reduced  viscosity  is  the  ratio  of  solution  viscosity  over  the  solvent 
viscosity  at  one  atmosphere.  The  effect  of  polymers  and  resin  on  the  viscosity- 
pressure  characteristics  of  the  solvent  is  the  same,  up  to  a reduced  viscosity 
level  of  twelve.  For  a reduced  viscosity  level  above  twelve,  the  effect  of 
polymers  on  the  viscosity-pressure  characteristics  is  a function  of  polymer  type. 

The  only  extensive  work  on  the  viscosity-pressure  effect  of  binary  mixtures 
lias  been  done  by  Dow  (15).  For  mixtures  of  molecules  with  linear  structure,  the 
viscosity-temperature  and  viscosity-pressure  coefficients  of  the  components  are 
additive.  However,  for  complex  molecules,  the  viscosity  isobars  of  the  mixtures 
display  variations  with  concentration.  The  deviation  from  the  additive  rule 
always  increases  with  the  pressure. 

Griest,  Webb  and  Schiessler  (16)  studied  three  hydrocarbon  mixtures.  The 
viscosity  of  these  compounds  is  some  additive  function  of  their  constituent 
groups  whether  these  groups  are  combined  in  the  same  or  different  molecules,  as 
long  as  the  basic  molecular  symmetry  is  unchanged. 

b.  Available  Prediction  Techniques.  Several  attempts  have  been  made 
to  calculate  the  viscosity  of  a liquid  under  pressure  in  the  absence  of  experi- 
mental data.  All  have  met  with  very  little  success.  Therefore,  the  major 
predictive  effort  has  been  directed  toward  empirical  correlations.  Since 
pressure-viscosity  data  are  difficult  to  obtain,  any  correlations  requiring  such 
data  have  very  limited  use.  This  survey  includes  only  those  correlations  which 
do  not  require  pressure-viscosity  data. 

A definition  of  a is  given  in  Equation  [1]. 

a = i log  ^ [1] 

o 

where: 


a = viscosity-pressure  coefficient, 
cm2/kg, 

p = pressure,  kg/cm2, 

hp  = viscosity  at  pressure  p,  centipolse,  and 

H = viscosity  at  atmospheric  pressure, 
centipolse . 


By  assuming  that  a is  constant  for  a particular  liquid  and  is  not  a function  of 
pressure,  Worster  (17)  proposed  the  following  linear  relationship  between  a and 
log  Do  for  naphthenic  mineral  oils. 

a = 0.006  + 0.00097  log  n 

o 


[2] 


where : 


a = viscosity-pressure  coefficient, 
cm 2 /kg,  and 

ri  = viscosity  at  atmospheric  pressure, 
centipoise. 

The  scatter  for  a is  ± 20  percent. 

Hartung  (18)  proposed  a graphical  correlation  to  predict  the  pressure 
viscosity  isotherms  of  liquid  lubricants  at  100°  and  210°F  up  to  100,000  psig. 
Data  required  to  make  the  prediction  are  atmospheric  pressure  densities  and 
viscosities  at  100°  and  210°F.  This  method  has  an  accuracy  of  ± 20  percent. 
However,  errors  in  this  method  are  increased  for  oils  which  have  less  than  four 
centistokes  viscosity  at  210°F  or  have  viscosity  indices  of  less  than  -100. 

Clark  (19)  assumed  that  for  all  oils  the  pressure  change  required  to  dupli- 
cate the  effect  of  any  specified  temperature  change  is  the  same.  As  a result, 
the  ASTM  Kinematic  Viscosity-Temperature  Chart  D341  may  be  adapted  for  viscosity- 
pressure  properties  by  replacing  each  temperature  value  on  the  scale  with  its 
equivalent  pressure  value.  This  method  usually  gives  results  reliable  to  only 
± 50  percent  even  after  classifying  oils  as  paraffinic,  etc. 

Lockhart  and  Lenoir  (20)  presented  a viscosity-pressure  correlation  for  high 
molecular  weight  pure  hydrocarbons.  Lubricating  oils,  neutral  stocks,  bright 
stocks,  and  distillates.  This  method  requires  a knowledge  of  the  atmospheric 
pressure-viscosity  at  the  temperature  of  interest  and  of  the  Watson  character- 
ization factor.  This  factor  is  defined  as: 


where : 

K = Watson  characterization  factor, 
w 

T^  = mean  average  boiling  point  of  the 
petroleum  fraction  in  °R,  and 

= density  at  60°F  and  atmospheric 
pressure  in  grams/cm3. 
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The  accuracy  of  the  correlation  decreases  as  the  pressure  increases.  At  a 
pressure  of  10,000  psig,  and  for  viscosities  of  less  than  1000  centipoises,  the 
average  error  is  ± 9 percent. 

Chu  and  Cameron  (21)  gave  a correlation  which  relates  pressure  coefficient, 
a,  to  atmospheric  viscosity  and  temperature.  The  correlation  predicts  that  at 
constant  temperature,  a decreases  with  viscosity  level.  This  contradicts  experi- 
mental data. 

Dixon  and  Webb  (22)  presented  a monograph  for  constant  viscosity  as  a 
function  of  temperature  and  pressure.  This  method  requires  the  knowledge  of 
viscosity  as  a function  of  temperature  to  predict  the  pressure-viscosity.  The 
average  error  is  reported  to  be  ± 20  percent  at  pressures  up  to  25,000  psig. 


The  first  quantitative  correlation  between  viscosity-temperature-pressure 
dependence  and  the  chemical  constitution  of  mineral  oils  was  obtained  by 
Roelands,  Vlugter,  and  Waterman  (23).  The  correlation  can  be  expressed  as 
fol lows: 


where 


viscosity  at  pressure  p,  centipoise, 
viscosity  at  atmospheric  pressure,  centipoise, 
pressure,  psig, 

percentage  of  carbon  atoms  in  aromatic  ring  structure, 
percentage  of  carbon  atoms  in  naphthenic  ring  structure,  and 
an  empirical  constant. 


The  carbon  percentages,  C,. , are  determined  by  the  n-d-M  method  (2A) 


where 


n = refractive  index 
d = density,  and 
M = molecular  weight 


The  exponent  y has  been  empirically  related  to  the  carbon  percentages  by 
the  following  equation: 


log  (y  - 0.890)  = 0.00855  (C.  + 1.5  CH)  - 1.930 


A maximum  deviation  of  five  percent  is  reported  for  values  of  log  h0  ranging 
from  about  -0.1  to  1.5,  for  temperatures  from  25°  to  90°C,  and  for  pressures  up 
to  1000  atmospheres. 


By  combining  the  Slotte-Fahrenheit  equation  for  the  viscosity-temperature 
relationship  and  the  Barus  equation  for  the  viscosity-pressure  relationship, 
Appeldoorn  (25)  obtained  a four-constant  viscosity-temperatur e-pressure  equation 


where 


viscosity  at  any  temperature  and  pressure, 

viscosity  at  atmospheric  pressure  and  reference  temperature,  t , 
pressure, 

any  temperature,  °F, 
reference  temperature,  °F, 

viscosity-temperature  coefficient  at  atmospheric  pressure, 
pressure-viscosity  coefficient  at  the  reference  temperature,  and 


c = constant  which  determines  how  much  the  viscosity-temperature 
coefficient  increases  with  pressure  or,  conversely,  how  much 
the  viscosity-pressure  coefficient  decreases  with  temperature. 

This  equation  was  claimed  to  be  accurate  (<  10%  error)  for  95  percent  of  the 
data  tested  over  the  temperature  range  of  70°  to  250°F  and  up  to  15,000  psi. 

in  order  to  eliminate  the  two  pressure-viscosity  data  points  required  to 
evaluate  the  constants  for  the  equation,  Appeldoorn  assumed  that  the  relationship 
between  constants  c and  b and  between  constants  b and  a are  linear  and  obtained 
the  following  two-constant  empirical  correlation. 

log  = a (log  - 0.0234  p + 0.0302  p log  ■—)  [7] 

o o o 


The  two  constants,  a and  r|  , can  be  determined  from  two  atmospheric  viscosity 
data  points  at  temperatures  t and  t.  The  two-constant  equs.un  is  not  as 
accurate  as  the  four-constant  equation. 


Kouzel  (26)  assumed  that 


fldn 


ver  sus 


in  dpjp=o 

log  r)|p-o,  and  log  n versus  p are  all  linear. 


1 dn  , 

In  dt J p*o ’ 08 


1 dn) 

— 1 versus 

n dtjp^o 


The  resulting  equation  is 


nD  p 0.278 

i08  f = Tcfeo  (0.0239+0.01638  nQ 
o 

where: 

Hp  = viscosity  at  pressure  p and  temperature  T,  centipoise, 
nQ  = viscosity  at  atmospheric  pressure  and  temperature  T,  and 
p = pressure,  psig. 


[8] 


The  equation  is  correlated  for  the  viscosity  range  of  0.5  to  200  cer  ipoises  and 
up  to  5000  psig  and  425  = F.- 


Fresco,  Klaus  and  Tewksbury  (27)  developed  a graphical  relationship  involving 
the  logarithm  of  viscosity,  the  pressure  coefficient,  and  the  viscosity- 
temperature  properties  (ASTM  Slope)  of  the  fluid.  Second-degree  polynomial 
equations  were  also  developed  to  try  to  reproduce  the  graph.  The  pressure 
coefficient  is  defined  as: 


log  ^ = Paa  (10  4) 

where : 

Up  = kinematic  viscosity  in  centistokes  at  pressure  p, 

UQ  * kinematic  viscosity  in  centistokes  at  atmospheric  pressure, 
p = pressure,  psig, 

a = pressure  coefficient,  psig-1,  and 


[9] 


a = temperature  correction,  560/ °R 
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When  the  resultant  general  correlation  was  tested  with  data  from  mineral  oils 
pure  hydrocarbons  and  nonhydrocarbons,  the  average  deviation  was  ± 11  percent 


Roelands  (28)  made  the  assumption  that  sixty  percent  of  the  solvent  viscosity 
pressure  index  contributes  toward  the  blend  viscosity-pressure  index. 


where 


viscosity-pressure  indexes  of  the  base  oil 
and  blend  respectively,  and 

hypothetical  viscosity-pressure  index  of 
the  blend  read  from  the  chart. 


The  above  equation  requires  the  knowledge  of  a total  of  four  atmospheric 
viscosities;  two  for  the  blend,  and  two  for  the  base  oil.  For  some  base  oils 
density  and  refractive  index  at  104°C  may  also  be  required. 


After  calculating  the  viscosity-pressure  index  of  the  blend,  the  following 
equation  is  used  to  calculate  the  pressure-viscosity. 


log  (log  n + 1.2000)  = Z log  1 + 


dynamic  viscosity  at  pressure  p and  atmospheric  pressure 
respectively  at  the  temperature  of  interest,  centipoise. 


The  pressure-viscosity  studies  discussed  in  Annual  Report  AFML-TR-70-304 
Part  II,  separated  the  pressure  coefficient  of  the  polymer  solution  into  two 
parts. 


viscosity  of  polymer  solution  at  pressure  p,  and  at 
atmospheric  pressure  respectively,  centistokes, 

pressure,  psig, 

pressure  coefficient  of  solvent,  psig  ^ , and 
pressure  coefficient  of  polymer,  psig 


This  work  correlated  the  pressure  coefficients  of  polymer  solutions 
graphically  as  a function  of  the  relative  viscosities.  The  relative  viscosity 
is  the  ratio  of  the  viscosity  of  the  polymer  solution  to  the  viscosity  of  the 
solvent  at  atmospheric  pressure.  The  graph  can  also  be  represented  by  second 
degree  polynomial  equations.  The  value  cc  the  solvent  pressure  coefficient  can 
be  predicted  by  one  of  the  available  methods. 
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2.  Test  Apparatus  and  Procedures.  The  design  and  operation  of  the  PRL 
pressure  viscometer  have  been  described  in  detail  in  former  Annual  Reports  and 
in  a paper  by  Klaus,  Johnson,  and  Fresco  (29).  The  overall  viscometer  system 
is  shown  in  Figure  1. 

a.  Pressure  Supply  and  Measurement  Equipment.  The  pressurized  nitrogen 
gas  bottle  (A)  is  used  to  obtain  pressures  from  0 to  2000  psig.  Additional 
pressures  up  to  10,000  psig  are  obtained  with  the  hydraulic  booster  system 
(B,C,D).  Two  pressure  gages  (E,  5000  psig,  and  F,  10,000  psig)  are  connected  in 
parallel  to  provide  adequate  sensitivity  over  the  entire  pressure  range.  The 
5000  psig  gage  has  a precision  of  ± 5 psig.  Whereas,  the  10,000  psig  gage  has  a 
precision  of  ± 10  psig.  Both  gages  are  calibrated  using  a dead  weight  tester 
which  has  an  accuracy  of  ± 1 psig.  The  calibration  data  are  shown  in  Table  1. 

b.  Constant  Temperature  Apparatus.  A Pyrex  glass  jar  twelve  inches  in 
diameter  and  twenty  inches  deep  is  used  as  the  constant-temperature  water  bath 
(1),  shown  in  Figure  1,  for  determinations  at  100°  and  130°F.  For  32°F  deter- 
minations, an  ice-water  slurry  is  used. 

For  100  and  130CF  determinations  the  bath  is  controlled  by  a magnetic  roto- 
set  type  low  voltage  thermoregulator  in  conjunction  with  a low  voltage  relay 
system  which  controls  the  power  input  to  an  electric  immersion  heater.  Air  is 
continuously  bubbled  through  the  water  to  insure  even  temperature  distribution 
in  the  bath.  The  temperature  of  the  bath  can  be  controlled  to  within  ± 0.05°F. 

An  air  driven  stirrer  connected  above  the  bath  is  used  to  stir  the  ice-water 
slurry.  Screens  are  set  up  to  prohibit  the  flow  of  ice  past  the  sight  glass  of 
the  Jerguson  gage  without  prohibiting  the  flow  of  water.  Bath  temperature  can 
be  maintained  to  within  0,2rF  from  the  desired  32°F  temperature. 

c.  Pressure  Vessel.  The  pressure  viscometer  (I  in  Figure  1)  is  a 
modified  21-T-50  Jerguson  level  gage.  The  gage  has  a pressure  rating  of  10,000 
psig. 


The  capillary  unit  (Figure  2)  used  to  measure  the  viscosity  of  fluids  is 
attached  to  the  top  closure-  (K)  by  a Teflon  compression  cone  (see  Figure  3). 

The  whole  capillary  unit  is  placed  inside  the  pressure  vessel  (J). 

The  series  of  valves  and  passages  on  the  viscometer  are  designed  to  permit 
the  manipulation  of  the  fluid  from  reservoir  to  the  efflux  bulb.  By  proper  valve 
sortings  the  fluid  can  then  be  allowed  to  flow  from  the  efflux  bulb  under  the 
head  of  test  fluid  represented  by  the  difference  in  height  between  the  liquid 
levels  in  the  efflux  bulb  and  the  reservoir.  The  liquid  levels  in  the  capillary 
and  the  reservoir  can  be  observed  through  the  series  of  one-half  inch  diameter 
"bull's  eye"  windows  on  both  sides  of  the  gage  as  indicated  in  Figure  1. 

d.  Viscometer  Capillaries.  A detailed  drawing  of  a viscometer  capillary 
is  shown  in  Figure  2.  The  upper  etched  line  above  the  storage  bulb  (VJ)  is 
designated  as  the  fill  line.  The  viscometer  capillary  is  filled  to  this  line  for 
each  viscosity  determination  to  insure  a constant  drainage. 

The  bottom  bulb  (X)  is  denoted  as  the  efflux  bulb  and  the  efflux  time  is  the 
time  required  for  the  fluid  to  pass  from  the  etched  line  on  the  top  of  bulb  (X) 
to  the  one  below  the  bulb  (X). 
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The  small  bulb  (W)  is  a storage  bulb.  Its  purpose  is  to  allow  time  for 
the  manipulations  required  at  the  start  of  a test  and  also  to  insure  steady- 
state  conditions  before  the  measured  flow  begins. 

The  bottom  of  the  capillary  is  trumpet  shaped  to  minimize  the  kinetic 
energy  effects.  A series  of  capillaries  ranging  in  inside  diameter  from  0.9  to 
5.0  millimeters  are  available  to  obtain  efflux  times  for  viscosities  ranging 
from  0.4  to  10,000  centistokes. 

e.  Calibration  of  a Viscometer  Capillary.  The  viscometers  are 
calibrated  at  atmospheric  pressure  with  fluids  of  known  viscosity  and  are  used 
subsequently  as  relative  instruments  to  determine  the  viscosity  of  unknown 
liquids  in  the  general  fashion  of  modified  Ostwald  viscometers.  The  flow 
through  the  capillary  viscometer  can  be  adequately  described  by  Poiseuille's 
law  expressed  in  terms  of  efflux  time. 


where : 


t = 


8LV 

"P1 

Tgr4h 

prpg 

prp8 


K = 
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4, 

TTgr  h 
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efflux  time,  sec, 

capillary  length,  cm, 

3 

volume  of  liquid,  cm  , 

gravitational  constant,  980.6  cm/sec, 

capillary  radius,  cm, 

height  of  liquid  from  reservoir  level 
to  the  top  of  the  efflux  bulb,  cm, 

3 

density  of  liquid,  gm/cm  , 

3 

density  of  gas,  gm/cm  , 
kinematic  viscosity,  cm,  and 

3 

function  of  viscometer  geometry,  cm  /sec. 


When  the  viscometer  is  calibrated  at  atmospheric  pressure,  the  gas  density, 
pg,  can  be  neglected  since  it  is  on  the  order  of  0.1  percent  of  the  liquid 
density,  p^.  Then  the  Equation  [14]  becomes: 

t = K u [15] 

a 

where : 

t = efflux  time,  sec, 


function  of  viscometer  geometry,  cm  /sec,  and 

apparent  viscosity  from  atmospheric  calibration 
centistokes . 


Therefore,  at  atmospheric  pressure,  the  apparent  viscosity  is  the  atmospheric 
viscosity.  In  the  case  of  pressure  viscosities,  however,  the  density  of  the 
gas,  pgp  becomes  a substantial  factor.  In  order  to  use  the  atmospheric  cali- 
bration in  the  pressure-viscosity  determination,  the  following  relationship  is 
used . 


where 


true  viscosity  at  pressure,  p,  centistokes, 

apparent  viscosity  from  atmospheric  calibration 
centistokes , 

3 

density  of  liquid  at  pressure,  p,  gm/cm  , and 

3 

density  of  gas  at  pressure,  p,  gm/cm  . 


Values  of  the  density-pressure-temperature  relationship  for  nitrogen  were 
obtained  by  Sage  and  Lacy  (30).  Wright  (31)  developed  a series  of  charts  to 
predict  the  dens ity-pressure- temperature  relationships  for  petroleum  oils  and 
pure  hydrocarbons.  The  correlation  covers  a temperature  range  of  0°  to  500°F 
and  pressures  up  to  100,000  psig.  An  average  error  of  less  than  one  percent  is 
obtained  from  this  method  The  only  required  data  are  the  density  of  the  liquid 
at  atmospheric  pressure  and  the  desired  temperature.  These  data  have  also  been 
curve  fitted  into  equations  which  show  an  average  error  of  about  1.7  percent. 

As  the  critical  point  is  approached,  the  error  may  be  as  high  as  five  percent 
(32).  In  viscosity-pressure  studies  of  polymer  solutions  by  this  laboratory, 
these  same  relationships  are  used  to  predict  the  density-pressure-temperature 
propertiet  for  the  polymer  and  high  molecular  weight  petroleum  resin  solutions 
used  (AFML-TR- 70-304,  Part  II). 


At  constant  temperature  and  pressure,  calibration  constant,  K,  is  only  a 
function  of  the  height  of  the  liquid  from  reservoir  level  to  the  top  of  the 
efflux  bulb,  h.  Since  the  calibration  constant,  K,  is  a strong  function  of 
capillary  radius,  r,  which  is  a weak  function  of  temperature  and  pressure, 
calibration  constant,  K,  is  also  a weak  function  of  temperature  and  pressure. 


From  Equation  [14] 


-T 


...  dK  A 

AK  = — Ar 

dr 


= - 4 


ngr5h 


= - 4 k v* 

c [_r J 

K - K_  = - 4 K — 
2 c [_r 

K ■ Kc  [t]) 


When  Ar/r  is  expressed  as  the  glass  expansion  coefficient  and  is  assumed  to  be 
constant , 


K = K ll  - 7.47  x 10  (Tc-T)J  [20] 

When  Ar/r  is  expressed  as  the  glass  compressibility  and  is  assumed  to  be  constant, 


K = K [1  - 8.33  x 10  (p-p  ) 

c c 


where : 


K = real  function  of  viscometer  geometry,  cm  /sec, 

Kc  = function  of  viscometer  from  calibration,  cm^/sec, 

L = capillary  length,  cm, 

V = volume  of  liquid,  cmd, 

g = gravitation  constant,  980.6  cm/sec^, 

r = capillary  radius, 

h = height  of  liquid  from  reservoir  level  to  the  top  of  the 
efflux  bulb,  cm, 

Tc  = calibration  temperature,  °F, 

T = real  temperature,  °F, 

pc  = calibration  pressure,  psig,  and 

p = real  pressure,  psig. 

It  is  obvious  from  Equations  [20]  and  [21]  that  the  calibration  constant,  K, 
can  be  considered  constant  for  a wide  range  of  temperature  and  pressure. 

Figure  4 shows  a typical  calibration  curve  with  data  at  100°  and  130°F.  The 
calibration  curve  has  a curvature  rather  than  being  straight,  because  the 
radius  of  the  reservoir  varies  with  liquid  height. 
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f.  Viscosity  Measurement.  The  details  of  the  viscosity  measurement 
procedures  have  been  described  in  previous  Annual  Reports.  Tn  brief,  for  each 
fluid,  the  pressure  gage  is  disassembled,  cleaned  and  reassembled.  The  test 
fluid  is  degassed  before  being  put  into  the  pressure  gage.  The  pressure  gage 
is  kept  in  the  constant-temperature  bath  fcr  two  hours  to  ensure  thermal 
equilibrium.  Then,  the  system  is  pressurized  slowly.  The  test  fluid  is  raised 
from  the  reservoir  to  the  top  etched  line  of  the  viscometer  capillary  slcwly  by 
proper  manipulation  of  the  valve  settings. 

The  system  is  allowed  to  stand  for  another  hour  to  ensure  thermal  equilibrium. 
Then,  the  fluid  is  allowed  to  run  down  the  capillary  under  its  own  head.  The 
efflux  time  is  measured  with  a stop  watch  to  ± 0.1  second.  The  height  of  liquid 
from  reservoir  level  to  the  top  of  the  efflux  bulb  is  measured  by  cathetometer 
to  ± 0.005  cm.  The  reproducibility  of  the  pressure-viscosity  measurements  is 
shown  in  Table  2 in  terms  of  the  pressure  coefficient  which  is  defined  as: 


where : 


1 . Ip 

a = — In  — 

P 1* 


[22] 


J -1 


a = pressure  coefficient,  psig  , 

p = pressure,  psig, 

Pp  = viscosity  at  pressure  p,  cs. , and 

p0  = atmospheric  viscosity,  cs. 


g.  Effect  of  Variables  on  Pressure  Viscosity  Determinations.  Many  of 
the  variables  affecting  pressure  viscosity  determinations,  such  as  misalignment, 
drainage,  kinetic  energy,  surface  tension,  and  heat  of  compression  of  the  gas, 
have  been  studied  and  have  been  presented  in  detail  in  previous  reports  and  by 
Klaus,  et  al . (29).  By  means  of  proper  equipment  design  and  calibration 
procedures,  each  of  these  items  has  been  shown  to  have  negligible  effect  under 
normal  variation  on  the  pressure-viscosity  determination. 

The  effect  of  temperature  and  pressure  on  the  calibration  constant,  K, 
can  be  corrected  by  using  Equations  [20]  and  [21].  The  effect  of  dissolved  gas 
in  the  fluid  film  on  the  wall  of  the  viscometer  capillary  is  eliminated  by 
using  a dry  viscometer  capillary  for  every  run. 

Since  one  hour  in  the  constant  temperature  bath  is  required  to  eliminate 
the  effect  of  heat  of  compression,  a minimum  gas-liquid  contact  time  of  one 
hour  at  the  pressure  of  interest  is  required.  Data  on  Figure  5 show  the  effect 
of  gas-liquid  contact  time  on  the  pressure  coefficients. 

For  a close  approximation,  the  effect  of  nitrogen  gas  saturation  on  the 
pressure-viscosity  of  fluids  is  to  reduce  the  pressure  viscosity  to  its 
atmospheric  level.  On  this  basis  the  average  percent  saturation  of  di-2- 
ethylhexyl  sebacate  after  one  hour  at  both  3000  and  5000  psig  is  about  0.94. 

This  corresponds  to  0.29  and  0.59  percent  deviation  in  the  pressure  viscosity 
at  3000  and  6000  psig,  respectively.  For  the  polyolefin  oil  MLO  7754  after  one 
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hour  gas-oil  contact  time  at  3000  and  6000  psig,  the  percent  saturations  are  on 
the  order  of  0.35  and  3.52,  respectively.  This  corresponds  to  0.15  and  2.95 
percent  deviation  in  the  pressure  viscosity  at  3000  and  6000  psig,  respectively. 


i 


I * 


The  study  showed  that  the  gas  solubility  is  the  major  contribution  to  experi- 
mental error  for  pressures  as  low  as  3000  psig.  The  quantity  of  gas  dissolved  in 
the  fluid  increases  with  test  time  and  the  solubility  increases  with  system 
pressure.  Therefore,  PRL  Viscometer  should  be  operated  under  low  pressure  and 
with  minimum  contact  time  to  minimize  the  gas  solubility  problem.  The  limiting 
slope  of  tbe  viscosity-pressure  curve  at  0 psig  can  then  be  used  to  predict  the 
viscosities  at  higher  pressures. 
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Shear  stress  has  no  effect  on  the  viscosity  level  of  a Newtonian  fluid. 
For  most  non- Newtonian  fluids  the  shear  effect  is  avoided  by  operat  ng  the 
system  at  a shear  stress  below  100  dynes  per  square  centimeter.  It  has  been 
shown  that  for  shear  stresses  of  up  to  approximately  3 x 10'  dynes  per  square 
centimeter,  non-Newtonian  properties  of  the  types  of  polymer  solutions  used  in 
the  pressure  viscosity  study  are  not  observed  (33,  10). 


3.  Development  of  Variables.  Based  on  the  concept  that  the  -nobility  of 
the  molecules  in  a liquid  is  governed  by  their  free  volume,  different  models 
have  been  derived  to  describe  the  molecular  transport  phenomena  such  as  viscous 
flow  and  diffusion  The  basic  concept  underlying  free-volume  theories  implies 
that  each  molecule  in  a liquid  is  confined  to  a "cage"  bounded  by  its  immediate 
neighbors . 


Eyring  and  co-workers  (34)  interpreted  the  "cage"  as  an  energy  barrier. 
Viscous  flow  is  defined  in  this  system  as  the  result  of  the  movement  of  an 
activated  molecule  from  one  equilibrium  position  to  another  in  the  preferred 
direction  of  shear.  This  movement  is  treated  as  a jump  over  the  energy  barrier 
and  is  called  activation  energy  for  viscous  flow.  As  the  result  of  an 
elaborate  theoretical  treatment  of  this  model,  a rate  equation  of  the  following 
form  is  proposed. 
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E 4 
v is 

RT 


+ D 


[23] 


where: 


E 

vis 

n 

R 

T 

D 


activation  energy  for  viscous  flow, 

absolute  viscosity, 
molar  gas  constant, 
temperature,  and 
a constant. 


Both  E^g  and  D are  charac  ter  ist  ic  of  a given  liquid.  Evis  is  also  a 
function  of  pressure. 

Cohen  and  Turnbull  (1,35)  Introduced  the  concept  that  the  statistical 
redistribution  of  the  free  volume  occasionally  opens  up  a void  within  the  "cage" 
sufficiently  large  to  permit  considerable  displacement  of  the  molecule 
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contained  by  it.  For  the  molecular  displacement  to  be  diffusive,  the  void 
created  should  exceed  a characteristic  minimum  volume  which  approaches  the 
molecular  volume,  and  should  be  just  large  enough  to  enable  another  molecule 
to  jump  in  after  this  displacement. 

By  treating  the  free  volume  of  the  individual  molecules  of  a liquid  as  a 
random  variable,  and  by  describing  the  distribution  ot  the  free  volume  among 
the  numerous  molecules  by  a simple  exponential  distribution  function,  Cohen  and 
Turnbull  obtained  the  following  relationship. 


In  n 


+ In  He 


where : 


vf 

v*/vf 


absolute  viscosity, 
a constant, 

overlap  factor  between  1/2  and  1, 
minimum  required  specific  free  volume, 
specific  free  volume,  and 
is  a function  of  temper,  cure  and  pressure. 


[24] 


Following  this  concept  Roeland  (28)  used  the  Weibull  distribution  function 
instead  of  the  simple  exponential  distribution  function  to  approximate  the 
distribution  of  free  volume  and  developed  a new  correlation: 


In  n * 
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[25] 


where: 


n 


^e 

vf 

Vf* 


n 


Vf*/vf 


absolute  viscosity, 
fictitious  viscosit''  for  Vf  = 
specific  free  volume, 

minimum  required  specific  free  volume, 

shape  parameter  which  is  essentially  independent  of 
temperature,  and  which  predominately  determines  the 
overall  form  of  the  statistical  distribution  of  the 
free  volume  of  liquid  concerned  as  a function  of 
Vf*/Vf,  and 

is  a function  of  temperature  and  pressure. 


a.  Atmospheric  Viscosity.  From  the  concept  that  the  mobility  of  the 
molecule  in  a liquid  is  governed  by  its  free  volume,  kinematic  viscosity  can  be 
considered  as  a function  of  three  variables:  (1)  free  volume  which  is  defined 

as  the  volume  per  cage  available  for  free  redistribution,  (2)  flow  unit  which 
is  considered  to  be  the  volume  completely  occupied  by  an  individual  molecule 


or  ar.  effective  segment  of  a molecule,  and  (3)  average  formula  weight  per  atom 
which  is  molecular  weight  divided  by  the  total  number  of  atoms  in  the  molecule. 
Then,  at  atmosphere  pressure  and  at  the  temperature  of  interest, 


whe  r e : 


Uo  T ' ' lV  <*•  Vl.P  1 

o 

atmospheric  kinematic  viscosity, 
average  formula  weight  per  atom, 
flow  unit, 
free  volume, 

atmospheric  pressure,  and 
temperature  of  interest. 


b.  Atmospheric  Density.  The  differences  in  density  among  isomers  at 
constant  temperature  and  pressure  are  due  to  the  differences  in  packing.  It  is 
logical  to  assume  that  tighter  packing  corresponds  to  smaller  free  volume.  For 
liquids  with  the  same  average  formula  weight  per  atom,  the  differences  in 
density  at  constant  temperature  and  pressure  are  due  to  the  differences  in 
packing.  Density  can  be  used  to  measure  the  difference  of  free  volume  on  a 
relative  basis.  High  density  corresponds  to  small  free  volume.  Low  density 
corresponds  to  a large  free  volume.  At  atmospheric  pressure  and  the  temperature 
of  interest: 

T * f '"a-  (Vt'  1271 

where : 

p0  = atmospheric  density, 

= average  formula  weight  per  atom, 

V £ = free  volume,  and 

T = temperature  of  interest. 

c.  AS  I'M  Slope.  In  assessing  the  effect  of  temperature  on  the  ratio 
Vj*/Vf  in  Equation  [25|,  Roeland  (19)  assumed  the  following  expression: 


where: 


vf  (T-T) r 

= minimum  required  specific  free  volume, 
= specific  free  volume,  and 
= absolute  temperature. 


The  exponent  r approximates  unity  and  the  parameter  T0  denotes  the  absolute 
temperature  where  the  specific  free  volume,  Vf,  of  the  liquid  would  become  zero. 
Further,  the  constant  k represents  the  value  of  the  ratio  Vf*/Vf  that  would  be 
reached  when  the  absolute  temperature  T is  decreased  to  1 “K  above  the  limiting 
temperature  Tn. 


When  Equation  (28]  is  substituted  in  Equation  [25],  the  rearranged 
Equation  [29]  results: 

In  In  [— ] = - rn  in  (T-T  ) ■+•  In  [kr  (1  + — )] 

h on 


Assume  T0  = 0 


where : 


~ n + c 
e 


In  In  (n  -*■  c)  = - rn  In  (T)  + In  [kr  (1  + ~))r' 


r|  = absolute  viscosity, 

rie  = absolute  viscosity  for 
equal  infinite 

r = constant  approximately 
equal  to  unity, 

n * shape  parameter, 

T = absolute  temperature, 

k = Vf*/Vf  when  T equal  1 °K, 

c = a constant, 

Vf  » specific  free  volume,  and 

Vj*  = minimum  required  specific 
free  volume. 


MacCoull  - Walther  (36)  developed  an  emperical  viscosity-temperature 
equation  which  has  the  following  form: 


log  log  (po  + c)  = 


log  T + N 


where : 


pc  = kinematic  viscosity,  cs, 

T = absolute  temperature,  °R, 
mQ,N0  = empirical  constants,  and 
c = 0.6  for  y0  > 1.5  cs . 

0.65  for  1.0  < p0  < 1.5  cs. 

0.70  for  0.7  < p0  < 1.0  cs. 

0.75  for  0.4  < pQ  < 0.7  cs. 

On  the  basis  of  the  McCoull  - Walther  equation,  standardized  ASTH  kinematic 
viscosity-temperature  charts  (37)  in  which  the  ordinate  is  a log  log  scale  of 
viscosity  (cs.)  and  the  abscissa  is  a log  temperature  (“F)  scale,  were  developed. 
ASTM  slope  of  a liquid  may  be  defined  as  the  negative  "measured"  slope,  or  the 
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measured  linear  distance  on  the  ordinate  of  the  ASTM  charts  C and  D (D  341-39) 
between  the  viscosity  values  at  100°  and  210°F  and  divided  by  the  measured 
linear  distance  between  210°  and  100°F  on  the  abscissa.  Any  other  two 
temperatures  may  be  used  with  the  designation  of  ASTM  slope  between  two  specified 
temperatures  as  shown  in  Figure  6.  Therefore,  the  ASTM  slope  is  directly 
proportional  to  the  slope  mQ  in  Equation  [31].  For  kinematic  viscosities  above 
1.5  cs. , slope  mQ  is  equal  to  4.5  times  the  ASTM  slope.  Proportionality  constants 
for  kinematic  viscosities  below  1.5  cs.  are  listed  in  Table  3.  ASTM  slope  is 
also  tabulated  in  table  form  (38). 


A comparison  of  Equations  [30]  and  [31]  shows  that  since  c can  be  assumed 
to  be  a constant  for  kinematic  viscosities  above  1.5  cs.,  Mc  is  directly 
proportional  to  rn.  Since  r is  close  to  unity,  Mq  is  approximately  directly 
proportional  to  n.  This  means  that  ASTM  slope  is  also  approximately  directly 
proportional  to  the  shape  parameter  n in  Weibull's  function  for  describing  the 
statisti  al  distribution  of  free  v lume  amongst  the  numerous  molecules  of  a 
given  liquid.  Conversely,  the  ASTM  slope  may  now  be  interpreted  as  the  quantity 
which  predominantly  determinates  the  overall  form  of  the  statistical  distri- 
bution of  the  free  volume  of  the  liquid  concerned  as  a function  of  temperature. 
This  in  turn  means  that  ASTM  slope  can  be  considered  as  a measurement  of  the 
flexibility  of  the  molecule.  A rigid  molecule  will  have  a large  ASTM  slope, 
while  a flexible  molecule  will  have  a low  ASTM  slope. 


d.  Relating  the  Pressure  Coefficient  to  Atmospheric  Viscosity,  Atmo 
pheric  Density  and  ASTM  Slope.  As  long  ago  as  1893  Barus  (39)  established  an 
empirical  equation  to  describe  the  isothermal  viscosity-pressure  relationship 
for  a given  liquid. 


where 


viscosity  at  pressure  p, 

atmospheric  viscosity  at  the  temperature  of  interest,  and 

pressure  coefficient  which  for  a given  liquid  is  a function 
of  temperature  only. 


An  equation  of  the  same  form  was  later  obtained  as  a theoretical  approxi- 
mation by  the  Eyring  theory  of  viscous  flow.  The  isothermal  viscosity-pressure 
relationship  for  a given  liquid  can  be  expressed  as: 


where 


viscosity  at  pressure  p, 

atmospheric  viscosity  at  temperature  of  interest  T 
characteristic  constant  for  a given  liquid,  and 
the  molar  gas  constant. 


According  to  the  Eyring  Model  of  Viscous  Flow,  the  indicated  jump  of 
activated  molecules  over  a potential  energy  barrier  is  only  possible  if  a suit- 
able "hole"  has  been  created  in  the  neighborhood.  The  size  of  such  holes  (per 
mole)  is  represented  by  the  quantity  Vvis  • 

Vvis  should  increase  with  flow  unit  size  because  a larger  "hole"  is  required 
for  a larger  flow  unit.  A more  rigid  molecule  would  require  a larger  "hole"  than 
a flexible  one.  By  introducing  the  concept  of  molecular  interlocking,  the  increase 
of  Vvis  with  respect  to  ASTM  slope  can  be  further  explained.  An  increase  in  ASTM 
slope  corresponds  to  an  increase  in  molecular  rigidity.  Molecular  interlocking 
is  more  effective  for  more  rigid  molecules,  which  in  turn  increases  the  flow  unit 
and  the  W The  effect  of  free  volume  on  molecular  interlocking  is  more 
abstract.  When  the  free  volume  is  large,  the  amount  of  molecular  interlocking 
introduced  by  applying  external  pressure  is  larger  than  when  the  free  volume  is 
small.  In  other  words,  when  the  free  volume  is  small,  the  molecules  are  already 
tightly  packed.  Molecular  interlocking  already  exists.  Applying  pressure  just 
forces  the  already  interlocked  molecules  to  interlock  tighter.  Thus  the  effect 
on  the  flow  unit  and  Vvis  is  much  smaller  than  when  the  free  volume  is  large 
and  external  pressure  introduces  interlocking  of  free  molecules. 

The  same  type  of  logic  can  also  apply  to  the  statistical  distribution  of  the 
free  volume  concept  introduced  by  Cohen  and  Turnbull  (1,  35).  From  Equation  [32] 
the  pressure  coefficient  can  be  interpreted  as  the  rate  of  log  viscosity  change 
with  respect  to  pressure  or  simply  the  pressure  effect  on  the  opportunity  to 
locate  a proper  size  free  volume.  The  pressure  effect  on  the  free  volume 
decreases  the  opportunity  for  a large  flow  unit  to  find  the  right  size  free 
volume  over  that  of  a small  flow  unit.  The  reduction  of  free  volume  by  external 
pressure  is  also  more  significant  in  decreasing  the  opportunity  for  a rigid 
molecule  to  locate  a right  size  free  volume  than  for  a more  flexible  molecule. 
Again,  interlocking  is  more  effective  for  rigid  molecules,  which  in  turn 
increases  the  size  of  the  flow  unit  and  makes  it  more  susceptible  to  pressure 
effects.  Moreover,  the  amount  of  molecular  interlocking  introduced  by  external 
pressure  is  larger  for  large  free  volume  than  for  smaller  free  volume  fluids. 

As  a result,  the  pressure  effect  on  the  opportunity  to  find  a right  size  free 
volume  is  larger  for  a cluid  which  has  large  free  volume  than  for  a fluid  which 
has  small  free  volume. 

From  the  above  discussion,  it  is  clear  that  the  pressure  coefficient,  a,  is 
a function  of  atmospheric  viscosity,  atmospheric  density  and  ASTM  slope. 

4.  Correlation  Methods  and  Results.  In  order  to  study  the  pressure  effect 
on  the  viscosities  of  Pennsylvania  paraffinic  mineral  oils,  mineral  oil  blends, 
and  a resin  blend,  pressure-viscosities  of  a series  of  paraffinic  mineral  oils 
ranging  from  gas  oils  to  bright  stock,  two  mineral  oil  blends  and  a Pa.  resin 
blend  were  evaluated  at  pressures  up  to  3000  psig  at  100=F.  For  some  of  the 
fluids  the  pressure-viscosities  were  also  measured  at  32°  and  130°F.  As  shown  in 
Figures  7 through  13,  plots  of  the  log  of  kinematic  viscosity  versus  pressure 
are  straight  lines  for  pressures  up  to  3000  psig.  This  is  in  agreement  with  the 
results  obtained  in  previous  studies  (AFML-TR-67-107 , Part  I;  AFML-TR-70-304, 

Part  II).  These  data  also  showed  that  almost  all  fluids  when  plotted  in  terms 
of  log  kinematic  viscosity  versus  pressure  result  in  straight  lines  for  pressures 
up  to  10,000  psig.  Thus,  the  pressure  oefficient  a,  is  independent  of  pressure 
when  it  is  expressed  as  follows: 


► 
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where : 


log  -p  = a p log  (e)  or  In  — = a p [34] 

Uo 

4 = kinematic  viscosity  at  pressure  p,  centistokes, 

U0  = atmospheric  kinematic  viscosity  at  the  temperature 
of  interest,  centistokes, 

a = pressure  coefficient,  1/psig, 

p = pressure,  psig, 

e = 2.718, 

log  = base  10  logarithm,  and 
In  = nature  logarithm 


Data  for  the  mineral  oils,  mineral  oil  blends  and  Pa.  resin  blend  at  100°F 
produced  by  this  investigation  were  used  as  part  of  the  data  to  establish  the 
new  correlations. 

Data  from  previous  studies  at  this  laboratory  and  from  an  API  Project  42 
were  also  used  to  develop  a correlation  at  100°F.  The  data  for  this  correlation 
was  further  restricted  to  an  average  formula  weight  in  the  neighborhood  of  4.5, 
and  a viscosity  level  above  one  centistoke.  This  eliminated  some  of  the  non- 
hydrocarbon data  from  AFML-TR-67-107,  Part  I,  and  a few  fluids  from  the  API 
Project  42  (8).  All  of  the  polymer  blends  and  resin  blends  discussed  in  AFML-TR- 
70-304,  Part  II  were  included  to  develop  these  correlations.  Data  from  the  ASME 
Report  (7)  and  Roelands  (23)  were  not  included  because  of  the  low  level  of 
precision  in  the  low  pressure  region.  Brief  descriptions  and  physical  properties 
of  fluids  used  to  develop  the  correlations  are  listed  in  Tables  4 through  9. 

These  data  cover  a density  range  of  from  0.75  to  0.96  gm/cc,  and  a range  of  ASTM 
slope  of  from  0.28  to  0.93,  and  a viscosity  range  of  from  1.4  to  1984  centistokes. 

a.  Multiple  Linear  Regression  Analysis.  In  section  three,  development 
of  variables,  it  was  shown  that  the  pressure  coefficient  is  probably  a function 
of  atmospheric  density,  atmospheric  viscosity,  and  ASTM  slope.  However,  the 
functional  form  of  the  relationship  was  not  known.  It  was  decided  that  multiple 
linear  regression  analysis  would  be  used  as  a first  step  to  search  for  the 
functional  form.  Multiple  regression  analysis  is  defined  as  linear  when  in 
searching  for  the  best  fit,  only  the  constant  and  the  coefficient  of  the  variable 
are  being  changed.  The  best  fit  is  obtained  by  the  least  squares  method. 

Among  the  many  available  multiple  linear  regression  methods,  the  step-up 
multiple  linear  regression  of  the  statistical  package  program  from  the  Computation 
Center  of  The  Pennsylvania  State  University  was  used.  This  method  allows  the 
user  to  examine  the  significance  of  each  individual  variable  used  in  the  corre- 
lation. 

In  step-up  multiple  linear  regression,  intermediate  results  are  used  to  give 
statistical  information  at  each  step  in  the  calculation.  These  intermediate 
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answers  are  alro  used  to  control  the  method  of  calculation.  A number  of  Inter- 
mediate regression  equations  are  obtained  as  well  as  the  final  multiple  regression 
equation.  These  equations  are  obtained  by  adding  or  deleting  one  variable  at  a 
time  as  follows: 

y = b(o)  + b(l)  x (1)  [35] 

y = b'(o)  + b'(l)  x (1)  + b ' ( 2)  x (2)  [36] 

y = b"(o)  + b'(l)  x (1)  + b" (2)  x (2)  + b"(3)  x (3)  [37] 

The  variable  added  is  the  one  which  makes  the  greatest  improvement  in  "goodness 
of  fit."  The  coefficients  represent  the  best  values  when  the  equation  is  fitted 
using  the  specific  variables  included  in  the  equation.  An  important  feature  of 
this  procedure  is  that  a variable  may  be  indicated  to  be  significant  in  an  early 
stage  and  enters  the  equation,  but  after  several  variables  are  added  to  the 
equation,  the  initial  variable  may  be  indicated  to  be  insignificant.  The 
insignificant  variable  is  then  removed  from  the  regression  equation  before 
another  variable  is  added.  Therefore,  only  significant  variables  are  included 
in  the  final  multiple  regression  Equation  [44], 

The  three  variables,  atmospheric  viscosity,  atmospheric  density  and  ASTM 
slope,  were  evaluated  as  independent  variables  and  associated  variables.  A set 
of  pseudo  linear  variables  to  a maximum  of  fourth  order  was  created.  The  maximum 
of  fourth  order  was  used  because  it  is  desirable  to  give  the  regression  the 
maximum  freedom  but  also  to  minimize  the  possibility  of  obtaining  a correlation 
with  abrupt  fluctuations  which  are  very  unlikely  in  the  actual  physical  situation. 
The  variable,  atmospheric  viscosity  was  expressed  in  terms  of  log  atmospheric 
viscosity  to  reduce  the  range  of  its  variation.  The  pseudo  linear  variables 
are  shown  in  Table  10. 

In  curve  fitting,  as  the  number  of  variables  increases,  the  average 
deviation  (see  Table  11)  always  decreases  until  it  becomes  zero.  There  is  no 
unique  statistical  procedure  to  determine  the  best  regression  equation.  One  of 
the  most  common  indicators  is  standard  error  of  estimate  (see  Table  11).  This 
indicator  should  always  decrease  with  the  addition  of  more  variables.  If  the 
standard  error  of  estimate  increases  with  the  additional  variable,  this  indicates 
that  the  addition  of  the  variable  isn't  statistically  sound.  However,  the  over- 
riding factor  is  always  whether  the  model  described  by  the  regression  equation 
is  physically  reasonable.  Moreover,  in  order  for  the  regression  equation  to  be 
useful,  it  has  to  be  convenient  and  economical  to  use.  Figure  14  shows  the 
standard  error  of  estimate  versus  number  of  variables  from  the  multiple  linear 
regression  analysis.  The  regression  equations  with  one  to  four  variables  are 
shown  in  Table  12.  After  a detailed  examination  of  plots  of  the  regression 
equations  obtained,  it  was  found  that  none  of  the  regression  equations  was 
satisfactory.  When  the  number  of  variables  was  kept  low,  the  average  of 
deviation  (see  Table  11)  was  deemed  to  be  too  large.  When  the  number  of 
variables  was  increased,  the  model  obtained  was  considered  to  be  physically 
unreasonable.  It  was  then  decided  to  use  multiple  nonlinear  regression  analysis. 

b.  Multiple  Nonlinear  Regression  Analysis.  When  the  regression  analysis 
involves  multiple  variables  with  nonlinear  parameters,  it  is  called  multiple  non- 
linear regression  analysis.  The  type  of  nonlinear  model  being  considered  here  is 
the  polynominal  type  with  variable  powers.  One  of  the  advantages  of  multiple 
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nonlinear  regression  analysis  is  that  it  provides  maximum  orderly  flexibility 
with  a minimum  number  of  variables.  However,  multiple  nonlinear  regression 
analysis  requires  definite  functional  form  and  good  initial  values  for  the 
parameters.  If  the  initial  values  for  the  parameters  are  too  far  off,  the  system 
will  probably  diverge. 

Equations  [39],  [40],  and  [41]  from  the  multiple  linear  regression  analysis 
were  picked  as  the  functional  form.  The  multiple  nonlinear  regressional  analysis 
program,  NLIN-2,  from  the  Computation  Center  of  The  Pennsylvania  State  University, 
University  Park,  Pennsylvania,  was  used.  The  program  is  based  on  least  squares 
estimation.  Using  parameter  values  from  the  multiple  linear  regression  analysis 
as  initial  values,  all  three  of  the  equations  successfully  converged.  The 
regression  equations  are  shown  in  Table  13  as  Equations  [42],  [43],  and  [44] 
respectively.  Results  show  that  the  ASTM  slope  term  in  the  third  variable  of 
Equation  [43]  is  insignificant.  Therefore,  a new  function  was  created  by  omitting 
the  ASTM  slope  term.  The  new  regression  equation  is  shown  in  Table  13  as 
Equation  [45].  After  examining  Equations  [42],  [44],  and  [45]  in  detail  both 
from  the  physical  soundness  and  the  statistical  (standard  error  of  estimate) 
point  of  view.  Equation  [451  is  the  best  among  all  the  regression  equations 
tested. 

Since  only  ASTM  slope  and  atmospheric  viscosity  were  used  as  independent 
variables  in  Equation  [42],  it  was  assumed  that  the  differences  in  free  volumes 
were  being  accounted  for  by  the  use  of  kinematic  viscosity.  Equation  [42]  can 
be  used  as  a close  approximation  when  atmospheric  viscosity  data  are  not 
available.  As  it  will  be  shown,  Equation  [42]  is  the  better  equation  when  the 
average  formula  weights  per  atom  of  the  fluids  deviate  significantly  from  the 
4.5  average.  In  order  to  be  sure  that  an  additional  term  in  Equation  [42]  would 
not  improve  the  equation  significantly,  a new  group  of  pseudo  linear  variables 
made  up  of  the  independent  variables,  ASTM  slope  and  atmospheric  viscosity,  to 
the  maximum  of  fourth  order  were  constructed.  The  multiple  linear  regression 
analysis  was  again  used  to  pick  the  third  variable.  The  new  functional  form  was 
used  in  the  multiple  nonlinear  regression  analysis  to  obtain  Equation  [46],  As 
shown  in  Table  13,  Equation  [46]  has  a larger  standard  error  of  estimate.  In 
terms  of  physical  reasonableness.  Equation  [46]  is  also  inferior  to  Equation  [42]. 

c.  Discussion  and  Comparison  of  Equations  [45]  and  [42].  Data 
calculated  by  the  relationship  expressed  in  Equation  [45]  are  plotted  on  a semi- 
log scale  of  pressure  coefficient  versus  log  of  the  atmospheric  viscosity  with 
lines  of  constant  ASTM  slope.  The  fluids  of  constant  density  levels  of  0.75, 

0.85,  0.95  and  1.00  gm/cc  are  shown  on  Figures  15,  16,  17  and  18,  respectively. 

The  figures  show  that  all  constant  ASTM  slope  lines  converge  to  a single  point 
at  an  atmospheric  viscosity  level  of  one  centistoke.  This  is  due  to  the 
constraint  of  Equation  [45].  The  significant  message  from  Figures  15  through  18 
is  that  the  dependency  of  the  pressure  coefficient  on  ASTM  slope  decreases  with 
decreasing  atmospheric  viscosity  level.  This  result  is  shown  more  clearly  in 
Figure  19.  This  trend  can  be  explained  in  terms  of  molecular  interlocking.  For 
constant  free  volume,  low  viscosity  level  represents  a comparatively  small  flow 
unit  which  has  relatively  little  branching  and  coiling.  The  rigidity  of  the 
molecule,  ASTM  slope,  has  little  effect  on  the  amount  of  molecular  interlocking 
introduced  by  the  external  pressure.  Therefore,  the  effect  of  ASTM  slope  on  the 
pressure  coefficient  is  small.  For  low  ASTM  slopes,  the  molecules  are  relatively 
flexible,  and  the  molecular  interlocking  introduced  by  the  external  pressure  is 
not  very  effective.  These  effects  explain  the  lower  effects  of  fluids  with  low 
ASTM  slope  and/or  low  atmospheric  viscosity  on  the  pressure  coefficient. 
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The  increase  of  pressure  coefficient  with  respect  to  atmospheric  viscosity 
was  explained  in  detail  in  Section  3 on  the  development  of  variables.  Data  from 
Figures  15  through  18  show  a leveling  off  of  the  pressure  coefficient  at  high 
atmospheric  viscosity.  This  can  be  explained  by  segmental  flow.  Long  chain 
molecules  have  a high  degree  of  freedom.  The  flow  unit  is  no  longer  a function 
of  molecular  chain  length.  The  flow  unit  is  only  a segment  of  the  long-chain 
molecule.  The  concept  of  segmental  flow  has  been  suggested  by  many  authors  (41). 
For  some  density  levels  the  use  of  Equation  [45]  actually  predicts  a decrease  in 
pressure  coefficient  with  increasing  atmospheric  viscosity.  This  is  considered 
to  be  physically  unreasonable.  As  a result,  horizontal  dotted  lines  were  drawn 
in  the  regions  where  Equation  [45]  predicts  a decrease  in  pressure  coefficient. 
This  analysis  is  in  compliance  with  the  segmental  flow  concept.  This  functional 
form  was  chosen  by  the  statistical  analysis  based  on  the  data  available  where 
pressure  coefficient  appears  to  decrease  with  viscosity.  This  apparent  decrease, 
upon  analysis,  appears  to  represent  no  change  in  viscosity-pres3ure  coefficient 
with  viscosity  level  within  the  limits  of  error  of  the  measurement  of  the 
properties  measured.  Some  of  the  errors  in  measurement  come  from  the  data 
obtained  from  different  investigators. 

Considering  the  nature  of  the  data  obtained  from  Equation  [45],  Figure  20 
was  constructed  to  show  the  operational  regions  of  the  equation.  This  region 
covered  all  of  the  data  used  for  the'correlation.  The  maximum  deviation  of 
0.111  x 10  4 psig  from  the  horizontal  line  is  well  within  the  experimental  limits 
of  different  investigators.  There  are  four  data  points  inside  this  region.  The 
differences  between  the  predicted  values  from  Equation  [45]  and  from  Figures  15 
through  18  are  shown  in  Table  14.  The  results  show  that  the  differences  are 
insignificant. 

The  effect  of  atmospheric  density  at  constant  atmospheric  viscosity  and 
ASTM  slope  on  the  pressure  coefficient  based  on  Equation  [45]  is  shown  in 
Figure  21.  The  decrease  in  the  pressure  coefficient  with  increasing  atmospheric 
density  has  been  explained  in  detail  in  Section  3.  This  effect  is  due  to  the 
different  effects  on  the  size  of  the  flow  unit  introduced  by  the  newly  formed 
interlocking  molecules  and  the  tightening  of  the  already  interlocked  molecules. 
These  molecular  configurations  in  turn  affect  the  pressure  coefficient.  The 
data  on  Figure  21  show  that  the  effect  of  atmospheric  density  is  much  greater 
for  high  viscosities  than  for  low  viscosities.  At  high  viscosity  levels  there 
is  more  branching  and  coiling  on  the  molecules.  The  free  volume  has  a strong 
effect  on  the  molecular  interlocking  introduced  by  the  application  of  external 
pressure. 

Data  derived  from  the  relationship  in  Equation  142]  are  plotted  on  Figure  22 
as  pressure  coefficient  versus  the  log  of  the  atmospheric  viscosity  with  lines  of 
constant  ASTM  slope.  The  general  trend  of  these  data  is  the  same  as  those  in 
Figures  15  through  18.  The  biggest  difference  between  the  two  equations  is  that 
no  segmental  flow  is  indicated  from  the  data  on  Figure  22.  The  density  term  in 
the  kinematic  viscosity  has  less  effect  than  the  differences  in  free  volume. 

By  not  accounting  for  the  differences  of  free  volumes  between  fluids  at  high 
viscosity,  the  segmental  flow  effect  is  obscured. 

In  order  to  understand  the  significance  of  the  deviations  between  the  experi- 
mental values  and  the  predicted  values,  the  discrepency  of  the  experimental  data 
among  different  investigators  must  be  considered.  The  experimental  values  of  the 
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pressure  coefficients  for  the  same  fluids  as  obtained  by  different 
investigators  are  compared  on  Table  15.  In  order  to  show  the  effect  of  the 
variation  of  atmospheric  viscosity  and  density,  and  ASTM  slope  on  the  prediction 
of  the  pressure  coefficient,  the  deviations  of  the  predicted  values  among 
different  investigators  based  on  Equation  [45]  are  also  shown  in  Table  15.  From 
the  data  on  Table  15,  the  maximum  deviation  among  reported  experimental  values  on 


pressure  coefficient  is  0.154  x 10  4 psig  1 and  the  maximum  deviation  among  the 
predicted  values  based  on  Equation  [45]  is  0.098  x 10-4  psig-1.  Therefore,  the 
combined  maximum  deviation  can  be  as  large  as  0.252  x 10  4 psig. 


The  data  on  Figure  21  show  that  atmospheric  density  has  little  effect  on 
the  pressure  coefficient  at  low  viscosity.  As  the  viscosity  increases,  the 
effect  of  atmospheric  density  on  the  pressure  coefficient  increases  substantially. 
This  indicates  that  the  differences  between  Equation  [45]  and  [42]  should  be  small 
at  low  viscosities.  The  differences  should  increase  with  increasing  atmospheric 
viscosity.  The  comparisons  of  Equations  [45]  and  [42]  in  predicting  pressure 
coefficients  used  to  develop  the  equation  at  atmospheric  viscosity  range  of  1 to 
10  cs.,  1C  to  100  cs.  and  100  to  2000  cs.  are  listed  in  Tables  16,  17,  and  18. 

A summary  of  the  comparisons  is  shown  in  Table  19.  The  summary  indicates  that 
the  differences  between  Equations  [45]  and  [42]  are  mainly  in  the  atmospheric 
viscosity  range  of  100  to  2000  centistokes. 


d.  Temperature  Effects.  Combining  Frying' s theory  on  viscous  flow 
(Equation  [33])  and  Equation  [34] 


In  -2- 


Cr  ' v vis 


= In  f - 

lull  a 


i n I . , o i 

In  — + In  — — 


P I I a 


where 


absolute  viscosity  at  p,  cp, 

atmospheric  viscosity  at  temperature  of  interest,  cp, 
molar  gas  constant, 
temperature  of  interest, 

characteristic  constant  for  a given  liquid,  cm^, 
kinematic  viscosity  at  p,  cs, 

atmospheric  kinematic  viscosity  at  temperature  of  interest,  cs, 
density  at  p,  gm/cc. 


atmospheric  density  at  temperature  of  interest,  gm/cc, 
pressure  coefficient,  psig-^,  and 
natural  logarithm. 


Since  In  (pG/p)  is  on  the  order  of  2.5  percent  of  In  (U/M0),  It  can  be  assumed 
that 
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Since  V . 
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can  be  considered  independent  of  temperature. 


where: 


100°F 

T 


560 

aT  ai003F  Jt  + 460_ 

pressure  coefficient  at  temperature  T,  psig  ^ , 
pressure  coefficient  at  100°F,  psig-J-,  and 
temperature  of  interest,  °F. 


[49] 


The  pressure  coefficient  of  a fluid  at  100°F  is  controlled  by  atmospheric 
viscosity  and  density  at  100°F  along  with  the  ASTM  slope  which  can  be  considered 
to  be  independent  of  temperature.  The  pressure  coefficient  of  a fluid  at  the 
temperature  of  interest  should  then  be  controlled  by  the  atmospheric  viscosity 
and  density  at  the  temperature  of  interest  along  with  the  ASTM  slope.  A 
comparison  of  equations  [49],  [45]  and  [42]  at  several  temperatures  is  shown  in 
Tables  20  through  23.  A summary  of  the  comparisons  is  shown  in  Table  24. 
Differences  in  the  average  absolute  deviation  is  defined  as  the  average  absolute 
deviation  at  a given  temperature  minus  the  average  absolute  deviation  at  100°F. 

The  results  show  that  Equation  147],  which  was  developed  from  Erying's  theory  on 
viscous  flow,  can  only  be  used  as  a rough  estimation  of  the  effects  of  temperature 
on  the  pressure  coefficients.  Equation  [42]  with  ASTM  slope,  and  atmospheric 
viscosity  at  the  temperature  of  interest  is  a second  choice  if  the  information 
for  atmospheric  density  at  the  temperature  of  interest  is  not  available.  In 
general,  the  average  absolute  deviations  of  the  fluids  at  all  the  temperatures 
of  interest  other  than  100°F  are  about  the  same  as  the  average  absolute 
deviations  of  the  fluids  at  100°F.  The  exception  is  in  the  fluids  at  32'F. 

This  is  at  least  partially  due  to  the  fact  that  at  such  a low  temperature,  the 
ASTM  slopes  of  some  polymeric  fluids  and  dense  or  congested  center  molecules  do 
tend  to  change  from  their  100°-210l,F  ASTM  slope  values  (42,  43). 
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5.  Comparison  of  Equations  [45]  and  [42]  with  other  Methods  of  Prediction. 
According  to  the  analysis  of  the  Technical  Data  Book  - Petroleum  Refining  (44) 
and  previous  studies  at  this  laboratory,  the  correlations  provided  by  Kouzel  (26), 
Roelands  et  al.  (23),  and  the  data  reported  in  previous  Annual  Reports  are  the 
best  available  to  predict  the  effect  of  pressure  on  viscosity.  The  detailed 
discussions  of  the  above  correlations  are  in  Section  1 of  this  study.  A summary 
of  the  various  correlations  used  for  the  comparison  is  shown  in  Table  25.  Lists 
of  physical  properties  required  for  each  correlation  are  shown  on  Table  27. 
Pressure  viscosity  values  at  5000  psig  were  used  to  calculate  the  pressure 
coefficient  for  the  correlations  requiring  pressure.  Since  some  of  the  corre- 
lations were  based  on  absolute  viscosity,  the  ratio  of  density  at  5000  psig,  p, 
and  atmospheric  density,  p0,  was  required.  In  general  p/p0  only  varies  from 
0.9815  to  0,9840.  The  exact  value  can  be  estimated  from  available  data  and 
correlations  (31). 


Pressure  coefficient  data  from  four  different  sources  were  used  for 
comparison.  The  comparison  of  various  methods  for  predicting  the  pressure 
coefficient  at  100CF  using  the  data  from  the  current  investigation  and  that 
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taken  from  Annual  Report  AFML-TR-67-107,  Part  I,  shown  in  Table  28.  The  fluid 
types  include  mineral  oils,  mineral  oil  blends,  resin  blends  and  nonhydrocarbons. 
Brief  descriptions  of  the  fluids  are  shown  on  Tables  4 and  5.  Polymer  and  resin 
blends  at  100JF  are  compared  in  Table  29.  The  comparisons  or  the  PRL  pressure 
coefficient  data  at  32°  and  130°F  are  shown  on  Table  30.  The  above  comparisons 
are  summarized  in  Table  31  in  terms  of  average  absolute  deviation,  bias  and  RMS. 
The  definitions  of  these  terms  are  given  in  Table  11.  Data  on  Table  31  show 
that  the  predicted  pressure  coefficient  from  Fresco's  correlation  in  general 
tends  to  be  lower  than  the  experimental  values.  Roelands'  correlation  for 
polymer  blends  showed  a tendency  to  predict  a higher  pressure  coefficient.  The 
large  average  absolute  deviation  in  Roeland's  prediction  for  data  from  Fresco  at 
100°F  is  largely  due  to  its  inability  to  predict  the  pressure  coefficient  of 
hydrogenated  polybutene.  Kouzel's  correlation  predicts  pressure  coefficients 
consistently  below  the  experimental  values. 

The  comparison  of  various  methods  of  predicting  the  pressure  coefficient  of 
hydrocarbons  from  the  API  Project  42  (8)  is  shown  on  Tables  32  and  33.  The  data 
is  summarized  in  Table  34.  The  pressure  coefficient  data  for  pure  hydrocarbons, 
mineral  oils,  and  polymer  blends  from  the  ASME  Report  (7),  and  data  on  mineral 
oils  from  Roelands  et  al.  (23)  were  used  for  comparisons  in  Tables  36  and  37, 
respectively.  The  summaries  of  the  comparisons  are  shown  in  Tables  38  and  39. 

In  comparing  the  bias  between  the  different  data  sources,  it  is  interesting  to 
note  that  the  data  from  Roelands  et  al.  gives  larger  positive  bias  which 
indicates  that  the  experimental  data  from  Roelands  et  al.  tends  to  have  higher 
experimental  pressure  coefficients. 

In  order  to  understand  the  effectiveness  of  the  various  methods  in 
predicting  the  pressure  coefficients  of  different  fluid  types,  the  data  were 
separated  into  a mineral  oil  group,  a pure  hydrocarbon  group,  a resin  blend  and 
polymer  blends  and  a nonhydrocarbon  group.  The  resin  and  polymer  blends  were 
considered  as  one  group  since  the  data  reported  in  AFML-TR-70-304,  Part  II, 
showed  that  resin  and  polymer  blends  can  be  correlated  on  the  same  basis.  The 
summary  of  the  comparison  is  shown  in  Table  40. 

In  the  course  of  examining  the  difference  in  the  average  absolute  deviations 
based  on  different  methods,  it  should  be  noted  that  the  correlations  do  not 
cover  the  same  range  of  fluids.  The  method  of  Fresco  does  not  include  polymer 
blends  and  resin  blends.  On  the  other  hand,  the  data  from  AFM1 -TR-70-304 ; Part  II 
are  useful  only  for  polymer  and  resin  blends.  Roelands,  et  al.  proposed  separate 
methods  for  polymeric  fluids  and  for  nonpolymeric  fluids.  However,  their  method 
does  not  include  nonhydrocarbons.  The  method  of  Kouzel  includes  all  fluids,  but 
as  shown  in  Table  40,  its  predictive  effectiveness  is  the  poorest  of  the  five 
cot  relations  evaluated.  This  is  particularly  true  with  polymer  blends,  resin 
blends  and  nonhydrocarbons. 

Equations  [45]  and  [42]  from  the  current  investigation  are  the  only 
correlations  that  cover  all  four  groups  of  fluids.  For  mineral  oils,  Equation 
[45]  is  as  good  as  Roeland's  method  and  Equation  [42]  has  an  edge  over  both. 

The  other  methods  are  less  effective.  It  is  important  to  point  out  that  the 
method  of  Roelands  requires  four  physical  property  values  and  four  additional 
correlations.  Equation  [45]  requires  only  three  physical  property  values  and 
one  additional  correlation  as  shown  on  Table  27.  Only  two  physical  property 
values  and  one  additional  correlation  are  required  for  Equation  [42].  For  pure 
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hydrocarbons.  Equation  [45]  has  the  lowest  average  absolute  deviation. 

Equation  [42]  is  a close  second  in  predictive  accuracy.  The  average  absolute 
deviations  of  the  other  methods  are  at  least  1.9  times  higher  than  the  average 
absolute  deviation  obtained  with  Equation  [45]. 


For  polymer  and  resin  blends,  the  method  reported  in  AFML-TR-304,  Part  II, 
has  a slightly  lower  average  absolute  deviation  than  does  Equation  [45],  That 
method  was  used  to  predict  only  17  out  of  34  available  polymer  and  resin  blends. 
Moreover,  that  method  requires  viscosity  information  for  both  base  oil  and 
finished  blends.  Roeland's  method  ranked  below  those  of  AFML-TR-70-304 , Part  II, 
Equation  [45],  and  Equation  [42]  in  predictability.  The  method  of  Roelands  was 
limited  by  the  assumption  that  the  viscosity-  pressure  index  of  the  blend  consists 
of  a fixed  60  percent  of  solvent  viscosity-pressure  index  and  40  percent  of  a 
hypothetical  polymer  viscosity-pressure  index  as  shown  on  Table  25.  The  method 
of  Roelands  also  requires  viscosity  information  for  both  the  solvent  and  the 
blend.  Kouzel's  correlation  is  veiy  poor  for  polymer  and  resin  blends.  For 
nonhydrocarbons,  Equation  T45]  again  has  the  lowest  average  absolute  deviation. 
Equation  [42]  and  the  method  of  Fresco  can  be  considered  about  equal  but  less 
effective  than  Equation  [451.  However,  the  number  of  nonhydrocarbons  studied  is 
too  small  to  be  significant  relative  to  the  differences  in  the  effectiveness 
between  Equation  [45]  and  [42],  and  the  method  of  Fresco. 


For  correlations  which  can  be  used  for  a wide  variety  of  fluids.  Equation 
[45]  has  the  lowest  average  absolute  deviation.  Equation  [42]  represents  a 
close  second  choice.  Equation  [45]  is  the  best  general  relationship.  If 
atmospheric  density  data  are  not  available.  Equation  [42]  is  the  one  to  use. 
For  the  mineral  oils  studied.  Table  40  shows  that  Equation  [42]  is  slightly 
better  than  Equation  [45].  However,  the  weakness  of  Equation  [42]  shows  up 
for  polymer  and  resin  blends. 


Although  in  general,  Equations  [45]  and  [42]  are  better  than  any  other 
methods  in  predicting  the  pressure  coefficients  of  fluids,  they  also  have 
limitations.  First  of  all,  the  equations  were  developed  for  fluids  which  have 
an  average  formula  weight  per  atom  in  the  neighborhood  of  4.5.  The  equations 
tend  to  predict  too  low  a pressure  coefficient  for  fluids  with  much  higher 
average  formula  weights  per  atom.  This  is  because  the  higher  density  values  of 
these  fluids  are  due  partly  to  the  higher  average  formula  weight  per  atom  and 
not  due  entirely  to  tighter  molecular  pa .king.  Both  highly  aromatic  fluids  and 
fluids  with  heavy  elements  are  in  this  category. 


Secondly,  the  ASTM  slope  was  assumed  to  be  independent  of  temperature. 
However,  Klaus,  Hersh,  Pohorilla  and  Fenske  (42,  43)  showed  that  the  ASTM  slope 
decreases  at  low  temperature  for  dense  center  molecules  in  the  synthetic  fluids. 
They  also  pointed  out  that  for  waxy  mineral  oils,  polymeric  fluids,  and  polymer- 
thickened  blends  the  ASTM  slope  increases  at  low  temperature.  Silicones  and 
organic  esters  generally  show  a decrease  in  ASTM  slope  at  high  temperatures. 

The  ASTM  slope  of  a chlorinated  aromatic  hydrocarbon  and  tricresyl  phosphate 
were  shown  to  increase  at  high  temperature.  Therefore,  for  some  of  the  fluids, 
the  ASTM  slope  depends  on  the  temperature  range  used  to  calculate  the  ASTM 
slope.  This  will  have  an  effect  on  the  predicted  pressure  coefficient  value. 


Thirdly,  Equations  [45]  and  [42]  do  not  fully  account  for  the  differences 
in  the  smoothness  of  the  molecules^. jJnder  constant  temperature,  pressure  and 


1 


average  formula  weight  per  atom  one  can  increase  the  kinematic  viscosity  either 
by  increasing  chain  length  and  branching  or  by  introducing  fused  rings.  If  the 
molecular  chain  is  stiff  and  its  branches  are  rigid,  a highly  branched  molecule 
may  have  the  same  ASTM  slope  and  viscosity  level  as  a molecule  made  up  of  smooth 
fused  rings.  The  result  is  that  molecular  interlocking  introduced  by  external 
pressure  will  be  more  pronounced  for  highly  branched  molecules.  As  a result  the 
highly  branched  molecule  will  have  a larger  pressure  coefficient  than  the  fused 
ring  molecule.  These  differences  are  not  fully  accounted  for  by  the  current 
models. 

Since  the  above  three  limitations  have  different  effects  on  the  fluid,  the 
net  result  may  not  be  clear.  Moreover,  in  examining  the  deviations  of  the 
individual  fluids,  one  has  to  keep  in  mind  that  the  discrepancy  of  the  experi- 
mental values  for  pressure  coefficient  among  different  investigators  can  be  as 
high  as  0.25  x 10-4  psig  (see  Table  15).  The  deviations  in  the  comparisons  show 
that  for  Equation  [45],  in  general  the  fluids  with  large  negative  deviations  are 
of  the  dense  center  type  molecules.  The  ASTM  slope  of  the  dense  center  molecule 
decreases  with  temperature.  Therefore,  the  ASTM  slope  based  on  100°  and  210°F 
is  larger  than  the  ASTM  slope  based  on  a wider  temperature  range.  When  the 
ASTM  slope  is  higher  than  anticipated  for  a molecular  structure.  Equation  [45] 
predicts  a pressure  coefficient  larger  than  the  experimental  one.  This  causes  a 
negative  deviation. 

On  the  other  hand,  the  fluids  which  have  large  positive  deviations  in  the 
comparison  are  waxy  mineral  oils,  polymeric  fluids  and  polymer-thickened  blends. 
The  ASTM  slopes  of  these  type  of  fluids  tend  to  increase  at  low  temperature. 
Therefore,  the  ASTM  slope  based  on  100°  - 2 10CF  is  smaller  than  the  ASTM  slope 
based  on  a wider  temperature  range.  Because  the  ASTM  slope  is  smaller  than  it 
should  be,  predicted  pressure  coefficients  are  too  small.  This  causes  a positive 
deviation.  The  higher  average  formula  weight  per  atom  of  some  of  the  compounds 
also  causes  a positive  deviation. 

By  examining  the  comparison  of  the  various  methods  in  predicting  the 
pressure  coefficients  of  nonhydrocarbons  which  lie  beyond  the  limits  designed 
for  Equations  [45]  and  [42]  in  Table  41,  a better  understanding  of  the  limita- 
tions of  the  equations  is  found. 

All  the  fluids  in  Table  41  have  average  formula  weights  per  atom  above  4.5. 
It  is  interesting  to  note  that  Equation  [45]  predicts  some  of  the  fluids  quite 
accurately.  These  fluids  are  the  ones  with  low  enough  atmospheric  viscosities 
sc  that  the  difference  in  molecular  packing  doesn't  have  much  effect  on  the 
pressure  coefficient  (see  Figure  21).  Once  again  the  large  positive  deviations 
are  due  to  higher  average  formula  weight  per  atom  and  the  low  100-210°F  ASTM 
slbpe  for  polymeric  fluid.  The  large  negative  deviation  is  due  to  higher  100'- 
210°F,  ASTM  slope  for  dense  center  material.  The  summary  of  the  data  on  Table  41 
is  shown  in  Table  43. 

The  data  on  Table  43  show  that  Equation  [42]  has  the  lowest  average 
deviation.  For  these  fluids,  Equation  [45]  is  not  effective  for  predictions. 

This  is  understandable  because  Equation  [45]  has  a strong  dependence  on  density. 
The  high  average  formula  weight  per  atom  distorted  the  meaning  of  density  as  a 
comparative  measurement  of  the  differences  m free  volume.  Therefore,  for 
fluids  with  average  formula  weights  per  atom  above  4.5,  Equation  [42]  is  the 
best  relationship  to  use  for  an  approximation  of  the  pressure  coefficients. 
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Table  1 

CALIBRATION  OF  GAGES  FOR  THE  PRL  PRESSURE  VISCOMETER 


5,000  psig  Capacity  Gage 
Dead  Weight  Average  Gage 

Gage  Pressure*,  Reading, 

psig  psig 

10,000  psig 
Dead  Weight 
Gage  Pressure* 
psig 

Capacity  Gage 

Average  Gage 
Reading, 
psig 

530 

545 

1,030 

1,040 

1,030 

1,045 

2,030 

2,045 

1,530 

1,545 

3,030 

3,075 

2,030 

2,045 

4,030 

4,078 

2,530 

2,558 

5,030 

5,085 

3,030 

3,045 

6,030 

6,100 

3,530 

3,550 

7,030 

7,090 

4,030 

4,045 

8,030 

8,095 

4,530 

4,543 

9,030 

9,130 

4,930 

4,950 

10,030 

10,150 

■k 

Dead  Weight  Gage 

has  an  accuracy 

of  + 1 psig 

i 


Figure  1.  SCHEMATIC  DIAGRAM  OF  PRL  PRESSURE  VISCOMETER  SYSTEM 


Figure  2.  CAPILLARY  UNIT 


Storage  Bulb 
Efflux  Bulb 
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Pressure  = One  Atmosphere 

Calibration  Fluid  = MLO  7372  Paraffinic 
Neutral  (9.902  Centistokes  at  100  F,  6.114 
Centistokes  at  130  F) 

• 100°F 


Head,  cm 


Figure  4.  CALIBRATION  CURVE  FOR  CAP II  LARY  P-81 


Table  2 


DETERMINATION  OF  THE  REPRODUCIBILITY  OF  THE  PRESSURE 
COEFFICIENT  MEASUREMENTS  IN  THE  PRL  PRESSURE  VISCOMETER 


Test  Fluid:  MLO  7710  = Di-2-ethylhexyl  Sebacate 

MLO  7372  = Paraffinic  Neutral 


Test  Temperature  = 100°F. 


Test  Fluid 

Gas -Liquid 
Contact  Time,  hr 
at  3,000  psig 

Pressure 
Coefficient , 
(1/Psig)x(10  ; 

Maximum 
Deviation, 
(l/psig)x(10  ) 

MLO  7710 

1 

1,048 

+ 

1 

1.051 

0.008 

1 

1.050 

1 

1.056 

MLO  7710 

4.5 

1.026 

+ 

4.5 

1.021 

0.005 

4.5 

1.024 

+ 

MLO  7710 

7.5 

1.002 

+ 

0.003 

7.5 

1.005 

MLO  7372 

1 

1.326 

+ 

0.014 

1 

1.312 

+ 

Figure  5.  THE  EFFECT  OF  DISSOLVED  GAS  ON  THE  PRESSURE  COEFFICIENTS 


RELATION  BETWEEN  SLOPE  m OF  MacCouil 
WALTHER  EQUATION  AND  ASTM  SLOPE 


Relation 
Between  m and 


ASTM  Slope  (C) 


EFFECT  OF  PRESSURE  ON  THE  VISCOSITY 
OF  PARAFFINIC  NEUTRAL  CMLO  7375) 


! 

i 

/ 


! 


i 


O MLQ  7375  (Mineral  Oil) 

A PRL  4561  (Mineral  Oil  Blend) 

O PRL  4^60  (Mineral  Oil  Blend) 

V PRL  4559  (Resin  Blend) 

Temperature:  100 °F 
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Figure  13.  EFFECT  OF  PRESSURE  ON  THE  VISCOSITY 
OF  A MINERAL  OIL,  MINERAL  OIL  BLENDS 
AND  A RESIN  BLEND  AT  SIMILAR  VIS- 
COSITY LEVEL 
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Based  on  pressure  viscosity  measurements  up  to  10,000  psig. 

Fluids  used  in  studies  included  in  Annual  Report  AFML-TR- 70- 304 , Part  II. 


Fluid 

Designation 

Solvent 

Description 

Centistoke 
100  °F 

Viscosity 

210°F 

ASTM 

Slope 

7685 

Naphthenic  Gas  Oil 

3.27 

1.24 

0.865 

2859 

Naphthenic  Mineral  Oil 

8.80 

2.28 

0.849 

7516 

Naphthenic  White  Oil 

78.86 

8.20 

0.769 

7788 

Paraffinic  Mineral  Oil 

164.7 

15.26 

0.675 

7710 

Di-2-ethylhexyl  Sebacate 

12.57 

3.32 

0.700 

Table  7 


SOLVENTS  USED  IN  EARLIER* 
VISCOSITY-PRESSURE  STUDIES 


Table  8 


' 

I 

I 

1 


t i 


[ \ 


M i 


1 

* 


BRIEF  DESCRIPTION  OF  THICKENERS  USED  IN 
EARLIER*  VISCOSITY-PRESSURE  STUDIES 


Thickener  Designation 

Description 

PRL  2905 

High  iriolecular  weight  polymethacrylate  (50.4 
wt.  per  cent  Acryloid-160  in  a naphthenic  gas 
oil)  . 

PRL  2906 

Low  molecular  weight  polvmethacrylate  (47.7 
wt.  per  cent  Acryloid-/5  in  a naphthenic  gas 
oil)  . 

Ac- 30 7 

Low  molecular  weight  polymethacrylate  (52.3 
wt.  per  cent  Acryloid-25  in  di-2-ethylhexyl 
sebacate) . 

Paratone  cx-2 

Polysiobutylene . (33.6  wt . per  cent  in  a 

naphthenic  white  oil) . 

MLO  7760 

Super-refined  heavy  resin  with  molecular 
weight  of  approximately  2000  to  3000  and 
basically  paraffinic  hydrocarbon. 

*AFML-TR- 70-304,  Part 

II. 
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(Continued  on  next  page. 
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PSU  615  2-n-Butyl-3-n-hexyldecalin  15.28  0.885  0.8639 


First 


Second 


Third 


Fourth 


* **  *** 

ASTM  , V , D 

ASTM2,  V2,  D2,  (ASTM) (V) , (ASTM) (D) , (D) (V) 

ASTM3,  V3,  D3,  (ASTM)(V2),  (D) (V2) , (ASTM) (D2 ) , 
(ASTM2)(V),  (D2)(V),  (ASTM2)(D),  (ASTM) (V) (D) 

ASTM4,  V4,  D4,  (ASTM' 'V3) , (D) (V3) , (ASTM)  (D3)  , 
(ASTM2)(V2),(D2)(V2) , \STM2)(D2),  (ASTM3) (V) , (DJ) (V)  , 

(ASTM3) (D) , (ASTM2),  )(D),  (ASTM) (V2) (D) , (ASTM) (V) (D2) 


I«JJJ 


Table  11 

DEFINITION  OF  ABSOLUTE  DEVIATION  TERMS  USED  IN 
THE  ANALYSIS  OF  THE  DIFFERENT  CORRELATIONS 


Dev  = Experimental  Value  - Calculated  Value 


Avg  = — £ | Dev | 


Bias  = — £ Dev 
n 


Fy  T 
RMS  = /-I  Dev 
v n 


n = Number  of  points 


Sum  of  Residual  Squares  = £ Dev^ 


Standard  Error  of  Estimate  = / - \ v £ Dev^ 

/ (n-k) 


k = Number  of  parameters 


60 


- rr~~m 


Legend 


• Addition  of  Variable 


A Deletion  of  Variable 


Figure  14.  STANDARD  ERROR  OF  ESTIMATE  VERSUS  NUMBER  OF  VARIABLES 
FROM  THE  MULTIPLE  LINEAR  REGRESSION  ANALYSIS 


log  (Atmospheric  Viscosity) , log  (centistokes) 
Atmospheric  Density,  gm/c.c. 
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Numbers  on  curves  indicate  ASTM  Slope  0.95  0.85  0.75 
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Figure  16.  LIQUID  PRESSURE  COEFFICIENT  PREDICTION  CHART  BASED  ON  EQUATION  [45]  AT  ATMOSPHERIC  DENSITY  OF 
0.85  gm/cc  (Concluded  on  next  page.) 
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Figure  16.  LIQUID  PRESSURE  COEFFICIENT  PREDICTION  CHART  BASED  ON  EQUATION  [45]  AT  ATMOSPHERIC  DENSITY  LEVEL 


Figure  17.  LIQUID  PRESSURE  COEFFICIENT  PREDICTION  CHART  BASED  ON  EQUATION  [45]  AT  ATMOSPHERIC  DENSITY  LEVEL 
OF  0.95  gm/cc  (Concluded  on  next  page.) 


Figure  17.  LIQUID  PRESSURE  COEFFICIENT  PREDICTION  CHART  BASED  ON  EQUATION  [45]  AT  ATMOSPHERIC  DENSITY  LEVEL 
OF  0.95  gm/sec  (Concluded) 
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Atmospheric  Density  = 0.85  gm/cc 

Numbers  on  curves  indicate  atmospheric 
viscosity 

2000 I 


ASTM  Slope 

EFFECT  OF  ASTM  SLOPE  AT  CONSTANT  ATMOSPHERIC 
VISCOSITY  ON  PRESSURE  COEFFICIENT  BASED 
EQUATION  [45] 


Figure  20.  OPERATIONAL  REGIONS  OF  EQUATION  [45] 


Figure  22.  LIQUID  PRESSURE  COEFFICIENT  PREDICTION  CHART  BASED  ON  EQUATION  (42)  (Concluded  on  next  page. 
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Figure  22.  LIQUID-PRESSURE  COEFFICIENT  PREDICTION  CHART  BASED  ON  EQUATION  [42 ] (Concluded) 
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Equation  [45]  was  used  to  predict  the  pressure  coefficients  at  100  F,  and  Equation  149]  was  use 
correct  for  the  temperature  effects. 
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* Deviation  = Experimental  Value  - Predicted  Value 

**  Equation [45] was  used  to  predict  the  pressure  coefficients  at  100°F  and  Equation  49  was  used  to  correct 
for  the  temperature  effects 
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Deviation  = Experimental  Value  - Predicted  Value 

Equation  [45]  was  used  to  predict  the  pressure  coefficients  at  100°F,  and  equation  [49]  was  used  to  correct 
for  the  temperature  effects. 


-Defined  in  Table  11 

** Equation  [45]  was  used  to  predict  the  pressure  coefficients  at  100°F,  and  equation  [49]  was  used  to 
correct  for  the  temperature  effects 
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COMPARISON  OF  VARIOUS  METHODS  IN  PREDICTING  PRL  PRESSURE 
COEFFICIENT  DATA  TROM  CURRENT  INVESTIGATION  AND  FRESCO  (27) 


m 

vO 

vO 

co 

vo 

oo 

00 

00 

CM 

m 

CO 

CO 

•H 

CO 

o 

CM 

oo 

m 

CO 

00 

ON 

CM 

CM 

ON 

vo 

MT 

vO 

00 

r*^ 

CM 

CM 

vO 

00 

rH 

rH 

rH 

O 

CVI 

-3- 

rH 

rH 

CM 

VO 

CM 

o 

CM 

vO 

rH 

rH 

O 

1 

o 

1 

o 

1 

O 

1 

rH 

1 

o 

i 

o 

i 

o 

1 

o 

1 

<?c?c? 

rH 

1 

o 

1 

o 

1 

O 

1 

o 

O 

O 

i 

O 

<r 

r^. 

vO 

vO 

Ht 

>3- 

o 

vO 

oo 

o 

m 

<r 

VO 

o 

vO 

vO 

CO 

rH 

CM 

CM 

vO 

vO 

m 

ON 

CM 

o 

00 

CM 

o 

o 

O 

fH 

rH 

O 

O 

O 

O 

O 

o 

o 

CO 

rH 

CM 

1 

1 

CM 

1 

1 

o 

1 

o 

1 

O 

o 

1 

<? 

O 

O 

1 

i 

? 

o 

1 

o 

i 

o 

1 

o 

1 

? 

o 

1 

«H 

n* 

<r 

rH 

rH 

CM 

oo 

<r 

oo 

vO 

ON 

rH 

o 

O 

in 

CM 

vO 

rH 

co 

vO 

rH 

CO 

CO 

o 

rH 

rH 

CO 

O 

CM 

<r 

ON 

CM 

m 

CM 

rH 

O 

O 

o 

O 

O 

o 

rH 

o 

o 

rH 

o 

o 

O 

O 

o 

m 

VO 

o 

o 

oo^oomndsr 

NinfOHiONfOH 

OOOrHrHOrHO 


\OvonmHNNincMOOH(n 

HHrj^NNr^cM'd-oooooo 

oooHNOOd<r<rHd 


ONCMincocovocM<r 

ooaocMcor^cMcocM 

OOOrHrHO*HO 


oooooooo 


vor^vooNiHvOvomcMOvon- 

OOOCMCOOOCMCMCOOO 


nnvoocoaNMvo 

jnnoor^oosfooso 

fHcovor^ONvo^tm 


fncoHvor^mrNNNrNvoH 

mm*<j-r^oomooir><j-^o<T 

HHHHHHHH  CM  CS  H H 


OOOOOO'jnON 

ojo<rH<ffOom 

oocONr^vONrsvo 



oooooooo 


<fCOOvOrnOCNON<-OOrHv£) 

NNfS\OvOCONvOfsrsNl/) 

oooooooooooo 


a\vOMOvOHnN 
in  H vO  H 10  Nin  in 

H^mNMnmin 

oooooooooooooooo 


omoovonninvO'jinmin 
mrHvor^^-<rvoor^oOfHr^ 
cMcocoinr-^vor^cMcocMOm 
oooooooooooooooooooo  on  on 

oooooooooooo 


<f  o fn  cm  cm  *h 
OOOvCOCOCMCOOOvO 

'JOSONONOOONlO'J 

co  on  co  co  <r  m 
<r 


00  o CO  H a\  O'  vO 

NOH^vOOOvOO\'fNinvD 

N^OvOOCOCOO'<finn|CO 
H H N ^ CO  N N M H H 00 


CM  m O v£>  rH  o on 
N N (M  Sf  ^ vO  m 
co  co  »h  cm  in  m m 
vor^r^r^r^^r*^^ 


ONr^oo»H^aovo^-co«<rcM<r 

Hinmvoo'inHmcMcMHoo 

minm'OOoomNHHHvo 

NNNfNNCMrsr^NrsNN 


*Deviation  = Experimental  value  - Predicted  value  (b)  Equation  [ ^» 2 ] (d)  Roelands,  et . al.  (23) 

(a)  Equation  [ 45 ] (c)  Fresco  (27)  (e)  Kouzel  (26) 
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Deviation  = Experimental  Value  - Projected  Value  (b)  Equation  [42]  (d)  Roelands,  et  al.  (23) 

(a)  Equation  [45]  (c)  Fresco  (27)  (e)  Kouzel  (26) 
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(e)  Roelands,  et . al.  (23) 

(f)  Kouzel  (26) 


COMPARISON  OF  VARIOUS  METHODS  IN  PREDICTING  PRESSURE 
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Table  42 


BRIEF  DESCRIPTIONS  OF  NON-HYDROCARBONS  IN  TABLE  41 


Fluid 

Designation 


Description 


MLO  7742 
MLO  7723 
MLO  7428 
MLO  7017 

MLO  7756 
MLO  7743 
MLO  7741 
MLO  7668 
MLO  7639 
MLO  7061 


Halocarbon  Oil  14-25 
DuPont  PR  143 

4- ring  polyphenyl  ether 

Improved  Lubricity  Silicone  Fluid  (G.E.  Silicone 
81406) 

Halocarbon  Oil  208 
Halocarbon  Oil  11-21 
Halocarbon  Oil  11-14 

5- ring  polyphenyl  ether 

Phosphate  Ester  Hydraulic  Fluid  (Skydrol  500  1) 
Chlorinated  Biphenyl 


C.  The  Effects  of  Chemical  Reactions  in  Boundary  Lubrication.  The  term 
"boundary  lubrication"  was  first  introduced  by  Sir  William  Hardy  (45)  in  1922  in 
describing  a region  where  "the  solid  faces  are  near  enough  to  influence  directly 
the  physical  properties  of  the  lubricant,  and  the  friction  depends  not  only  on 
the  lubricant  but  on  the  chemical  nature  of  the  solid  boundaries."  Recently, 
Dowson  (46)  defined  boundary  lubrication  as  the  "surface  interaction  between  a 
monomolecular  layer  of  boundary  lubricant  and  the  solids  dominate  the  operation 
of  the  contact.  Hydrodynamic  effects  are  negligible,  and  there  is  considerable 
asperity  contact,"  While  the  words  "monomolecular  layer"  are  debatable,  this 
definition  is  more  precise  than  Hardy's. 

Campbell  (47)  further  clarified  the  term  by  stating,  "Boundary  lubrication 
is  lubrication  by  a liquid  under  conditions  where  the  solid  surfaces  are  so 
close  together  that  appreciable  contact  takes  place  between  the  asperities.  The 
friction  and  wear  are  influenced  predominantly  by  interaction  between  the  lubri- 
cant and  the  solid,  and  the  bulk  flow  properties  of  the  lubricant  play  little  or 
no  part  in  friction  and  wear  behavior."  The  word  "interaction"  will  be  qualified 
as  "chemical  interaction"  in  this  work. 

The  object  of  this  study  is  to  investigate  the  effects  of  the  chemical 
interactions  between  the  solid  surfaces  and  the  lubricant.  The  lubricant  will 
be  restricted  to  mineral  oils. 

Lubrication  is  a very  complex  phenomenon.  It  is  the  reduction  of  friction 
and/or  wear  by  a solid,  liquid,  or  gas  film  existing  between  two  surfaces  moving 
against  one  another.  In  this  work,  only  lubrication  between  metal  surfaces  by 
liquid  hydrocarbons  will  be  considered. 

Lubrication  can  generally  be  classified  into  three  regions:  hydrodynamic, 

elastohydrodynamic,  and  boundary  region.  A classical  Stribeck  (48)  diagram,  as 
shown  in  Figure  23,  illustrates  these  regions.  At  the  minimum  of  the  curve 
corresponding  to  point  A,  the  coefficient  of  friction  is  about  0.002  and  yN/P 
quotient  (where  y is  the  oil  viscosity  in  centipoises,  N the  journal  speed  in 
rpm,  and  P the  normal  load  per  projected  bearing  area  in  psi.)  varies  from  1 to 
10  according  to  the  oil  type,  bearing  geometry  and  material,  and  surface  finish. 
For  values  of  uN/P  above  the  minimum,  the  curve  is  quite  linear.  In  this  region, 
the  oil  film  is  thick  enough  so  that  the  solid  surfaces  are  completely  separated 
from  one  another.  The  coefficient  of  friction  depends  on  the  fluid  viscosity, 
sliding  speed,  and  beating  geometry.  This  region  is  the  most  desirable  state  of 
lubrication  due  to  its  inherent  protection  of  the  solid  surfaces  by  the  fluid. 
There  is  no  wear  on  the  bearings.  The  lubrication  mechanism  in  this  region  is 
fairly  well  explained  by  the  hydrodynamic  theory  which  is  summarized  by  Hersey 
(49).  Its  development  can  be  traced  from  Sir  Issac  Newton  (1687)  to  Petroff 
(1883),  Tower  (1883)  and  Reynolds  (1886).  The  journal  bearing  in  an  automobile 
operates  hydrodynamically  at  high  speeds. 

According  to  Stribeck,  the  lubricating  film  begins  to  break  down  and 
asperity  contact  between  the  surfaces  occurs  at  point  A (Figure  23) . The  region 
between  points  A and  B in  Figure  23  has  many  names:  elastohydrodynamic,  quasi- 

hydrodynamic, partial  fluid,  and  thin  film  lubrication.  All  these  names  indicate 
one  thing,  that  is  the  solid  boundaries  begin  to  exert  an  Influence  on  the 
lubrication  at  point  A.  Classically,  it  has  been  said  that  the  load  in  region  AB 
is  supported  partly  by  hydrodynamic  action  of  the  lubricant  and  partly  by  the 


boundary  films.  Recently,  Dowson  (46)  proposed  that  in  this  region,  the  local 
elastic  distortion  of  the  solids  provided  coherent  hydrodynamic  films;  and 
asperity  interaction  was  prevented.  The  thin  oil  film  was  subjected  to  very 
high  pressure  and,  as  a result,  the  local  viscosity  increased  greatly.  Experi- 
mental data  from  Archard  and  Kirk  (50)  was  used  to  demonstrate  this  theory.  As 
N/P  decreases  below  A,  the  proportion  of  hydrodynamic  support  decreases  until  at 
and  below  point  B,  lubrication  is  pure  boundary.  In  region  BC,  the  coefficient 
of  friction  is  Independent  of  the  load;  this  is  called  Amonton's  Law.  It  should 
be  stressed  that  the  curve  shown  in  Figure  23  is  only  a representation  of  a set 
of  data.  For  any  other  system,  the  points  A and  B may  change  their  positions. 

The  regions  discussed  above  are  not  very  clear  cut.  In  the  context  of  this 
study,  the  boundary  lubrication  region  has  been  previously  defined  as  the 
conditions  under  which  the  distance  between  the  solid  boundaries  is  close  enough 
that  chemical  interactions  (as  well  as  physical)  between  the  solids  and  the 
lubricant  play  an  important  role  in  the  lubrication  mechanism  Thus  the  region 
ABC  in  Figure  23  would  be  included  in  this  study. 

1.  Solid  Surfaces  and  Dry  Friction.  Since  the  solid  surface  plays  an 
important  role  in  the  mechanisms  of  boundary  lubrication,  the  nature  of  the 
surface  will  be  examined  On  a microscopic  scale,  all  metal  surfaces  used  in 
lubrication  are  extremely  rough  no  matter  how  finely  they  have  been  polished. 

An  indication  of  the  roughness  generated  by  the  various  machining  operations  is 
shown  in  Table  44.  Roughness  characteristics  have  been  studied  by  Shaw  and 
Macks  (51).  With  the  most  careful  polishing,  flatness  can  be  carried  to  about 
5 y in.  (1270A).  The  atomic  spacing  in  a metal  is  about  0.01  y in.  (2.5A). 

There  have  been  many  attempts  both  theoretically  and  experimentally  to 
characterize  the  surface.  The  resulting  knowledge  is  valuable  but  a satisfactory 
model  has  yet  to  emerge.  As  Moore  (52)  pointed  out,  most  models  failed  m 
characterizing  the  surface  using  only  one  parameter  such  as  the  average  height 
of  the  asperities  or  the  root  mean  square  (RMS)  of  the  distance  from  a mean  plane. 
In  order  to  describe  the  surface  fully,  one  has  to  take  into  account  the  height 
of  the  asperities;  the  void  volume  of  the  surface  (the  density  of  the  asperities); 
and  the  sharpness  of  the  peaks  given  by  the  second  derivative  of  the  surface 
profile.  There  are  many  experimental  techniques  that  can  be  used  to  determine 
a surface  profile.  Excellent  reviews  in  this  area  are  given  by  Green  (53),  and 
Brown  et  al.  f54).  An  indication  of  the  surface  roughness  of  a metal  has  been 
determined  by  Furey  (55)  using  a stylus  instrument  on  52-100  steel. 


When  two  dry  surfaces  meet,  the  contact  area  therefore  is  not  equal  tc  the 
corresponding  surface  area  of  the  metal.  The  problem  of  determining  contact 
areas  for  metal  surfaces  with  idealized  asperity  shapes  and  distributions  has 
been  investigated  by  many  researchers.  Greenwood  (56/  showed  that  if  the  real 
area  of  contact  between  surfaces  was  determined  by  ideal  plastic  flow  of  the 
microcontacts,  the  area  of  contact  was  proportional  to  load.  For  elastic 
deformation,  the  average  size  of  a microcontact  was  almost  constant  and 
independent  of  the  load  if  the  statistical  distribution  of  asperity  heights  were 
considered. 

The  nature  of  the  surface  also  depends  on  its  past  history  (57).  A surface 
being  polished  is  different  from  a surface  being  ground  using  the  same  material. 
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One  can  study  the  surface  topography  with  LEED  (low  energy  electron  diffraction) 
method  using  electrons  whose  range  in  the  solid  may  be  only  a few  tenths  of  an 
Angstrom  unit.  The  resulting  diffraction  pattern  is  primarily  that  of  the  surface 
region.  Using  this  method,  one  finds  grinding  leads  to  primarily  mechanical 
attrition  of  the  surface  without  changing  its  molecular  crystallinity  while 
polishing  leaves  a fairly  deep  and  nearly  amorphous  surface  layer,  known  as  the 
Beilby  layer,  Beilby  (58)  showed  that  such  a layer  appeared  amorphous  under  the 
microscope  and  had  the  general  appearance  of  a film  of  viscous  liquid  that  not 
only  had  covered  the  surface  smoothly  but  had  flowed  into  the  surface  irregular- 
ities such  as  cracks  and  scratch  marks.  Raether  (59)  concluded  that  the 
Beilby  layer  was  actually  microcrystalline.  The  formation  of  this  layer  may  be 
due  to  the  high  temperature  generated  from  the  polishing  operation  and  the  sub- 
sequent metal  softening  if  not  actually  melting. 

Surface  defects  constitute  another  aspect  of  the  surface  condition.  Actual 
crystals  are  not  perfect;  they  may  consist  of  microcrystalline  domains  bounded 
by  slip  planes,  or  they  may  have  lattice  irregularities  surrounded  by  good 
crystals.  The  irregularities  may  involve  a screw  or  spiral  dislocation  (60,  61). 

It  is  not  uncommon  for  20  percent  of  the  volume  of  a crystal  to  be  occupied  by 
dislocations  (62). 

a.  Thermodynamics  of  Surfaces.  Surface  tension  is  defined  as  the  work 
spent  in  forming  a unit  area  of  surface  (alternatively  called  surface  free  energy). 
The  surface  stress  is  defined  as  the  work  spent  in  stretching  the  surface.  In  a 
liquid,  the  surface  mobility  of  the  atoms  is  great  and  the  surface  tension  is 
equal  to  the  surface  stress.  In  a solid,  these  two  quantities  are  not  necessarily 
equal  due  to  the  relative  immobility  of  the  atoms  at  the  surface. 

It  is  helpful  to  imagine  that  the  process  of  forming  a fresh  surface  of  a 
monatomic  substance  is  divided  into  two  steps:  first,  the  solid  or  liquid  is 

cleared  so  as  to  expose  a new  surface,  keeping  the  atoms  fixed  in  the  same 
positions  that  they  occupied  when  in  the  bulk  phase;  second,  the  atoms  in  the 
surface  region  are  allowed  to  rearrange  to  their  final  equilibrium  positions. 

In  the  case  of  liquid,  the  two  steps  occur  as  one,  but  with  the  solid,  the  second 
step  may  occur  only  slowly  through  surface  diffusion.  Thus  It  may  be  possible 
to  stretch  or  to  compress  the  solid  surfaces  without  changing  the  number  of  atoms 
in  it,  only  their  distances  apart.  Under  this  condition,  a surface  stress  would 
be  present.  Dunning  (60)  postulated  that  the  dilatant  stress  in  small  crystals 
could  be  relieved  by  the  appearance  of  dislocations.  A cube  of  ideal  crystal 
subjected  to  compressive  surface  stress  suffers  deformation  equivalent  to  applying 
traction  to  each  edge  of  the  cube.  This  stress  could  be  relieved  by  rows  of 
dislocations  in  which  lines  of  atoms  are  missing,  or  by  means  of  vacancies  in 
surface  planes  (63).  It  thus  appears  that  with  a non-equilibrium  surface  whose 
surface  stress  and  surface  free  energy  values  are  not  equal,  there  is  a source 
of  surface  Imperfections  in  addition  to  those  representing  the  surface  eruption 
of  bulk  imperfections  The  ability  of  a surface  to  relieve  a nonequilibrium 
stress  may  be  very  temperature  dependent.  At  temperatures  near  the  metal's 
melting  point,  both  bulk  and  surface  diffusion  generally  become  appreciable.  For 
copper  (58)  at  725‘C,  it  takes  0.1  sec.  for  a copper  atom  in  the  bulk  to  diffuse 
100  A,  At  room  temperature,  it  will  take  102  7 sec.  The  field  emission  microscope 
has  shown  that  surface  migration  can  be  observed  at  as  low  as  half  the  evaporation 
temperature  of  the  solid. 
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At  the  surface,  not  all  atoms  possess  the  same  energy  level.  Those  atoms 
or  molecules  at  the  rugged  asperity  tips  or  the  edge  of  a surface  crack  or 
imperfections  of  the  crystal  lattice  possess  higher  energy  and  surface  mobility 
than  the  average  surface  atom  which  has  a normal  number  of  nearest  neighboring 
atoms  These  high  energy  atoms  are  much  more  reactive  and  provide  the  most 
active  adsorption  sites  for  polar  molecules. 

In  a situation  where  constant  surface  contact  and  rubbing  take  place,  large 
shear  stresses  are  always  present.  The  resulting  high  temperatures  at  the 
asperity  tips  from  the  frictional  energy  will  increase  a surface  atom's  mobility 
and  accelerate  the  rate  to  relieve  the  surface  stress  by  creating  more  dislocations 
or  imperfections.  This  results  in  more  chemically  active  sites  for  adsorption  and 
reaction. 


b.  Dry  Friction  of  Metals.  When  two  solid  surfaces  rub  against  one 
another,  frictional  heating  is  generated.  Depending  on  the  load  and  sliding 
speed,  very  high  temperatures  can  occur  at  the  surface,  Because  of  the  surface 
roughness,  the  real  contact  area  is  smaller  than  the  surface  area  of  the  two 
surfaces.  Therefore  the  heating  is  further  localized  around  the  asperities.  The 
surface  temperature  had  been  measured  (64)  by  measuring  the  electrical  potential 
difference  between  the  rubbing  surfaces.  In  the  case  of  a steel  rider  gliding 
over  different  metals,  it  can  be  shown  that  a temperature  rise  of  several  hundred 
degrees  can  occur.  The  temperature  rise  levels  off  at  the  melting  temperature  of 
the  lower  melting  metal.  If  one  of  the  surfaces  is  transparent,  l.e,,  glass,  the 
local  temperature  can  be  determined  by  infrared  spectrophotometry  In  this  way, 
local  temperatures  as  high  as  1200cC  have  been  measured  for  a steel  on  glass 
system  (65). 

The  oxide  film  and  the  Beilby  layer  are  the  top  most  surface  layers  of  a 
typical  metal  surface  (66).  These  layers  actually  define  the  surface  properties 
of  the  solid  surface,  The  oxide  film  is  usually  of  low  mechanical  strength  and 
can  be  sheared  off  under  high  shear  stress.  The  Beilby  layer  may  be  harder  than 
the  bulk  material.  Moore  (67)  measured  the  thermo-electric  output  between  a wear 
pin  (ingot  iron)  and  a rough  a-silicon  carbide  surface.  The  temperature  at  the 
junction  was  found  to  range  from  300'  to  900  C.  The  peak  temperatures  found  are 
high  enough  to  cause  ferrite-austenite  transformation  In  a peatlitic  steel  and 
the  subsequent  cooling  rate  fast  enough  to  cause  austenite  to  transform  to 
martensite,  which  Is  harder  than  steel.  Due  to  che  fact  that  the  hot  spot 
temperature  exists  only  at  the  site  of  real  contact  and  the  different  temperatures 
reached  at  different  asperity  contacts,  the  martensite  formation  exists  only  at  a 
few  spots.  This  may  explain  plowing  and  scratching  sometimes  observed  in 
boundary  lubrication  in  wear  tests  under  severe  conditions. 


There  have  been  many  attempts  to  predict  theoretically  the  hot  spot 
temperatures  in  dry  sliding.  Blok  (68)  and  Jaeger  (69)  analyzed  the  heat  transfer 
problem  and  derived  several  expressions  to  estimate  the  hot-spot  temperature  at 
the  junction.  They  assumed  adiabetic  heating  in  different  geometries  with 
different  simplifying  assumptions  . Arihar-'  (70)  summarized  the  equations  and 
popularized  their  usage.  However,  due  to  the  extreme. y short  duration  of  the 
contact  and  the  small  contact  region  involved,  accurate  experimental  measurement 
has  been  very  difficult  to  obtain  The  accuracy  of  the  theoretical  prediction 
is  still  open  to  question.  Sevetal  assumptions  used  in  the  theoretical  analysis 
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are,  at  best,  quite  crude.  The  temperatures  are  calculated  on  the  assumption 
that  the  heat  is  generated  at  the  area  of  true  contact  (i.e,,  a single  area  of 
contact  is  regarded  as  a plane  source  of  heat).  The  estimation  of  the  true  area 
of  contact  again  depends  on  whether  plastic  deformation  or  elastic  deformation 
is  assumed,  On  the  other  hand,  the  theoretical  result  indicates  a maximum 
possible  temperature  rise,  at  least,  qualitatively.  The  predicated  temperatures 
are  always  much  higher  than  the  experimentally  measured  temperatures  in  dry 
sliding.  The  temperature  rise  is  predicted  to  be  from  700°C  to  1200°C  (70)  for 
the  dry  wear  of  steel  (0,52  percent  carbon  steel,  V.P.N.  250). 

Archard  (70)  also  pointed  out  the  effect  of  a surface  film  of  low  thermal 
conductivity  on  the  flash  temperatures.  A protective  film  such  as  an  oxide  or 
a boundary  lubricant  would  increase  the  surface  hot-spot  temperature  and  if  the 
surface  film  thickness  was  large  compared  with  molecular  dimensions,  the 
temperature  gradient  existing  in  the  film  could  be  easily  the  largest  transien*- 
temperature  in  the  contact  zone  and  might  be  several  times  greater  than  the 
flash  temperature. 

In  the  unlubricated  wear  of  steel,  Quinn  (71,  72,  73,  74)  proposed  an 
oxidative  wesr  model.  He  assumed  all  of  the  mild  wear  of  unlubricated  steel-on- 
steel  was  due  to  oxidation  of  the  steel,  i.e,,  the  steel  surface  oxidized  to  a 
critical  thickness  at  which,  under  high  shear,  produced  wear  particles.  The 
temperatures  at  which  the  steel  oxidized  were  those  hot  spot  temperatures  at  the 
■junction.  Quinn's  wear  equation  is  shown  below: 

W = [d  A Ap  e"(Qp/RT0)l/[^2p2Vf2]  [53] 

where: 


VJ  = wear  rate,  volume  removed  per  unit  sliding  distance, 
d * distance  along  which  a wearing  contact  is  made, 

A = real  area  of  contact, 

Ap  = Arrhenius  constant  for  parabolic  oxidation  of  steel, 

Qp  * activation  energy  for  parabolic  oxidation  of  steel, 

R = gas  constant, 

T0  - temperature  of  oxidation, 

£ * critical  oxide  film  thickness  at  each  asperity  contact, 

p = density  of  the  oxide  film, 

V = speed  of  sliding,  and 

f = mass  fraction  of  oxide  film  which  is  oxygen. 

Detailed  derivation  of  this  equation  is  shown  in  Ref.  (71). 

Equation  [53]  can  be  criticized  that  it  contains  too  many  ill-defined 
parameters  such  as  d,  £,  A,  and  T0,  Quinn  (71)  recognized  this  but  he  agreed 
that  independent  estimates  could  be  made  on  d,  £,  and  A;  therefore,  T0  could  be 
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calculated  from  equation  [53]  with  some  wear  rate  data,  Ap  and  Qp  values  were 
obtained  from  the  static  oxidation  data  on  high  purity  steel  by  Kubaschewski  (75] 
Quinn  compared  the  oxidation  temperature  calculated  from  equation  [53],  830°C 
(1526  F)  with  Archard's  (70)  flash  temperature  theory,  1250C,  and  also  with 
Furey's  (77)  experimental  measurement  of  the  contact  temperature,  350^0  (662'F) 
from  dynamic  thermocouples  which  Quinn  (73)  himself  also  obtained  experimentally 
using  x-ray  diffraction  patterns  analyzing  wear  product  composition  (Fex0  ) as 
function  of  temperatures. 

The  discrepancy  among  the  temperatures  is  not  surprising.  Static  oxidation 
data  may  not  be  compatible  witn  dynamic  wear-oxidation.  In  the  wearing  process, 
new  surfaces  and  the  thickness  of  the  oxide  film  are  constantly  changing.  With 
three  parameters  d,  £,,  A cxudeiy  estimated,  one  is  surprised  to  see  that  the 
calculated  temperature  falls  within  the  reasonable  limits  (3505C  - 1250°C) . 

The  steel  Quinn  used  was  a low  alloy  medium  carbon  steel  similar  to  the 
steel  used  in  large  marine  turbine  shaft  bearings.  The  steel  composition  is  as 
follows: 


Carbon 

0.42% 

Nickel 

2.54% 

Silicon 

0.28% 

Chromium 

0 71% 

Manganese 

0.  62% 

Molybdenum 

0 . 55% 

Sulfur 

016% 

Copper 

0.16% 

Phosphorus 

0.015% 

Earles  (78,  79,  80,  81)  and  his  coworkers  also  studied  the  wear  character- 
istics of  steels  in  relation  to  surface  temperatures  and  oxidation.  While  the 
linear  sliding  speeds  in  Quinn's  work  ranged  from  1 m/sec.  to  8 m/sec.,  Earles 
(81)  used  much  higher  sliding  speeds,  20-100  m/sec.  Earles  showed  that  at  the 
same  load,  when  the  sliding  speeds  varied,  the  coefficient  of  friction  and  wear 
went  through  a transition  from  one  type  of  behavior  to  another,  which  were 
related  to  temperature  dependent  changes  in  the  metal  and  oxide  properties.  He 
showed  that  at  higher  sliding  speeds,  hence  higher  contact  temperatures,  FeO 
instead  of  Fe203,  FejO^  was  formed  at  the  surface , Fe203  and  Fe304  are  known 

to  have  good  friction  and  wear  pr_perties,  while  FeO  does  not.  The  properties 
of  the  oxides  are  presented  in  Table  45. 

Tenwick  and  Earles  (80)  derived  an  equation  correlating  the  average  suriace 
temperature  (not  the  contact  or  flash  temperature  as  in  Quinn's  work),  load, 
sliding  speed  with  oxidative  wear  of  steels  Also  instead  of  using  a parabolic 
oxidation  model  to  des.ribe  the  oxidation  mechanism,  they  used  a linear  oxidatior 
model.  However,  as  Berry  (82,  83)  pointed  out,  the  two  models  were  essentially 
identical  except  ter  the  interpretation  of  the  temperatures  at  which  oxidation 
occurs » 

Recently,  Amsallem  et  ai.  (84)  studied  the  wear  of  sintered  iron  in  dry 
sliding  both  in  air  and  in  a nitrogen  atmosphere.  The  result  showed  that 
different  oxides  having  dulerent  wear  characteristics  were  formed  as  the  load 
varied  under  constant  sliding  speed  The  wear  debris  were  analyzed  with  x-ray 
and  electron  diffraction  techniques.  From  the  composition  of  the  wear  debris, 
the  temperature  was  estimated  to  be  570  C at  80  N loading 


2,  Chemistry  of  Lubrication  The  wear  resulting  from  lubricated  contacts 
differs  from  that  of  dry  frictional  wear.  The  nature  and  the  composition  of  the 
lubricant  play  a major  role  in  the  wear  characteristics  of  the  contact  junction 


The  wear  also  depends  on  other  parameters  such  as  the  relative  speed  of  the  moving 
parts,  load,  metal  compositions,  oxygen  concentration  in  the  oil  and  the  environ- 
ment, water  concentration,  operating  temperature  and  geometric  configuration  of 
the  contact  surfaces.  In  this  section,  the  plausible  or  probable  chemical 
reactions  that  could  take  place  in  boundary  lubrication  are  discussed. 


There  are  hundreds  of  chemical  species  in  a typical  lubricating  oil.  Many 
of  them  have  been  identified  only  in  the  past  twenty  years.  In  an  investigation 
which  started  in  1927,  the  American  Petroleum  Institute  Research  Project  6 has 
isolated  264  chemical  compounds  from  one  petroleum  oil.  When  the  lubrl.ating  oil 
is  subjected  to  high  local  temperatures  and  extreme  pressures  in  boundary  lubri- 
cation, the  possible  chemical  reactions  that  could  result  are  literally  numerous. 
In  addition  to  the  severe  environment  that  the  oil  encounters,  the  clean  metal 
surface  generated  in  the  wear  process  and  the  chemical  additives  usually  added  in 
commercial  oils  further  complicate  the  situation.  Therefore  it  would  be  futile, 
if  not  impossible,  to  study  the  elementary  mechanisms  of  the  chemical  reactions 
that  are  taking  place  in  boundary  lubrication.  However,  some  groups  of  reactions 
that  play  a key  role  in  wear  could  be  examined  and  the  overall  chemical  effects 
in  boundary  lubrication,  hopefully  can  be  assessed  accordingly. 


a.  The  Chemical  Composition  of  Lubricants.  Lubricating  oils  from 
petroleum  are  composed  mainly  of  hydrocarbons  together  with  small  amounts  of 
organic  compounds  of  sulfur,  nitrogen,  and  oxygen  and  traces  of  organometallic 
compounds  of  vanadium,  nickel,  iron  and  copper.  The  hydrocarbons  can  be  classi- 
fied into  three  groups:  paraffins  (saturated  straight  or  branched  chains), 

naphthenes  (saturated  cycloparaffins),  and  aromatics  (with  benzene  rings).  Some 


typical  structures  in  lube  oils  are  illustrated  in  Figure  24.  The  petroleum  oil 


fractions  used  in  lubrication  comprise  principally  compounds  containing  from  20 


to  40  carbon  atoms  per  molecule,  especially  from  C20  to  C25. 


The  number  of  possible  chemical  compounds,  according  to  stereochemistry  for 
C20  is  almost  infinite  (20' 0 possible  variations).  Chemical  composition  analysis 
on  a number  of  lubricants  in  a recent  study  (86)  indicate  the  approximate 
composition  according  to  the  classes  of  compounds.  There  are  many  molecular 
species  in  each  of  the  classes  As  is  indicated,  molecules  sometimes  have 
paraffinic  structures  together  with  aromatic  rings  or  one  single  molecule  can 
have  all  three  classes  of  compounds.  Therefore,  from  a chemical  point  of  view, 
it  is  important  to  realize  that  in  a "paraffinic"  oil,  all  three  classes  of 
compounds  exist,  although  the  physical  properties  are  influenced  most  strongly  by 
the  paraffinic  portions. 

During  the  past  twenty  years,  high  temperature  mass  spectrometry  has  become 
a powerful  tool  in  analyzing  complex  petroleum  fractions.  Some  of  these  studies 
were  done  by  Hood  (89),  Melpolder  (90),  Lumpkin  (91,  92).  In  another  approach, 


developed  largely  by  Waterman  and  Associates  (93,  94)  in  Holland  and  Kurtz  and 
Associates  (95)  in  this  country,  analyses  of  petroleum  oils  were  made  by  corre- 
lating physical  properties  of  pure  components  to  the  oil  being  analyzed.  The 
physical  properties  chosen  are  refractive  index,  density,  and  molecular  weight 


chains.  Ali  of  the  information  obtained  with  these  methods  is  in  accord  with 
the  picture  of  lubricating  oils  presented  above.  The  individual  hydrocarbon 
species  and  their  properties  identified  thus  far  have  been  compiled  by  the 
American  Petroleum  Institute  (API)  Research  Projects  6 and  42. 

After  carbon  and  hydrogen,  sulfur  is  the  next  prominent  element  present  in 
base  oils.  The  most  complete  investigation  of  the  sulfur  compounds  in  petroleum 
have  been  done  by  the  API  Research  Project  48,  which  has  isolated  166  sulfur 
compounds  from  four  crude  oils.  The  types  identified  are  alkyl,  cyclic  and 
aromatic  thiols;  alkyl  sulfides,  alkyl-cycloalkyl  sulfides  and  cyclic  sulfides; 
thiophenes  and  aromatic  thiophenes.  The  content  of  thiols  decreases  rapidly 
with  increasing  boiling  point  material.  The  alkyl  sulfides,  as  the  boiling 
point  increases  are  replaced  by  cyclic  sulfides.  From  a mass  spectrometric 
examination  of  a lubricating  oil  fraction.  Lumpkin  and  Johnson  (91)  concluded 
that  the  majority  of  sulfur  compounds  in  the  higher  boiling  fraction  were  of  the 
condensed  aromatic  type.  Later,  Lumpkin  (92)  found  that  the  trinuclear  aromatic 
portion  from  the  347c-460°C  (bSh'-bSC^F)  boiling  range  of  a high-sulfur  crude 
contained  the  following*  indanothiophenes  4.2%;  indenothiophenes,  1.8%; 
dibenzothiophenes,  39.7%;  acenaphthenothiophenes , 1.5%;  and  benzodlthiophenes, 

2.8%. 

I 

The  early  work  of  Bailey  and  Lochte  (1955)  showed  that  the  nitrogen  compounds 
which  were  extractable  with  dilute  mineral  acids,  arbitrarily  classified  as 
"basic",  contain  both  quinolines  and  pyridines  carrying  alkyl  substituents.  In 
1954,  a study  on  the  identity  of  the  "non-basic"  nitrogen  compounds  (cannot  be 
titrated  with  perchloric  acid)  in  petroleum  was  undertaken  by  the  API  Research 
Project  52.  The  results  indicate  that  pyrroles,  indoles,  and  carbazoles  are 
the  principal  constituents.  Phenazines,  benzonitriles,  and  amides  have  been 
detected  in  cracked  oils  The  results  were  summarized  by  Lathan  (97) . 

The  oxygen  compounds  in  petroleum  can  be  divided  into  two  classes:  acidic 

which  can  be  extracted  with  bases,  and  neutral  non-extractable  components.  The 
acids  include  normal  and  branched  chain  acids,  and  acids  containing  a cyclic 
group  in  the  molecule,  e.g.,  2,  3-dimethylcy clohexane  carboxylic  acid.  The 
presence  of  myristic,  palmitic,  stearic,  and  arachidlc  acids  has  also  been 
reported  (98).  Four  neutral  oxygen-containing  compounds,  1-,  2-  or  3-,  and  4- 
methyldibenzof urans  and  4,  6-  dimethyld ibenzof uran,  have  been  Isolated  from  the 
reference  petroleum  of  API  Research  Project  6. 

Excellent  reviews  in  the  compositions  of  petroleum  oils  are  given  by  Mair 
(99),  Bondi  (88). 

b.  Oxidation  Reactions.  The  oxidation  of  the  metal  surfaces  and  the 
lubricants  plays  a prominent  role  in  boundary  lubrication.  Analysis  of  wear 
debris  in  the  boundary  regime  indicates  the  presence  of  metal  oxides  and 
oxidation  products  of  hydrocarbon. 

1.  Oxidation  of  Metals,  Especially  Iron:  Let  the  clean  surface 

of  a meta.1  or  an  alloy  be  exposed  to  oxygen,  chemical  reaction  will  begin  at  the 
metal-gas  interface  If  the  oxides  formed  are  not  volatile,  they  immediately 
form  an  intermediate  layer  between  the  alloy  and  the  gas  The  layers  may  be 
compact  or  interpersed  with  pores  or  cracks,  so  that  further  reaction  may  involve 
atomic  diffusion  or  simply  passage  of  gas  molecules  through  the  pores-  Diffusion 
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mcy  be  through  the  grains,  along  grain  boundaries  or  surfaces  of  the  metal.  When 
an  alloy  reacts,  either  one  or  all  the  component  metals  may  form  compounds  with 
oxygen,  the  resulting  layer  may  be  a conglomerate  of  different  particles  or  clearly 
subdivided  into  individual  layers.  Therefore,  the  characteristics  of  the  oxidation 
layers  for  the  metal-gas  pair  are  vital  in  studying  the  oxidation  of  the  metal. 
Since  the  energy  states  at  an  interface  or  surface  differ  from  those  in  the  bulk 
material,  metallic  compounds  that  are  otherwise  unstable  may  be  formed  at  the 
interface.  Gulbransen  (100)  found  that  Wustite,  FeO  was  normally  unstable  below 
570°C  (1058°F)  but  an  iron  oxidized  below  this  temperature  formed  a thin  film  of 
FiO  underneath  an  Fe304  scale.  The  thickness  of  the  film  became  less  as  the 
temperature  of  formation  was  lowered  but  it  was  still  detected  at  400°C  (762°F). 

The  energy  required  for  four  moles  of  FeO  according  to  the  equation 


Fe304  + Fe  + 4Fe0  (2) 

amounts  to  2500  cal.  at  400°C;  it  was  suggested  that  this  energy  was  supplied  by 
the  surface  of  the  iron. 

If  a metal  has  several  oxidation  states  and  forms  a number  of  compounds  with 
oxygen,  these  compounds  are  arranged  in  layers  with  the  one  richest  in  oxygen  at 
the  oxide-air  interface,  and  the  one  richest  in  metal  in  contact  with  the 
metallic  base.  Paidassi  (101)  has  studied  the  oxide  layers  of  cobalt  and  iron 
and  discusses  the  observed  positions  of  the  various  oxides  formed. 

If  two  metals  of  an  alloy  react  with  oxygen,  the  reaction  products  may  exist 
as  a heterogeneous  mixture  on  the  metal  surface  or  as  separate  layers.  The 
formation  of  solid  solutions  may  also  be  present,  e.g.,  Fe304  - Mn304 , FeO  - MgO, 
Fe304  - MgFe204,  FeO  - MnO. 

Iron  forms  three  stable  oxides:  Wustite,  FeO,  magnetite,  Fe304  and  hematite, 

Fe203.  FeO  is  a p-type  conductor,  has  an  unusually  high  concentration  of  lattice 
defects  which  consist  of  vacant  cation  sites  and  an  equivalent  number  of  electron 
defects  represented  chemically  by  trivalent  Fe  or  Fe^+  ions.  Consequently 
diffusion  is  mainly  cationic  via  vacant  cation  sites.  Magnetite,  Fe304  has  a 
lower  defect  concentration.  Diffusion  is  both  cationic  and  anionic.  Hematite, 

Fe203  is  an  n-type  conductor  in  which  aniens  largely  diffuse.  Their  chemical 
reaction  equations  are  as  follows: 


Fe  + % O2 

- FeO 

[55] 

3Fe  + 2 O2 

Fe304 

[56] 

2Fe  + 3 O2 

2 Fe203 

[57] 

Quinn  (74)  in  studying  the  dry  wear  of  steel  in  air,  determined  the  composition 
of  the  wear  debris  in  a pin  and  disc  system  using  a x-ray  diffraction  technique. 
His  work  shows  clearly  the  volume  percentages  of  a-lron,  Fe203  and  FejC^  for 
externally  induced  temperature  experiments.  The  work  was  conducted  at  a rela- 
tively low  sliding  speed  of  five  centimeters  per  second  and  at  a high  sliding 
speed  cf  625  centimeters  per  second.  From  these  results,  Quinn  demonstrated 
that  at  low  temperature,  the  major  reaction  product  was  Fe203  between  the  two 
rubbing  surfaces  in  air.  As  the  temperature  goes  up,  the  amount  of 
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Fe203  decreases  and  the  amount  of  FejC^  increases.  FeO  is  shown  to  be  absent  in 
this  temperature  region.  The  amount  of  a-iron  i.e.  due  to  mechanical  wear,  is 
low  and  drops  to  an  insignificant  amount  when  the  bulk  temperature  rises  above 
100°C. 


Tao  (102,  103)  proposed  that  all  wear  was  due  to  chemical  corrosion,  i.e., 
the  growth  of  an  oxide  layer  and  its  subsequent  removing  by  rubbing.  He  found 
that  a mathematical  model  based  on  this  assumption  performed  quite  well  in 
correlating  his  wear  test  results  with  a ball-on-cylinder  device  (56  cm/sec). 

The  analytical  results  showed  that  the  rate  of  metal  oxidation  was  the  important 
factor  in  determining  the  wear  rate.  Since  the  oxidation  rate  depends  heavily 
on  the  temperature  at  the  contact  zone,  Tao's  results  also  imply  the 
importance  of  temperature  in  boundary  lubrication. 

2 Oxidation  oi  Lubricating  Oils.  Very  little  is  known  about  the 
oxidation  mechanism  of  lubricating  oils,  largely  because  of  the  extremely  complex 
chemical  compositions  involved.  Voluminous  literature,  however  exists  on 
oxidation  of  simpler  hydrocarbons.  Several  excellent  reviews  on  the  oxidation 
of  lubricating  oils  have  been  written  by  Zuidema  (104)  and  Bondi  (88). 

lubricating  oils  are  composed  of  such  a complex  mixture  of  hydrocarbons  that 
it  is  extremely  difficult  to  identify  specific  compounds  in  their  oxidation 
products,  except  for  such  degradation  products  as  water,  carbon  dioxide,  and  some 
of  the  lower  carboxylic  acids.  Since  lubricating  oils  are  composed  of  mainly 
three  groups  of  hydrocarbons,  paraffin,  naphthene,  and  aromatic,  it  would  be 
appropriate  to  examine  the  nature  of  the  chemical  reactions  of  pure  hydrocarbons 
of  these  three  types  under  oxidation  conditions. 

Chavanne  and  coworkers  (105,  106,  107)  studied  the  oxidation  of  a few 
paraffins  and  naphthenes.  They  oxidized  n -decane,  n-nonane,  and  n-octane  with 
oxygen  at  atmospheric  pressure  and  a temperature  of  120  C,  The  gaseous  oxidation 
products,  which  account  for  10  percent  or  so  of  the  total,  were  similar  in  all 
three  cases:  30-40%  C02;  5-7%  H2;  1-2%  saturated  hydrocarbons.  The  liquid 

products  contained  water,  succinic  acid,  and  formaldehyde.  The  three  hydrocartons 
produced  a significant  amount  of  methyl-octyl,  methy 1-hepty 1 and  methyl-hexyl 
ketones,  respectively,  together  with  series  of  carboxylic  acids  ranging  from 
formic  to  Cn_3  (n  is  the  number  of  carbon  atoms  in  the  hydrocarbon).  Octanol 
was  formed  from  n-octane.  From  these  data,  one  would  conclude  that  the  beta 
carbon  was  attacked  primarily,  the  gamma  secondarily,  and  so  on  toward  the  center 
of  the  molecule.  The  view  was  substantiated  by  other  workers  later. 

Several  investigators  have  shown  that  water  is  one  of  the  principal  oxidation 
products  of  lubricating  oils  (108,  109,  110).  Zuidema  (104)  proposed  the 
following  mechanism  for  the  oxidation  of  paraffins  at  temperatures  between  100 
to  200°C : 

H OOH 

CH  (CH  ) C-CH.  +0,  * CH-(CH-)  C-CH, 

32  n,  3 2 3 2 nt  3 

H | H 

I OH 
* i 

CH,(CH.,)  -c-ch^+o 

3 l n | 3 

H 


CH0(CH,)  COCK-  + H.O 
Jin  3 1 

l°2 

V 

CH_ (CH„)  C00H  = HCH0 
J L n 


126 


Since  Englei  (111)  and  Bakh  (112)  first  proposed  the  formation  of  a peroxide, 
this  scheme  of  oxidation  has  been  called  the  Engler-Bakh  peroxide  theory. 


Chavanne  and  Bode  (105)  made  an  extensive  study  of  the  oxidation  products 
of  1,4-  dimethylcyclohexane  under  atmospheric  pressures  and  at  100°C.  The 
principal  oxidation  products  were  1,4-  dimethylcyc lohexanol , with  water,  carbon 
dioxide,  3~methyl-6-acetylvaleric  acid,  glycolic  acid,  acetic  acid,  B-methyl- 
valeric  ac id,  dimethylcyc lohexanediol , and  acetonylacetone.  Small  amounts  of 
hydrogen,  carbon  monoxide,  methane,  ethane,  and  formic  acid  were  also  detected. 
The  authors  postulated  the  following  mechanism: 
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Stevens  and  Roduta  (113)  made  a study  of  the  slow  oxidation  of  a series  of  benzene 
derivatives.  They  bubbled  oxygen  through  the  hydrocarbon  for  a number  of  days 
at  temperatures  from  SO3  to  1405C.  Their  results  are  summarized  in  Table  46. 
Zuidema  (104)  again  proposed  the  formation  of  hydroperoxide  to  explain  Stevens' 
data.  He  pointed  out  that  oxidation  always  centered  about  a carbon  atom  attached 
to  the  benzene  ring.  The  methylbenzenes  all  formed  aldehydes  with  the  same  number 
of  carbon  atoms  as  the  parent  molecule.  The  substituted  benzenes  containing 
larger  groups  all  formed  ketones  containing  a phenyl  group  plus  the  other  R group 
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originally  attached  to  the  benzene  ring.  Acids,  resulting  from  the  further 
oxidation  of  aldehydes  or  ketones,  were  also  formed.  Acids  were  formed  from  the 
fragments  split  off  in  the  formation  of  ketones  from  secondary  compounds. 

Zuidema  proposed  the  following  mechanism: 
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Larsen,  Thorpe  and  Armfield  (114)  studied  the  rate  of  oxygen  absorption  of  thirty- 
nine  compounds.  The  average  values  for  the  products  of  oxidation  of  five  classes 
of  compounds  are  given  in  Table  47.  Conditions  of  oxidation  were  one  atmosphere 
of  oxygen  and  110°C  (except  naphthalene  compounds  at  150cC) . The  numbers  in 
Table  47  cannot  be  regarded  as  absolute  values,  they  only  indicate  the  average, 
qualitative  results.  Larsen  also  found  that  the  peroxide  content  rose  sharply 
in  the  early  stages  of  oxidation  and  then  fell.  The  tree  acids,  alcohol, 
carbonyl,  water  and  volatile  acid  are  probably  formed  by  decomposition  or 
condensation  reactions.  Aromatics  are  found  to  give  condensation  products 
which  darken  the  oil  and  precipitate,  whereas  paraffins  and  naphthenes  remain 
homogeneous  and  light  in  color  upon  oxidation. 


Larsen  et  al.  (110)  also  determined  functional  groups  in  the  oxidation 
products  of  a number  of  lubricating  oils  under  the  same  conditions,  i.e.,  one  j 
atmosphere  of  oxygen,  and  HO1  or  150^0.  The  results  for  nine  oils  are  shown  in 
Table  48.  The  first  seven  are  lubricating  oil  fractions,  while  the  eighth  is  a 
synthetic  oil  made  by  polymerization  of  cracked-wax  olefins.  The  last  is  a white 
oil.  Both  of  the  last  two  oils  are  aromatic-free.  The  first  six  are  California 
furfural  raffinates  of  various  stages  of  extraction,  as  indicated  by  the 
progressively  increasing  viscosity  index.  It  can  be  seen  that  as  the  aromatics 
are  removed,  there  is  a marked  tendency  toward  a lower  asphaltene  content  in  the 
oxidized  oil.  This  observation  is  consistent  With  the  conclusion  based  upon  work 
on  pure  hydrocarbons  that  aromatics  give  dark  insoluble  condensation  products. 
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Recent  theory  on  the  oxidation  of  hydrocarbons  proposes  free  radical  chain 
reactions.  Excellent  reviews  have  been  made  by  Emanuel  (114,  115,  116),  Boss 
(117),  Agabekov  et  al.  (118),  and  Boudart  (119).  The  radical  chain  mechanism 
is  consistent  with  Zuidema's  hydroperoxide  theory  but  more  sophisticated  in 
concept.  The  chain  mechanism  can  be  roughly  classified  into  four  stages: 
initiation,  propagation,  branching,  termination. 
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RH  + HR 
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RH  + RH  + 0„ 
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Termination 

R + R 

-*  R-R 

[77] 

R + ROO  

- ROOR 

[78] 

RO  + RO  

-+■  ROOR 

[79] 

Equations  [58]  and  [79]  represent  some  of  the  chain  reactions  that  can  take  place 
in  hydrocarbon  oxidation.  Different  starting  chemical  composition,  different 
temperatures  and  pressures  will  make  some  of  the  reactions  dominate  the  system. 

In  summary,  the  oxidation  reactions  of  lubricating  oils  are  numerous  and 
complex.  Peroxides  or  hydroperoxides  play  an  important  role  in  the  early  stage 
of  oxidation.  For  paraffins,  the  main  oxidation  products  are:  alcohols,  alde- 

hydes, ketones,  and  acids.  For  naphthenes,  alcohols,  ketones,  dicarboxylic  acids 
are  formed.  For  aromatics,  phenol,  aromatic  monocarboxylic  acids,  dicarboxylic 
acids,  and  unsaturated  dicarboxylic  acids  are  formed.  Comparatively,  aromatic 
hydrocarbons  are  more  stable  toward  oxidation  than  naphthenes,  which,  in  turn, 
are  more  stable  than  paraffins. 

In  this  section,  the  nature  of  the  oxidation  reactions,  intermediate 
compounds,  final  products  of  the  lubricating  oils  have  been  examined.  The 
possible  chemical  effects  of  these  reaction  products  on  boundary  lubrication 
will  be  discussed  in  a later  section. 

c.  Thermal  Decomposition  Reactions  of  the  Lubricants.  Petroleum  hydro- 
carbons, upon  heating,  decompose  to  form  smaller  or  more  reactive  species. 
Secondary  reactions  from  these  initial  products  may  result  in  the  formation  of 
larger  than,  smaller  than,  or  the  same  size  as  the  molecules  present  in  the 
starting  material.  The  temperature  at  which  this  phenomenon  occurs  depends  on 
the  bond  energies  of  the  weakest  bond  in  the  molecule.  In  this  respect,  the 
carbon-hydrogen  bond  is  stronger  than  the  carbon-carbon  bond.  Recent  studies 
(120,  121)  have  shown  that  triple  bends  are  more  thermally  stable  than  the  double 
bonds  which  are  more  stable  than  the  single  bonds.  Within  the  carbon-carbon 
bonds,  a naphthenic  ring  is  stronger  than  a straight  chain  hydrocarbon,  Among 
the  paraffins,  the  straight  chain  molecules  are  more  stable  than  the  branched 
ones.  In  a given  class  of  hydrocarbons,  the  susceptibility  to  cracking  increases 
with  the  molecular  weight . 

The  breakdown  of  paraffinic  hydrocarbons  occurs  at  any  carbon-carbon  bond 
with  some  predominance  at  the  middle  of  the  chain. 
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The  rupture  of  naphthenic  rings  takes  place  only  under  severe  conditions. 
Temperatures  of  1000' F are  required  to  achieve  extensive  ring  rupture  In 
complex  organic  molecules,  the  decomposition  of  paraffinic  side  chains  attached 
to  naphthenic  rings  constitute  the  major  reactions.  The  alkylated  naphthenes 
crack  to  yield  a paraffin,  naphthene  with  olefinic  side  chain,  or  an  olefin  and 
a naphthene  with  paraffinic  side  chain.  Naphthenes  with  very  short  side  chains 
are  more  stable  thermally  than  those  with  long  side  chains. 
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Olefins,  upon  cracking,  decompose  as  well  as  polymerize  through  condensation 
reactions.  Under  moderately  high  temperature,  both  types  of  reactions  occur 
simultaneously.  Generally,  however,  under  conditions  of  gradually  increasing 
temperatures,  the  polymerization  or  formation  of  high  molecular  weight  hydro- 
carbons would  be  the  first  reaction.  With  slightly  mere  heat,  isomerization  and 
decomposition  reactions  set  in.  The  simple  mono-olefins  are  found  to  polymerize 
to  higher  hydrocarbons  at  300°C  with  sufficient  contact  time.  High  pressure 
enhances  the  polymerization. 

The  thermal  stability  of  aromatic  hydrocarbons  depends  on  the  nature  of  its 
side  chains.  The  benzene  ring  by  itself  is  similar  to  a naphthenic  ring  requiring 
high  temperatures  (1000JF)  to  obtain  substantive  ring  breakage.  High  boiling 
fractions  from  petroleum  contain  arom.  tic  hydrocarbons  having  long  paraffinic  side 
chains  These  side  chains  exhibit  stabilities  similar  to  those  of  paraffinic 
hydrocarbons.  A number  of  studies  (123,  124)  have  been  conducted  on  the 
stability  of  alkylated  aromatic  hydrocarbons.  The  methylated  aromatics  are  the 
most  thermally  stable  compounds  at  the  same  carbon  number  as  compared  to  other 
hydrocarbons . 

The  most  widely  accepted  mechanism  of  the  thermal  decomposition  of  hydro- 
carbons is  the  free  radical  theory.  Rice  and  Herzfeld  (125)  proposed  that  the 
reactions  were  initiated  by  the  generation  of  free  radicals  by  splitting  the 
molecule  at  its  weakest  bond.  One  of  these  radicals  abstracted  the  hydrogen 
atom  from  the  parent  compound  to  form  a small  saturated  molecule  and  a new  free 
radical.  Both  radicals  could  react  this  way.  The  rates  of  removing  the  hydrogen 
atom  ate  1:3  2:10,3  for  primary,  secondary,  and  tertiary  hydrogens  at  600‘C 
respectively  A comprehensive  monograph  has  been  written  on  free  radical 
reactions  by  Steacie  (126). 

The  effects  of  the  thermal  dec ompos it  ion  on  the  rate  of  oxidation  are  not 
clear.  The  two  kinds  of  reaction  usually  occur  simultaneously.  The  generation 
of  free  radicals  by  thermal  degradation  would  certainly  increase  the  amount  of 
free  radical  concentration  The  free  radicals,  being  very  active  chemically, 
would  seem  to  increase  the  oxidation  rate. 

d.  Boundary  Surface  Reactions.  In  hydrodynamic  lubrication,  the 
surfaces  in  relative  motion  are  separated  by  a lubricant  layer  of  considerable 
thickness  and  under  ideal  conditions,  there  is  no  wear  of  the  solid  surface. 

If  the  sliding  speeds  are  low  and  the  loads  are  high,  the  lubricant  layer  breaks 
down  and  the  surfaces  are  separated  by  lubricant  films  of  only  molecular  dimensions. 
Under  these  conditions,  the  friction  is  influenced  by  the  nature  and  the  chemical 
constitutions  of  the  boundary  surface  and  the  lubricant.  The  bulk  viscosity 
plays  little  or  no  part  in  the  frictional  behavior. 
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Since  the  effectiveness  of  the  lubricant  depends  on  the  nature  of  the 
boundary  film,  there  have  been  many  attempts  (127  through  132)  to  determine  the 
chemical  compositions  and  the  physical  characteristics  of  the  film.  Tabor  and 
Willis  (131)  studied  the  lubrication  of  heavily-loaded  sliding  copper  contacts 
(100  kg/mm2  at  contact)  by  a low  viscosity  (10  cs.)  dimethylsilicone  fluid 
containing  0.5  percent  by  weight  of  stearic  acid.  At  temperatures  above  100cC, 
electrical  resistance  measurements  showed  that  the  film  was  nonconducting. 

Infrared  spectroscopic  analysis  showed  the  film  to  consist  of  a high  concentration 
of  copper  stearate  embedded  in  a cross-linked  polysiloxane  matrix.  The  molecular 
structure  of  the  metal  soap  was  very  complex.  They  further  suggested  that  the 
load-bearing  properties  of  the  film  were  dependent  on  the  high  concentration  of 
copper  stearate,  the  polysiloxane  network  simply  served  as  a matrix.  In  the 
absence  of  the  fatty  acid,  the  polymer  was  readily  worn  away  by  the  sliding 
process . 

Materials  such  as  fatty  acids  which  have  a polar  OH  group  at  one  end  of  the 
molecule  have  been  studied  extensively  in  boundary  lubrication.  They  have  been 
considered  to  adsorb  on  the  metal  surfaces  forming  well-defined  "boundary  films" 
which  are  oriented  monomolecular  adsorbed  layers  attached  to  the  metal  via  the 
polar  head  (OH  group).  These  films  are  considered  to  break  up  and  lose  their 
effectiveness  when  their  melting  points  are  exceeded  Recent  researchers  using 
radioactive  tracers,  electron  microscope,  and  x-ray  diffraction  have  shown  the 
main  constituents  of  the  film  to  be  a mixture  of  hydrocarbons,  surfactants  (polar 
molecules),  metals  and  their  oxides.  The  complexity  and  the  possible  inter- 
actions among  these  constituents  of  the  film  have  been  recognized  only  in  the 
last  decade. 

Summets-Smith  (132)  studied  the  chemical  reactions  cn  bearings  in  lubricated 
industrial  machines.  The  temperatures  reached  in  the  bearings  of  high  speed 
machines  are  such  that  chemical  reactions,  at  times  leading  to  bearing  failure, 
can  occur.  A summary  of  the  chemical  reactions  involving  sulphur  is  presented  in 
Table  49.  In  item  1,  the  metals  involved  are  tin-rich  alloys  (typical  composition 
80%  Sn,  15%  Sb,  5%  Cu)  There  are  two  points  worthy  of  comment  in  these 
examples-  The  reaction  products  were  identified  by  x-ray  diffraction  as  a 
mixture  of  a basic  copper  sulphate,  C0SO4  2Cu(OH)2,  and  a copper  sulphide, 
4Cu2S.CuS.  Both  of  these  were  confirmed  as  reaction  products  by  the  deficiency 
of  copper  in  the  underlying  white  metal  and  evidence  of  intergranular  attack. 

The  chemical  reaction  appears  to  be  closely  temperature  dependent.  On  a thrust 
pad,  the  thickness  of  the  reaction  layer  followed  the  temperature  profile  on  the 
pad.  The  second  point  is  that  after  numerous  bearing  failures  with  the  machines 
the  bearing  metal  was  changed  to  a copper  free,  lead-rich  white  metal  (10%  Sn, 

15%  Sb,  75%  Pb) . Nc  failure  occurred  within  the  normal  useful  life  of  the  bearing 
in  the  same  environment 

Godfrey  (128)  reported  that  the  major  constituent  of  the  boundary  lubricating 
film  in  a steel-on-steel  system  was  Fe304  with  small  amount  of  FeS  and  FeO.  The 
tests  were  conducted  with  a SAE  extreme-pressure  machine  and  a sulfurized  mineral 
oil  (sulfur  saturated  DSP  mineral  oil)  under  200  pounds  or  more  loading.  An  iron 
carbide  (77.4%  Fe,  and  17.4%  C)  was  identified  in  a tew  cases.  The  test  cups  were 
made  of  SAE  4615  steel  (60-62  Rockwell  hardness). 

Ever  since  Hermance  and  Egan  (133)  described  the  amorphous  organic  solid, 
they  found  on  the  palladium  contacts  as  "friction  polymer",  the  term  has  been 


used  to  describe  the  resinous-like  wear  debris.  Fein  and  Kreuz  (134)  analyzed 
the  compositions  of  this  amorphous  solid  In  the  contact  region  with  infrared 
spectroscopy  and  found  72.7%  C,  6.6%  H and  4.1%  Fe.  The  solid  analyzed  was 
obtained  from  a four-ball  machine  with  52100  steel  ball  bearings  at  0.35  cm/sec 
sliding  velocity  with  benzene  as  the  lubricant.  This  material  was  insoluble  in 
common  solvents  and  appeared  to  gradually  char  on  heating.  The  infrared  spectrum 
as  well  as  its  chemical  analysis,  showed  that  it  was  organic  in  character.  When 
cyclohexane  vapor  was  used  instead  of  benzene,  the  "polymer"  generated  was 
partially  soluble  in  cyclohexane  and  completely  soluble  in  chloroform.  Both 
"polymers"  indicated  the  existence  of  carbonyl  (C=0)  and  probably  carboxylate 
soap  (COO)  bonds.  Besides  these  resinous  solids,  Fe203  and  FejO^  were  also 
present.  They  also  found  that  if  oxygen  was  excluded  from  the  system,  the  wear 
product  consisted  of  black  ferromagnetic  particles,  identified  as  "FeC"  by  x-ray 
diffraction. 

Oxygen  seems  to  play  an  important  role  in  these  reactions  involving  the 
boundary  surfaces.  Brown  and  Burton  (135)  studied  the  frictional  behavior  of 
52100  steel  as  a function  of  oxygen  partial  pressure  and  found  at  extremely  low 
oxygen  concentration,  1C-9  torr,  the  friction  was  very  high  and  severe  wear 
resulted.  At  high  oxygen  pressure,  the  oxide  formed  acted  as  a boundary  lubri- 
cant, hence  reducing  the  frictioh.  Appeldoorn,  Goldman  and  Tao  (136,  137)  studied 
the  effects  of  water  and  oxygen  on  scuffing  and  found  that  water  and  oxygen  caused 
significant  increases  in  wear.  The  tests  were  conducted  on  a four-ball  wear 
machine  with  52100  steel  at  1200  rpm  and  77°F.  The  authors  explained  the 
observations  with  a corrosive  wear  mechanism  of  forming  and  subsequent  removal 
of  metal  oxides.  The  white  oil  used  was  acid-treated  to  remove  all  polar  and 
aromatic  impurities  They  also  contended  that  air  or  moisture  was  not  necessary 
for  successful  additive  action  as  stated  by  Tingle  (138),  Godfrey  (139,  140), 
and  Barcroft  and  Daniel  (141).  This  apparent  contradiction  is  probably  due  to 
different  operating  conditions  as  will  be  explained  later. 

The  results  of  these  studies  demonstrate  that  the  sliding  behavior  of  metal 
lubricated  by  hydrocarbons  under  boundary  conditions  can  be  related  to  chemical 
reactions  at  the  sliding  surfaces  involving  metal,  hydrocarbon  and  the  available 
oxygen  and  water.  The  following  sections  will  examine  the  organometallic 
chemistry  and  the  effect  of  temperature  and  pressure  on  such  reactions. 

All  chemical  reactions  involve  either  the  breaking  or  forming  of  a bond. 

The  strength  of  a bond  is  directly  proportional  to  the  distance  between  the  two 
molecules.  The  activation  energy  indicates  how  easy  a particular  reaction  would 
occur,  however,  there  ate  very  few  data  available  in  the  literature.  An  indirect 
indicator  of  the  relative  energy  requirement  of  the  formation  or  repture  of  a 
particular  bond  would  be  the  heat  of  formation,  AHj . Table  50  lists  some  AHf 
data  on  some  bonds  at  25  C. 

Ftom  Table  50,  it  can  be  seen  that  the  iron-oxygen  b^r.d  has  a large  negative 
heat  of  formation,  indicating  that  the  iron  is  easily  oxidized.  Carbon-oxygen 
bond,  /hich  has  a AHf  of  -12  kcal^mcle,  is  smaller  than  the  CO-O  bond  forming 
olefinir  peroxide  (AHf  = -50  kcal/mcle).  the  direct  bonding  between  iron  and 
carbon  is  relatively  difficult  but  sulphur  combines  readily  with  iron.  The  low 
energy  requirement  of- the  lron-sul  bond  explains  why  some  sulfur -containing  oils 
exhibit  beneficial  effects  in  boundary  lubrication.  Ftom  these  energy  require- 
ments of  the  various  bends,  organic  acids,  peroxides  (typical  oxidation  products) 


i 
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would  be  relatively  easy  to  react  with  iron.  A hypothetical  RCOO-Fe  bond  would 
require  a minimum  amount  of  energy  to  form.  Ln  this  respect,  the  hydrocarbons 
could  be  oxidized  to  give  peroxides  or  acids  which  would  react  with  the  iron  or 
iron  oxides.  If  oxygen  is  not  available,  the  iron  may  react  with  the  more  polar 
molecules  like  phosphate,  chlorides  or  long  chain  hydrocarbons  at  a higher 
temperature  to  meet  the  higher  energy  requirements. 

There  are  two  types  of  organometallic  compounds:  one  involves  the  direct 

bonding  of  metal -carbon  bond;  the  other  type  contains  both  the  metal  and  carbon 
atoms  but  does  not  have  a direct  carbon-metal  bonding.  The  latter  type  consists 
mainly  metal  alkoxides,  organometallic  complexes  and  organometallic  polymers. 
This  type  is  quite  difficult  to  analyze  and  some  are  very  unstable.  The 
postulated  RCOCM  (M=metal)  bonding  will  be  this  type  of  compounds.  The  vast 
literature  or  organometallic  chemistry  considers  mainly  the  compounds  having  a 
direct  carbon-metal  bonding. 


Carbor.  is  a fairly  electro-negative  element,  hence  might  be  expected  to  form 
ionic  bonds  only  with  the  most  electropositive  elements,  like  sodium  and  potassium, 
and  tc  form  electron-pair  covalent  bonds  with  other  elements,  like  iron,  copper, 
nickle,  and  chromium.  In  boundary  lubrication,  only  metals  of  the  d-block 
transition  elements  in  tne  periodic  table  are  of  interest.  The  properties  of 


these  organometallic  compounds  are  heavily  influenced  by  their  organic  ligands: 
alkyl,  aryl,  carbonyls,  or  hydrides.  Most  of  the  organometallic  compounds  would 
exhibit  the  following  behavior:  rapid  oxidation;  facile  cleavage  by  protlc 

agents;  and  complex  formation. 


There  are  several  relatively  well-known  iron  organometallic  compounds. 
Ferrocene,  (C2H5)2Fe,  and  its  derivatives  can  be  generated  by  heating  cyclo- 
pentadiene  with  iron  to  lOCT'C  ( 5 7 2 c F ) (145)  in  an  inert  atmosphere.  Relatively 
speaking,  ferrocene  has  exceptional  thermal  stability  (300cC)  as  compared  to 
most  other  organometallic  compounds  (100'-  200°C) . It  had  been  used  as  an  anti 
knock  additive  in  automotive  engines  but  was  abandoned  after  it  was  found  that 
the  wear  in  the  piston  tings  increased  by  as  much  as  eight  times  (146)  compared 
to  the  other  anti-knock  additives  such  as  tetraethylead.  Iron  carbonyls, 
Fe(C0)5,  can  be  formed  by  passing  carbon  monoxide  gas  over  heated  iron.  It  is 
an  active  compound  reacting  with  butadiene  to  give  a complex.  It  is  used  as  a 
catalyst  in  carbonylat ion  (process  in  which  carbon  monoxide  is  introduced  into 
a molecule)  The  hydrocarbonyls  of  cobalt  and  iron  are  efficient  homogeneous 
hydrogenation  catalysts.  Also,  iron  carbonyls  can  isomerlze  olefins. 


In  summary,  the  oxganometallic  compounds  can  be  very  complex.  Their 
chemistry  Is  not  well-understood , many  of  them  are  being  tested  as  a catalyst  ln 
hydrogenation  and  polymerization  reactions  of  hydrocarbons.  In  general,  they 
tend  to  react  more  readily  with  unsaturated  hydrocarbons.  Due  to  the  catalytic 
properties  of  the  transitional  metals,  the  reaction  products,  sometimes,  are  not 
predictable.  Whether  these  products  are  formed  in  boundary  lubrication  systems 


is  not  delimte.  A reasonable  guess  would  be  that  two  types  of  organometal lie 
complexes  are  formed  in  boundary  lubrication  depending  on  the  temperature  at 
the  junction  and  the  oxygen  concentration.  If  temperatures  ate  high  and  oxygen 
is  in  shert  supply,  carbon-metal  bonds  may  form,  if  oxygen  is  plentiful,  then  some 
kind  ol  metal  alkoxides  with  hydrocarbon  ligands  may  form  preferentially  due  to 
the  low  energy  requirement . 

3.  Experimental  Technique  Development.  Since  the  chemical  reactions 
in  boundary  lubrication  are  numerous  and  complex,  it  is  not  practical  to  study  all 
of  the  reactions.  In  reality,  chemical  reactions  that  ultimately  will  affect  the 
lubrication  are  those  involving  the  boundary  surface.  The  bulk  fluid  oxidation/ 
degradation  reactions  are  important  only  when  their  reaction  products  are  adsorbed 
onto  the  surface  and  dominant  the  surface  region.  Therefore,  the  chemical  reactions 
between  the  surfaces  and  the  adsorbed  molecules  are  the  key  reactions.  Since  the 
surface  consists  mainly  of  iron  molecules,  the  key  reactions  to  be  considered  can 
be  further  narrowed  to  chemical  reactions  between  iron  and  polar  hydrocarbons, 

Ot  the  many  parameters  affecting  these  reactions,  temperature  is  the  most 
important  variable.  The  rates  of  chemical  reactions  sometimes  vary  exponentially 
with  temperature.  Also  entirely  different  reactions  may  result  if  the  temperature 
is  changed  over  a large  range.  Therefore,  the  first  step  will  be  to  detine  the 
temperature  under  which  the  chemical  reactions  occur. 

The  temperature  at  the  lubricated  junction,  at  present,  is  ill-defined. 

There  have  been  many  attempts  to  experimentally  measure  the  junction  temperature. 

But  as  can  be  seen  from  the  previous  discussion  on  the  hot-spot  temperature  in 
dry  sliding,  large  discrepancies  exist  among  various  calculations  and  experimental 
results.  In  lubricated  contact,  the  picture  is  further  complicated  by  the 
presence  of  lubricant,  which  is  being  pumped  through  the  junction  by  the  relative 
motion  of  the  moving  surfaces.  From  the  standpoint  cf  heat  transfer,  the  overall 
average  temperature  at  the  junction  will  probably  be  lower  than  the  dry  sliding 
case  since  the  frequency  and  contact  ct  the  surface  asperities  are  reduced  by 
the  lubricant,  hence  less  heat  is  generated  by  friction,  Ihe  circulating  lubri- 
cant acts  also  as  a coolant  to  reduce  the  temperacure.  On  the  other  hand,  there 
will  be  an  additional  heat  so-rce  due  to  the  viscous  dissipation  of  heat  due  tc 
shear.  The  critical  tamification  due  to  the  presence  oi  the  lubricant,  however, 
is  the  localization  of  the  hot-spots  in  the  junction.  Petroleum  oils  usually 
have  very  low  thermal  conductivity  compared  t:  metals.  In  the  case  cf  dry  sliding, 
the  heat  generated  at  the  junction  is  quickly  conducted  away.  In  lubricated 
contact,  the  lubricant  acts  like  insulation.  As  asperity-asperity  contact  occurs, 
the  heat  generated  will  be  dissipated  at  a much  slcwer  rate  providing  localized 
hot-spot,  temperatures.  For  the  adsorbed  molecules  at  the  vicinity  of  the  heat 
source,  'he  localized  hot-spot  temperature  provides  the  energy  for  chemical 
react  ions . 


With  these  considerations  in  mind,  and  with  the  aid  of  the  results  of  the 
previous  work  carried  out  at  this  laboratory,  an  experimental  program  was  put 
together  to  estimate  the  surface  temperature  in  the  lubricated  junction 

a.  Review  of  Previous  Work  Boundary  lubrication  studies  using 
activation  analysis  are  discussed  in  Annual  Report  AFML-1R-67-107,  Part  11.  The 
wear  debris  from  a Shell  Four-ball  Wear  Tester  was  separated  inro  three  fractions; 
namely,  organic,  insoluble  metallic,  and  soluble  metallic  debris.  Ihe  method 


of  separating  the  organic  and  organometallic  debris  from  the  metallic  particles 
in  the  sludge  was  developed  using  pyridine  as  the  extracting  solvent.  Activation 
analysis  was  used  to  determine  the  metal  content  of  both  fractions.  By  assuming 
iron  had  the  same  reactivity  as  manganese,  the  amount  of  iron  in  the  sludge  could 
be  determined  by  counting  the  gamma-rays  emitted  from  the  manganese-56  produced 
by  the  thermal  neutrons  reacting  with  maganese-55. 

A continuation  of  studies  of  this  type  included  in  Annual  Report  AFML-TR-67- 
107,  Part  III,  indicates  that  the  assumption  of  equal  reaction  rates  of  iron  and 
manganese  was  probably  not  valid.  Atomic  absorption  spectroscopy  was  used  to 
determine  the  metal- content  of  the  various  fractions  from  wear  debris. 

Additional  studies  using  atomic  absorption  spectroscopy  are  included  in 
Annual  Report  AFML-TR- 74-201 , Part  1.  The  wear  debt  is  was  analyzed  by  washing 
the  four  ball  assembly  with  naphtha  ana  collecting  the  solids  on  a filter  paper. 

The  filter  paper  was  dissolved  in  pyridine  and  the  solution  was  filtered  again 
to  yield  the  organometallic  fraction.  The  metal  particles  left  on  the  second, 
different  filter  paper  were  dissolved  in  strong  hyurochloric  acid  and  the  filtrate 
contained  the  insoluble  metals  in  the  wear  debris. 

Although  these  investigations  were  exploratory  in  nature,  they  demonstrated 
positively  that  successive  solvent  extraction  coupled  with  atomic  absorption 
spectroscopy  was  a viable  way  to  analyze  the  wear  debris  and  gained  insight  Into 
the  mechanism  of  lubrication. 

b.  Experimental  Equipment  and  Operating  Procedures.  Two  Four -Ball  Wear 
Testers  were  used  in  this  study.  One  is  a Shell  Four-ball  Wear  Tester  which  has 
been  described  in  detail  in  previous  reports.  The  majority  of  the  data  was 
obtained  using  the  second,  new  four-ball  wear  tester  which  has  been  modified  as 
proposed  by  Brown  (147)  of  General  Electric  Company.  This  second  wear  tester 
will  be  referred  to  as  the  GE,  Brown  modified  four-ball  tester  in  this  work. 

1.  Four-ball  Wear  Tester  (GE/Brown  Modification).  A diagram 
illustrating  the  essential  working  parts  of  the  wear  tester  is  shown  in  Figure  25. 
It  consists  of  an  oil  cup  fitted  with  clamps  to  hold  three  balls  in  a planar 
configuration,  and  a precision  chuck  to  hold  the  fourth  spinning  ball  on  top  of 
the  other  three.  The  load  on  the  ball  contact  points  is  controlled  by  air  pistons 
which  provide  0.1  to  180  kg.  loads.  Two  interchangeable  pistons  were  used.  The 
low  load  piston  supplied  0.1  to  18  kg.  force,  the  high  load  1 to  180  kg.  force. 

An  air  bearing  system  supports  the  heating  block  (D)  and  ball-cup  (H) . This 
enables  accurate  measurement  of  the  shear  strain  by  a Statham  force  transducer 
with  D.C.  power  supply.  Both  the  transducer  output  and  the  temperature  control 
system  are  continuously  recorded  on  a Ml nneapolis-Honeywe 1 1 strip  chart  recorder. 
The  rotational  speed  of  the  top  ball  is  provided  by  an  eddy  current  motor  connected 
by  belt  drive  to  the  spindle  giving  stepless  variable  speed  control  under  full 
torque  from  30  to  3600  r.p.m.  A speed  reducer  which  supplied  60:1  reduction 
interchanges  with  the  direct  drive.  The  practical  speed  range  of  the  machine  is 
from  2 to  3300  r.p.m.  Heat  for  the  system  is  provided  by  a three-cartridge 
heating  block  which  provides  controlled  temperatures  to  lGGO'F.  The  thermocouple 
leads  and  heater  connections  dip  intc  the  mercury  pools  to  give  lrictlonless 
contacts.  Provision  for  controlled  atmosphere  studies  has  been  built-in  in  the 
system.  As  shown  in  Figure  25,  the  point  P denotes  the  connection  for  controlled 
atmosphere  gas  inlet.  In  this  study,  either  dry  air  or  helium  was  supplied  to 
the  ball-pot  all  the  time  while  running. 


The  wear  tester  is  made  by  Roxana  Machine  Works,  1036  Hansom  Drive,  St. 

Louis  37,  Missouri. 

The  following  procedures  ware  used  in  operating  the  GE/Browr  modified  four- 
ball  wear  tester: 

1.  The  52100  steel  h"  diameter  bail  bearings  (PRL  Batch  15)  were  washed  in 
naphtha  and  acetone  and  then  wiped  clean  with  a lintless  tissue. 

2.  The  oall-pot  assembly  was  rinsed  in  naphtha,  benzene  and  acetone 
successively.  Then  it  was  dissembled  and  each  individual  piece  (ball  pet,  bottom 
piece,  lock  ring,  clamp)  were  wiped  with  a white  lint  less  tissue  paper  wetted  with 
acetone.  If  any  stain  was  detected  on  the  tissue  paper,  the  step  was  repeated 
until  clean, 

3.  Three  clean  bail-bearings  were  then  clamped  into  the  ball-pot  and  locked 
tight . 

4.  The  fourth  ball  was  placed  into  the  ball  chock  which  was  cleaned  as 
described  in  step  2. 

5.  The  spindle  which  was  to  hold  the  ball  chuck  was  cleaned  with  wet  tissue 
paper  with  acetone  until  no  stain  could  be  observed  on  the  white  tissue  paper. 

6.  The  lubricant  to  be  tested  was  poured  into  the  ball-cup  and  its  volume 
recorded . 

7 . The  ball-pot  was  placed  on  the  heater  assembly  and  the  chuck  was  placed 
in  the  spindle  and  tightened  with  drawbar  (G) . 

8.  Temperature  was  set  on  the  Honeywell  strip  chart  recorder.  The  thermo- 
couple cold  well  was  checked  tc  see  if  additional  ice  was  needed.  The  cold  well 
was  maintained  at  0 C with  an  ice-dist ilied  water  mixture. 

9 The  valve  on  the  controlled  atmosphere  inlet  (P)  was  opened.  The  gas 
flow  rate  was  set  at  0.25  liter  per  minute  for  ail  runs. 

10.  The  air  beating  switch  was  turned  on  and  its  pressure  was  adjusted  to 
15  psi  after  loading  the  ball-pot. 

11.  The  load  switch  was  turned  on  and  the  air  pressure  was  adjusted  to 
desired  pressure  to  give  the  proper  load.  (With  the  high  lead  piston,  1 psi  = 

3 kg.;  low  load  piston,  1 psi  = 0.3  kg.). 

12.  The  automatic  timer  was  set  to  desired  running  time,  usually  100 
minutes . 

13.  The  timer,  heater,  and  motor  switches  were  turned  on 

14.  The  rheostat  of  the  heater  was  set  according  to  the  following  table: 
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RHEOSTAT  SETTING  FOR  VARIOUS  TEMPERATURES  ON 
THE  FOUR  BALL  MACHINE 


Temperature,  °F  Setting 


100° 

30 

167° 

35 

200° 

40 

300° 

50 

400° 

55 

500° 

60 

600° 

70 

700° 

80 

15.  The  Honeywell  chart  recorder  switches  and  the  D.C.  current  supply 
switch  were  turned  on 

16.  The  chain  from  the  Statham  transducer  was  hooked  on  the  arm  of  the 
ball-pot . 


17.  When  the  ball-pot  reached  the  set  temperature,  the  motor  was  auto- 
matically turned  on  and  the  timer  started.  At  the  end  of  the  set  time,  the 
motor  would  automatically  be  turned  off. 

18.  With  high  temperature  runs  (500;F  or  up),  the  ball-pot  should  be 
retightened- when  hot  to  prevent  rolling  of  the  stationary  balls. 

19.  At  the  end  of  the  run,  the  hall-pot  was  unloaded.  The  air  bearing 
and  the  switches' turned  off. 

2.  Atomic  Absorption  Spectrophotometer . The  principle  of  atomic 
absorption  spectroscopy  and  its  application  to  analysis  of  elemental  metals  are 
described  in  detail  by  Slavln  (148).  Briefly  stated,  a liquid  sample  is  converted 
into  an  atomic  vapor  by  a hot  flame  (acetylene-air  mixture)  and  irradiated  by  a 
monochromatic  light  source  whose  wave  length  is  set  to  coincide  with  the  emission 
lines  of  the  metal.  The  absorption  of  the  light  by  the  vaporized  sample  is 
correlated  to  the  metal  concentration. 

The  atomic  absorption  spectrophotometer  used  in  this  study  is  a model  303 
from  Perkin-Elmer  Corporation,  Norwalk,  Connecticut.  It  was  used  in  conjunction 
with  the  Perkin-Elmer  automatic  null  recorder  readout  and  a Servo/Riter  II 
Potentiometric  recorder,  model  PS01W64  by  Texas  Instruments,  Inc.,  Houston,  Texas. 

The  operating  procedures  used  in  this  study  are  similar  to  those  presented 
in  the  manual.  Some  modifications  were  made  to  improve  sensitivity  and  precision. 
The  following  procedures  were  used. 

1.  "Intensltron"  source  lamp  from  Perkin-Elmer  for  iron  was  used,  It  has 
approximately  10  percent  higher  sensitivity  than  an  ordinary  cathode  lamp.  The 
lamp  was  Installed  in  the  instrument . 
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2 Fewer’ switch  was  turned  on.  Half  an  hour  was  allowed  tc  warm  up  the 
lamp.  The  lamp  current  was  set  at  30  ma.  The  wave  length  was  set  at  2483. 

3.  The  lamp  position  was  adjusted  both  vertically  and  horizontally  to  give 
maximum  energy  as  indicateu  by  the  energy  meter. 

4 A three-slot  burner  head  was  used  instead  of  the  single-slot  burner  head. 
This  gave  higher  sensitivity.  The  burner  head  was  inspected  for  fouling.  A 
razor  blade  was  used  to  scrap  it  clean. 

5.  The  burner  head  was  adjusted  so  that  the  light  would  pass  above  the 
center  of  the  burner  head  with  a %"  clearance. 


* 


6.  Air  was  turned  on.  Its  pressure  set  at  30  psi.  Initial  rotameter 
reading  should  be  about  9 (at  the  center  of  the  ball). 

7.  Acetylene  was  turned  on  after  the  exhaust  fan  of  the  hood  was  switched 
on.  Initial  rotameter  set  at  9 (.the  top  plastic  ball), 

8.  The  flame  was  ignited  and  distilled  water  was  aspirated  for  15  minutes 
to  warm  up  the  burner  head 

9.  The  aspiration  rate  for  distilled  water  was  adjusted  to  give  8 e.c./min. 
At  this  setting,  maximum  sensitivity  was  achieved  for  low  concentrations. 

10.  The  wavelength  was  adjusted  to  give  maximum  energy  level  as  indicated 
by  the- energy  meter  by  turning  the  fine  control  knob  to  about  2480  t 10. 

11.  The  energy  meter  indicating  needle  should  always  be  inside  the  black 
region.  During  the  course  of  analysis,  if  it  fell  outside  the  region,  the  wave- 
length could  be  adjusted  to  get  it  back  into  the  black  region  without  any 
adverse  effect.  However,  if  it  was  not  adjusted,  then  the  energy  level  would 
continuously  fall  off  giving  erroneous  results. 

12.  The  energy  needle  should  be  kept  at  the  lowest  possible  level  within 
the  black  region  so  that  minimum  gain  could  be  used.  This  wcuid  reduce  noise 
level  significantly. 

13.  Depending  on  the  sample  fluid,  the  iiame  was  adjusted  by  changing  the 
alr/tuel  ratio  to  give  a lean,  hot  flame.  Caution  should  be  taken  not  to  make 
it  too  hot  lest  blow-back  would  occur  damaging  the  burner.  The  following 
rotameter  meter  reading  were  used  m this  study: 

ACETYLENE  AIR  RATIO  FOR  VARIOUS  FLUIDS  ON  THE  PERKIN-ELMER  ATOMIC  ABSORPTION 

SPECTROPHOTOMETER 


Sample  Fluid 
25%  HN03 
50%  HC1 
2:1  pyridine 
oil  mixture 
pyridine 


Air  Rotameter  Reading 
10 
10 


Acetylene  Rotameter  Reading 
9 
9 


, -t 

N 


NOTE:  All  readings  were  taken  at  the  renter  of  ball.  The  acetylene  reading 

was  tor  the  top  plastl:  ball 


minutes  before  shutdown.  This  would  clear  the  burner  chamber. 

I 1 


i 

j 

1 

! 
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15.  After  every  30  hours  of  operation,  the  burner  chamber  should  be  taken 
apart  and  washed  in  warm  water  and  soap  and  worn  parts  replaced.  For  seme  highly 
oxidized  oils  whose  viscosity  is  high,  the  cleaning  procedure  should  be  followed 
after  each  analysis.  Often  times  a widely  fluctuating  base-line  was  due  to  a 
dirty  burner  chamber . m 


3.  Calibration  of  Standards.  A calibration  curve  was  developed  by 
using  a series  of  carefully  prepared  standard  solutions  whose  metal  concentration 
was  known-  These  known  standards  were  aspirated  to  give  different  levels  of 
absorption.  The  level  of  absorption  depended  on  such  variables  as  the  flame 
temperature,  aspiration  rate,  electrical  noise  level,  gas  flowrate  and  recorder 
gain  setting.  It  was  difficult  to  reproduce  consistently  the  same  absorption 
level  of  a standard  repeatedly.  Therefore,  a peak-height  ratio  calibration  was 
developed.  Instead  of  plotting  the  absorption  levels  of  the  standards  versus 
concentrations,  the  ratio  of  the  absorption  level  with  respect  to  the  same 
standard  solution  versus  concentration  was  plotted.  Thus,  the  effects  of  operating 
parameter  variations  were  reduced  to  a minimum.  In  running  unknown  samples,  the 
same  standard  solution  was  used  after  each  unknown  so  that  the  ratio  of  the 
unknown  peak  height  to  the  standard  peak  height  could  be  measured.  A typical 
trace  is  shown  In  Figure  26  and  subsequent  calibration  chart  is  shown  in  Figure  27. 
By  using  the  ratio  of  two  peak  heights,  good  precision  was  obtained.  However, 
a new  calibration  was  needed  each  time  the  machine  was  used. 

c.  Wear  Debris  Analysis.  The  wear  debris  in  this  study  was  separated 
into  four  portions:  1.  Oil  soluble  or ganometallic  iron;  2.  Oil  insoluble  but 

pyridine  soluble  crganometalllc  iron;  3.  Solid  iron  particles  dissolved  in  3.92 
N HNO3;  4.  Iron  oxides  dissolved  In  5 . 06  N HCL  inhibited  with  0.014  M Acridine. 

The  filtration  apparatus  used  in  this  study  was  a 47  mm  diameter  Millipore 
Pyrex  filtering  unit.  The  filter  paper  used  was  a 0.25  micron  porosity  Teflon 
membrane  filter  paper  (Millipore  Fiuoropore)  which  contained  no  non  contaminants 
and  resisted  pyridine  and  acids.  Due  to  its  construction  (the  thin  Teflon 
membrane  was  sandwiched  between  two  polypropylene  grids)  leaks  developed  during 
filtration.  This  was  solved  by  inserting  a silicone  ruober  gasket  between  the 
filter  and  the  filtering  cup.  A clean  25  c.c.  sampling  bottle  was  used  to 
collect  the  filtrate  inside  the  filtering  flask.  A Duo  Seal  vacuum  pump  was  used 
to  supply  the  vacuum. 


1.  Oil  Soluble  Organometallic  Iron  Sample.  The  oil  from  the  ball-pot  after 
a wear  test  was  drained  into  the  filtering  apparatus.  The  oil  sample  collected 
was  diluted  with  double-distilled  pyridine-  The  mix  ratio  was  2:1  pyridine  to 
oil.  The  purpose  of  this  dilution  was  to  reduce  the  viscosity  of  the  oil  so  that 
samples  with  different  degrees  of  oxidation,  1 e.,  different  viscosities,  would 
end  up  approximately  the  same.  On  the  other  hand,  too  much  dilution  decreased 
the  concentration  of  the  metal  in  the  sample.  Since  there  was  very  little  iron 
to  start  with,  dilution  above  feur  would  make  the  iron  undetectable. 


2.  Pyridine  Soluble  Organ .meta 1 1 ic  Iton  Sample.  The  chuck  holding  the  top 
ball  was  unscrewed  from  the  spindle  and  put  into  the  ball-pot.  Fifteen  cubic 
centimeters  of  double-distilled  pyridine  w.  added  to  the  ball-pot  and  heated 
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until  the  temperature  reached  180JF.  This  was  followed  by  scrubbing  the  balls 
with  a rubber  policeman  and  the  contents  poured  into  the  filtering  cup.  The 
hot  pyridine  was  allowed  two  m'  ,utes  on  the  filter  paper  before  the  vacuum  was 
turned  on.  The  filtrate  was  collected  in  a 25  ml.  sample  bottle.  The  ball-pot 
assembly  and  the  chuck  were  then  carefully  rinsed  with  about  250  c.c.  of  pre- 
filtered naphtha  (0.25  y)  onto  the  filter  to  collect  the  solid  wear  debris.  The 
same  filter  paper  that  filtered  the  oil  now  should  have  only  a mixture  of  iron 
chips  and  iron  oxide  particles. 


3.  Iron  Particles  Sample  A method  to  separate  the  iron  oxides  from  the 
iron  particles  generated  during  the  wear  process  has  been  developed  during  the 
course  of  this  study.  Nitric  acid,  an  oxidizing  agent,  has  a very  low  solubility 
for  any  stable  oxides.  In  fact,  concentrated  nitric  oxide  would  not  dissolve  iron. 
Various  concentrations  of  nitric  acid  were  tested  and  a 3:1  mixture  of  concentrated 
nitric  and  distilled  water  (3.92  N)  was  selected.  The  following  data  were 
obtained  when  various  components  were  mixed  with  25  c.c.  3.92  N HNO3  for  two 
minutes  and  the  mixture  filtered  through  a 0.25  Millipore  Solvinert  filter  paper 
and  the  filtrate  analyzed  for  its  iron  content  with  the  Atomic  Absorption 
Spectrophotometer. 


As  can  be  seen  from  Table  51,  the  solution  rate  or  selectivity  between  Fe203 
and  iron  is  about  800  times;  80  times  for  Fe3C>4  and  about  4 times  for  FeO.  As 
discussed  previously,  FeO  usually  is  formed  at  about  570C  (1026°F) . Therefore 
unless  the  junction  temperature  reaches  that  high,  very  little  FeO  should  exist 
in  the  wear  debris.  Therefore  3.92  N nitric  acid  was  used  to  extract  iron  from 
the  wear  debris- 


The  filter  and  the  filter  holder  containing  the  debris  left  after  the 
pyridine  extraction  was  rinsed  with  20  c.c.  acetone  to  remove  entrained  cil  and 
to  wet  the  filter  paper  which  is  hydrophobic.  Ten  cubic  centimeters  of  3.92  N 
nitric  acid  was  measured  into  the  filtering  cup  and  was  allowed  10  minutes  contact 
time  At  the  end  of  ten  minutes,  the  vacuum  line  was  opened  and  the  filtrate 
collected  in  the  sample  bottle  ready  for  analysis 


4 Iron  Oxides  Work  reported  previously  in  AFML-TR- 74-201 , Part  1, 
involved  the  use  of  acridine-inhibited  hydrochloric  acid  to  separate  iron  and 
iron  oxides.  Those  data  are  shown  in  Figure  28.  It  can  be  seen  that  the 
selectivity  in  solution  is  fairly  good,  ranging  from  3 to  16  times  In  this 
work,  5.06  N HC1  (50:50  mixture  or  concentrated  HC1  and  distilled  water) 
inhibited  with  0.0i4  M acridine  was  used  to  dissolve  oxides  left  on  the  filter 


Before  the  addition  of  the  acid,  the  filter  paper  and  the  filtering  cup 
were  rinsed  with  20  c.c.  acetone  so  that  the  nitric  acid  with  its  iron 
concentration  would  not  mix  with  the  HC1  and  affect  the  result.  Ten  cubic 
centimeters  of  inhibited  HC1  was  introduced  cnto  the  filter  paper  and  ten  minutes 
contact  time  was  allowed. 


In  these  procedures,  only  one  filter  paper  was  used  and  the  various  solvents 
were  allowed  to  contact  the  debris  cn  a fixed  rime  basis.  No  claim  will  be  made 
that  all  of  the  iron  or  organomet a 11 1.  iron  extracted  from  the  wear  debris. 

Under  the  same  conditions,  the  iate  cf  solution  should  be  identical.  Therefore 
these  fixed  time  extractions  when  compared  relatively  tc  each  other,  provided  a 
valid  basis  for  comparison. 
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d.  Wear  Scar  Diameter  Determination.  Wear  s^ar  diameters  were  measured 
with  a microscope  containing  a calibrated  scale  in  the  ocular  lens.  The  measure- 
ments were  made  on  the  three  stationary  bails  both  parallel  and  perpendicular  to 
the  striations.  All  readings  were  averaged.  In  case  the  wear  scar  was  highly 
irregular,  either  an  imaginery  square  or  rectangle  would  be  superimposed  such  that 
the  wear  scar  area  outside  the  square  will  be  equal  to  the  unworn  area  inside  the 
square.  Then  the  following  equation  was  u=ed  to  calculate  the  diameter. 

d * 2/ab/r  [83] 

The  equation  can  be  derived  by  equating  Try'  = ab  where  a,b  are  the  length  of  the 
sides  of  a rectangle. 

e.  Precision  and  Consistency  of  Data.  The  wear  debris  analysis  pro- 
cedure has  gone  through  several  stages  of  development  in  this  work.  Accordingly, 
the  data  can  be  divided  into  three  sets.  Within  each  set,  the  data  are  internally 
consistent  Care  must  be  taken  when  comparing  runs  Horn  different  sets. 

One  set  consists  of  forty-one  runs.  These  runs  were  made  at  the  early  stage 
of  the  experimental  program,  some  ditticulties  in  analysis  were  encountered.  The 
operating  parameters  were  not  clearly  defined.  Hence  not  much  significant  data 
were  generated.  A second  set  consists  of  fifty  runs.  These  runs  were  made  in 
the  Shell  Four-Ball  Wear  Tester  with  various  additives.  Except  for  the  first 
couple  runs,  this  set  of  data  represent  a major  improvement  over  the  first  set. 

A third  set  ol  seventy-four  constitute  the  best  set  of  data.  By  this  time, 
the  analysis  procedures  were  firmly  set  and  the  operating  parameters  well-defined. 
Hence,  most  of  the  conclusions  arc  drawn  from  sets  two  and  three.  Set  one  and 
three  used  the  same  machine,  the  G.E. /Brown  Modified  Wear  Tester  while  set  two 
was  run  on  a different  machine,  the  Shell  Four-Ball  Wear  Tester.  Both  machines 
are  similar  in  construction  but  different  in  sophist icat Ion . Ihe  resulting  data 
from  these  two  machines  are  compatible  All  experimental  data  obtained  from  this 
study  together  with  some  pertinent  results  trum  previous  workers  have  been  used 
in  this  study.  Whenever  tnere  was  a change  in  procedure  or  material,  the  change 
is  noted. 

The  precision  and  repeatability  of  the  results  were  checked  by  repeating 
identical  runs  several  times.  The  results  are  shown  in  Iabie  52  This  actually 
represents  the  best  set  of  data  (set  chree;  The  variance  of  data  is  larger  in 
set  two  and  set  one.  Generally  speaking,  the  wear  scar  diameters  are  quite 
reproducible.  In  this  particular  example,  the  average  diameter  can  be  reproduced 
to  0.90  1 0.005  mm.  The  amount  ol  iron  recovered  as  oil  soluble  compounds  is  too 
small  in  this  series  to  be  determined  in  any  cettamty.  The  amount  of  pyridine- 
soluble  iron  compounds  as  determined  by  atomic  absorption  spectroscopy  is  quite 
consistent  In  the  order  of  magnitude.  Run  43  has  a higher  than  average  amount, 

21  pg.  The  standard  deviaticn,  0 , tor  this  set  ot  data  equals  4.387  and  the 
normal  error  limit  on  the  mean  would  be  om  = o//n  * 4 387/Z4  = 2.19,  i.e., 
pyridine  soluble  iron  = 13.5  ’ 2.19  og.  If  Run  21  is  excluded  according  to  the 
Chauvenet's  criterion  (149),  i.e.,  a reading  may  be  rejected  it  the  probability 
of  obtaining  the  particular  deviation  from  the  mean  is  less  than  ^ N,  the 
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standard  deviation  equals  0.816  and  the  average  will  be  11.0  i 0.47  ug  Therefore 
although  the  range  of  the  data  is  large  (10  to  21  ug) , the  normal  deviation 
certainly  is  much  smaller. 


The  amount  of  iron  and  iron  oxides,  due  to  the  manner  they  are  collected, 
is  much  less  consistent.  As  described  previously,  these  are  washed  out  of  the 
ball-pot  and  collected  on  the  filter  paper.  Complete  removal  of  wear  particles 
from  the  ball-pot  assembly  each  time  is  difficult  to  achieve,  if  not  impossible. 
Therefore,  although  these  samples  can  be  measured  with  much  higher  accuracy  by 
atomic  absorption,  they  have  a larger  standard  deviation,  24  tor  iron  and  34.1 
for  iron  oxides.  In  Figure  29,  the  amount  of  iron  extracted  by  nitric  acid  for 
several  series  of  experiments  (same  operating  conditions)  is  plotted  against  the 
reciprocal  of  absolute  temperatures.  The  general  trend  of  the  data  is  quite 
clear  in  spite  of  the  obvious  scatter.  In  Figure  30,  the  amount  of  iron  as  iron 
oxides  is  similarly  plotted,  the  increased  scatter  of  the  data  coincide  with  the 
higher  standard  deviation  calculated. 


4.  Lubricating  Characteristics  of  a Super-Refined  Paraffinic 
Mineral  Oil-  The  base  fluid  used  in  this  study  is  a super -refined  paraffinic 
mineral  oil,  designated  as  MLO  7789  Its  physical  properties  are  presented  in 
Table  53.  It  contains  no  aromatic  ring  structures  in  the  molecule  and  its 
chemical  properties  are  dominated  by  the  paraffinic  portions.  Generally  speaking 
the  paraffinic  side  chains  determine  the  thermal  stability  of  the  oil  even  In 
predominately  naphthenic  or  aromatic  oils. 


The  fluid  was  further  purified  by  percolation  through  a glass  column,  1" 
o.d.  and  5 ft.  long,  filled  with  activated  alumina  This  removed  oxygenated  polar 
impurities  in  the  oil.  Only  one  gallon  cf  oil  was  treated  by  one  cclunu,  Thus, 
base  fluid  with  the  same  level  of  purity  was  obtained.  For  comparison  purpose, 
untreated  oil  (MLO  7789),  treated  oil  (7789  A),  and  the  first  50  c.c.  of  oil 
from  percolation  (7789  B)  were  run  in  the  GE  /Brown  Modified  Four-Bail  Wear 
Tester  under  identical  conditions.  The  results  are  presented  in  Table  54-  As 
can  be  seen  from  the  table,  at  167  F,  wear  increases  with  fluid  purity.  This 
would  indicate  that  the  polar  impurities  removed  by  percolation  act  as  lubricity 
addit'ves.  Therefore  it  is  only  logical  that  purer  fluid  would  give  better 
additive  response  since  there  will  be  less  components  to  compete  for  the  surtace. 
This  point  is  shown  to  be  correct  experimentally  in  this  work.  I e difference 
between  7789  A and  7789  B is  small  indicating  that  there  Is  little  difference 
between  them. 


a Chemical  Composition  of  the  Ball-Bearing  The  half  inch  52-100 
steel  ball-bearings  are  supplied  by  SKF  Industries.  Its  tolerance  is  controlled 
to  0.000025  inches.  The  chemical  composition  of  the  ball  is  shown  In  the 
following  table. 


CHEMICAL  COMPOSITION  OF  THE  52- 100 •STEEL  BALL  BEARING 


Element 


Carbon  . . 

0 95  - 1.10 

Manganese 

• • • 

. 

0 30  - 0.5 

Phosphorus 

• . • 

• . 

0.025  max. 

Sulfur  . . 

. » • 

• 

0-025  max. 

Silicon 

• - , 

, 

0.2  - 0.35 

Chromium  . 

1 2 - 1.5 

Iron  ■ . 

. . . 

. . . 

96.5  - 97  3 

b.  Static  Experiments.  In  order  to  provide  a basis  for  comparison,  a 
series  of  "static"  tests  were  run  as  a function  of  temperature.  The  tests  were 
conducted  in  the  GE/Brown  Modified  Four-Ball  Wear  Tester.  Identical  procedures 
were  used  to  clean  and  assemble  the  ball-pot.  The  air  flow  rates  were  set  the 
same  as  the  wear  tests.  The  speed  was  sec  at  2.9  r.p.m.  (0.108  cm/sec).  In 
short,  exact  procedures  were  used  as  for  the  wear  te=ts  except  that  the  load  was 
set  at  zero  This  was  done  to  simulate  the  exact  environment,  the  metallurgy,  and 
the  mass  and  heat  transfer  characteristics  of  the  system.  Since  there  was  no 
wear,  only  the  oil  and  pyridine  phase  were  analyzed  for  iron  content.  The 
results  are  presented  in  Table  55  and  Figure  3x. 

As  can  be  seen  in  Figure  31,  both  the  amount  of  iron  in  oil  soluble  compounds 
and  in  pyridine  soluble  (oil  insoluble)  compounds  increase  with  temperature.  The 
iron  found  in  the  oil  phase  Is  the  product  of  the  direct  reaction  between  Iron 
and  hydrocarbon.  The  exact  nature  of  these  products  is  net  known.  Oxygen  seems 
to  play  a very  important  role  in  the  formation  of  these  products,  as  illustrated 
in  Table  56.  As  the  temperature  goes  higher,  the  effect  is  more  pronounced.  At 
SOO^F,  the  amount  of  oil  soluble  iron  is  ten  times  higher  than  the  run  with  helium 
and  at  700  F,  it  is  about  200  times  higher.  From  this  evidence,  it  can  be  deduced 
that  the  oil  oxidized  and  formed  the  various  oxidation  products,  e.g.,  ketones, 
alcohols,  carboxylic  acids,  aldhydes.  These  reaction  products  then  react  with 
the  iron  surface  via  adsorption-reaction-desorption  processes.  The  reaction 
products  may  be  metal  scaps,  ionized  metal  complexes,  or  soluble  metal  salts. 

From  the  bond  energy  consideration,  as  discussed  before,  direct  C - Fe  bonding 
requires  a high  activation  energy  The  lowest  energy  path  is  indirect  bonding 
through  oxygen,  e.g.,  RCOO  - Fe. 

The  iron-organic  compounds  extracted  by  pyridine  are  oil  insoluble.  Their 
formation  appears  to  depend  on  temperature,  oxygen  concentration,  chemicals 
present  and  the  nature  of  the  metal.  Then  oxygen  dependence  though  not  as  strong 
as  the  oil  soluble  crganometalllc= , is  quite  significant.  Preliminary  studies 
indicate  that  small  amounts  ot  polar  additives  also  affect  their  formation 
significantly.  A correlation  between  the  amount  of  oil  soluble  iron  and  pyridine 
soluble  iron  under  various  conditions  is  not  present.  It  seems  that  these  two 
classes  of  compounds  may  be  in  some  way  related,  or  may  actually  have  two 
distinctly  different  chemical  reaction  paths  (mechanisms)  In  their  formation. 

It  has  been  known  for  sometime  in  plastic  industries  that  iron  in  the  presence 
of  high  molecular  weight  hydrocarbons  may  react  to  form  oil  insoluble  resinous 
materials.  Due  to  their  high  molecular  weight  and  the  presence  of  extremely 
chemically  complex  mixtures  cf  different  moieties,  isolation  and  analysis  are 
not  practical.  Few  facts  are  known  about  their  structure  or  chemical  nature. 

In  lubrication  literatute,  reports  abound  about  a "friction  polymer"  (151,  152,  153) 
and  Insoluble  "metal  soaps"  (150,  152,  128,  131).  Generally,  they  have  been 
credited  with  the  effectiveness  ot  the  lubrication  system.  Most  recently, 

Zaslavsky  et  at.  (154)  and  Sanln  (155)  studied  the  "friction  polymers"  in  mineral 
oils  with  and  without  additives  and  both  found  iron  reacting  with  the  hydrocarbons 
and  forming  a part  in  the  resulting  pclymers.  Some  of  these  polymers  are  oil 
soluble  and  some  Insoluble.  Therefore  it  can  be  said  that  the  pytldme  soluble 
iron-organic  materials  are  some  kind  of  high  molecular  weight  cimplexes  of  iron 
and  oil  oxidation  products  and  are  similar  to  the  "friction  polymets"  found  in 
the  wear  processes. 
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c Dynamic  Oxidation-Wear  Study 

1.  Slow-Sliding  Wear  Tests.  To  assess  the  effects  of  a fresh  iron 
surface,  a series  of  experiments  was  conducted  under  conditions  similar  to  the 
static  tests  except  a 40  kg  load  was  applied.  To  avert  frictional  energy  build- 
up as  heat  in  the  conjunction,  the  sliding  speed  was  kept  at  2.9  r.p.m.  Two 
fluids  were  used,  the  MLO  7789  and  the  MLO  7789A.  These  are  super-refined  paraffinic 
mineral  oils.  MLO  7789A  is  the  oil  which  had  been  further  treated  by  percolating 
through  an  activated  alumina  column. 

The  results  for  PRL  7789  and  PRL  7789A  are  presented  in  Tables  57  and  58 
respectively.  As  can  be  seen  from  Table  57,  the  wear  increases  with  temperature 
as  indicated  by  the  mean  wear  scar  diameter.  As  wear  increases,  the  total  amount 
of  iron  recovered  also  increases,  but  the  percentage  of  the  iron  in  the  form  of 
organometallics  (oil  soluble  iron  + pyridine  soluble  iron)  increases  drastically 
from  3 to  46  percent.  This  indicates  that  as  temperature  rises,  chemical  wear 
becomes  increasingly  important  and  effective  lubrication  will  depend  more  heavily 
I on  the  nature  of  the  reaction  products  between  the  oil  and  the  metal. 
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In  Figure  32,  the  mean  wear  scar  diameter  is  ploLted  against  the  bulk  fluid 
temperatures.  At  a sliding  speed  of  2.9  r.p.m.  or  1 mm/sec,  very  little  heat 
builds  up  at  the  conjunction,  hence  it  can  be  reasonably  assumed  that  the 
temperature  at  the  junction  is  approximately  the  same  as  the  bulk  fluid.  As  can 
be  seen  from  Figure  32,  the  wear  scar  diameter  undergoes  two  transition  stages. 
There  is  a sudden  increase  in  wear  at  500°F;  it  then  stays  the  same  at  600°F. 
However,  at  700°F,  the  wear  increases  again.  The  7789A  data  are  plotted  on  the 
same  graph  for  comparison.  In  Figure  33,  the  iron  in  the  oil  soluble  and  oil 
insoluble  organometallics  is  plotted  against  temperature.  It  can  be  seen  that 
similar  transition  regions  are  present.  The  oil  soluble  iron  increases  rapidly 
with  temperature.  At  500JF,  there  is  a surge  in  its  formation  and  at  600°F,  it 
begins  to  level  off.  The  oil  insoluble  iron  has  the  same  Leant ition  point  at 
500°F,  but  continues  to  rise  to  700°F.  In  the  wear  scar  diameter  plot,  the  two 
fluids,  MLO  7789  and  7789A,  behave  the  same  in  the  overall  wear  as  shown  in 
Figure  33.  In  Figure  33,  the  amount  of  organic-iron  produced  from  the  two  fluids 
shows  different  behavior.  The  purer  fluid  7789A  (polar  impurities  removed  by 
percolating  through  activated  alumina)  generates  the  oil  insoluble  organic-iron 
compounds  at  a slower  rate  and  in  lower  amounts  at  high  temperatures.  This 
phenomenon  is  much  less  pronounced  in  the  case  of  the  oil  soluble  organic  iron. 
This  indicates  that  the  amount  of  polar  impurities  in  the  fluid,  though  not 
potent  enough  to  modified  the  wear  pattern,  do  exert  considerable  influence  on 
the  generation  of  oil  insoluble  organometallics,  and  to  a lesser  extent,  oil 
soluble  iron.  Therefore,  polar  impurities  must  have  adsorbed  on  the  surface. 

This  effect  demonstrates  that  the  generation  of  oil  insoluble  organic-iron 
compounds  is  more  surface  related. 

The  transition  regions  observed  in  Figure  33  are  similar  to  the  classical 
wear  scar  diameter  versus  load  plot  in  which  the  wear  goes  through  a sudden 
increase  at  certain  loads,  commonly  referred  to  as  the  transition  load.  The 
friction  level  usually  increases  in  this  region.  As  further  load  is  applied, 
depending  on  the  lubricant,  the  wear  either  stays  at  the  same  level  or  increases 
slightly  with  additional  load.  Eventually  as  load  increases  to  a certain  point, 
seizure  and  welding  will  take  place.  This  phenomenon  has  generally  been  explained 
as  follows.  Load  increases  the  contact  pressure  at  the  conjunction.  Higher 
pressure  decreases  the  lubricant  film  to  a minimum  and  more  asperity-asperity 
contact  are  produced,  hence  higher  friction.  Higher  friction  produces  more  heat 
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and  higher  temperatures  result.  At  this  point,  the  so-called  incipient  seizure 
with  recovery  takes  place.  This  is  the  beginning  of  the  transition  region.  The 
classical  theory  proposed  by  Blok,  Bowden  and  Tabor  is  that  the  melting  point  of 
the  lubricating  film  (reaction  product  of  polar  molecule  with  surface  or  pure 
chemical  species  adsorbed  on  the  surfaces)  has  been  reached,  hence  more  contact 
and  higher  temperature  result  from  the  film  breakdown  until  welding  takes  place. 

In  the  case  of  temperature  versus  wear,  a similar  situation  exist  except  in  this 
case,  the  heat  source  is  not  the  dissipation  of  friction. 

As  discussed  previously,  iron  and  iron  oxides  data  are  less  precise  due  to 
analytical  difficulties.  The  results  of  MLO  7789  and  MLO  7789A  are  presented  in 
Figure  34  and  Figure  35,  respectively.  For  both  fluids,  there  seems  to  be  an 
increasing  trend  of  iron  and  iron  oxide  generation  as  temperature  rises.  In  the 
case  of  MLO  7789A,  the  amount  of  iron  particles  is  always  larger  than  the  amount 
of  iron  oxides  indicating  less  effective  lubrication  and  higher  frequency  of  metal- 
metal  contact.  The  wear  scar  data  in  Figure  32  confirm  this  point.  Another 
interesting  observation  is  that  whenever  there  is  a large  increase  in  the  amount 
of  oil  soluble  and/or  oil  insoluble  iron,  there  is  a corresponding  drop  in  the 
amount  of  iron  particles  produced  or  a sudden  increase  in  the  oxides  formation. 

The  effect  is  more  pronounced  in  MLO  7789A.  This  is  reasonable.  MLO  7789A  is 
a purer  fluid  with  many  polar  impurities  removed.  Therefore  the  polar  oxidation 
products  produced  by  the  oil  at  certain  temperatures  have  less  competition  to 
adsorb/reach  the  iron  surface  and  react  with  iron  forming  oil  soluble  and 
insoluble  organo-iron  compounds  which  evidently  serve  as  lubricant  reducing 
metal-metal  contact.  In  MLO  7789,  the  polar  impurities  in  the  fluid  occupy  the 
surface  preferentially.  The  polar  oxidation  products  have  to  compete  with  them 
for  active  metal  surface.  Apparently  at  higher  temperatures,  the  oxidation  products 
from  the  oil  itself  are  able  to  reach  the  surface.  This  may  explain  the  fact 
that  for  MLC  7789,  the  rise  in  oil  soluble/ insoluble  iron  does  not  take  place 
until  500“F  while  for  MLO  7789A,  400°F  is  sufficient. 

2.  Correlation  Between  Static  and  Slow-Sliding  Wear  Test  Data. 

The  difference  between  the  static  data  and  the  slow-sliding  wear  test  data  is 
that  one  has  wear  and  the  other  hasn't.  Whether  the  freshly  exposed  metal  surface 
from  wearing  processes  exerts  a major  influence  on  the  chemical  reaction  products 
and  lubrication  is  the  question.  The  differences,  if  any,  of  the  two  sets  of 
data  could  also  result  from  a higher  temperature  zone  in  the  conjunction  of  the 
wear  tests.  It  will  be  difficult  to  distinguish  between  the  two  factors  from 
the  data. 

In  Figure  36  the  amount  of  oil  soluble  iron  is  plotted  on  a semi-log  scale 
versus  the  reciprical  of  absolute  temperature,  °R . A straight  line  is  drawn 
through  the  data  points  of  the  static  experiments,  indicating  that  the  oil 
soluble  Iron  in  the  static  case  obeys  the  Arrhenius  type  of  relationship  with 
temperatures  The  slow-sliding  wear  test  data  for  the  two  fluids  are  plotted  on 
the  same  graph,  and  as  can  be  seen,  fall  on  the  line  pretty  well.  In  fact,  the 
correlation  between  the  static  and  the  slow-sliding  wear  tests  for  MLO  7789A  is 
excellent.  MLO  7789,  on  the  ether  hand,  does  deviate  more  from  the  line.  This 
may  be  due  to  the  tact  that  MLO  7789  is  chemically  different  from  MLO  7789A. 

When  the  data  is  expressed  as  oil  soluble  iron  per  unit  area  and  plotted  against 
the  reciprocal  of  absolute  temperature  as  in  Figure  37,  a little  more  data 
scattering  is  observed.  The  definite  trend  of  rising  temperature  with  higher 
product  concentration  remains  the  same.  The  increased  data  scatter  may  be  due 
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to  the  inaccuracy  of  estimating  the  surface  area  or  the  chemically  different 
systems  of  MLO  7789  and  MLO  7789A.  In  spite  of  all  these,  and  considering  all 
the  factors  involved,  the  data  are  well  correlated.  This  implies  first  that 
the  static  experiments  are  chemically  compatible  with  the  slow-sliding  experi- 
ments. Secondly,  the  amount  of  reaction  products  is  an  exponential  function  of 
temperature.  Generally,  only  for  a relatively  simple  chemical  reaction  system, 
the  kinetic  rate  data  (initial  rate  of  formation  or  steady  state  rates)  can  be 
correlated  according  to  the  Arrhenius  equation  which  is  usually  expressed  in 
the  form  as  k = Ae~E/RT,  The  fact  that  the  amount  of  iron  (product  concentration  = 
o/t  rdt  for  each  temperature)  measured  as  oil  soluble  or gano-metallrc  iron 
compounds  obeys  the  Arrhenius  equation  suggests  that  these  reactions,  as  a class, 
in  forming  these  products,  have  rate(s)  of  formation  with  a constant  functionality 
with  the  product  concentration  or  the  steady  state  rate(s)  mask  the  rate  changes 
in  the  system  at  each  temperature,  Also,  a continuous  chemical  reaction  system 
is  indicated.  There  may  be  different  reactions,  each  more  predominant  at  different 
temperatures,  but  as  far  as  the  iron  involvement,  it  is  continuous. 

In  the  case  of  the  oil  insoluble  organic  iron  (pyridine  soluble),  similar 
treatment' yields  a little  different  results.  In  Figure  38  the  oil  insoluble 
iron  is  plotted  against  the  reciprocal  of  the  absolute  temperature.  The  raw 
data  indicate  that  MLO  7789  behaves  quite  differently  from  MLO  7789A  and 
apparently  there  are  some  differences  between  the  slow-sliding  runs  and  the 
static  experiments.  The  difference  is  small  for  the  same  fluid.  In  a broader 
view  point,  the  fluids  are  the  same  class  chemically  and  the  general  trend  should 
correlate  fairly  well  with  temperature.  This  is  shown  in  Figure  39,  where  the 
oil  insoluble  iron  per  unit  area  is  used.  The  differences  between  the  fluids 
and  the  load  still  exist,  yet  taken  together,  the  general  trend  is  quite  obvious. 
Straight-line  correlation  with  1/T  in  this  case  is  no  longer  observed.  In  the 
MLO  7789  slow-sliding  case,  the  data  can  be  correlated  with  two  straight-lines 
with  the  abrupt  change  of  slope  at  500°F  as  shown  in  Figure  38.  For  MLO  7789A 
slow-sliding,  the  first  big  incremental  change  also  occurs  at  500°F.  In  the 
static  system,  the  first  big  increase  is  at  400°F.  These  data  indicate  that  there 
may  be  two  different  reaction  regions  in  which  the  rates  of  formation  of  oil 
insoluble  organic-iron  compounds  are  different.  One  possible  theory  is  that  the 
observed  behavior  is  due  to  bulk  oil  concentration  effect.  At  500°F,  the  bulk 
oil  oxidation  reaches  a point  that  the  rates  of  oxidation  increase  drastically. 

This  could  have  two  effects:  the  higher  rates  of  oxidation  increase  the 
production  rates  of  the  active  species  for  the  formation  of  oil  insoluble  organo- 
metallics  involving  iron;  or  the  bulk  oil  concentration  of  polar  oxidation 
products  reaches  a stage  that  the  bulk  oil  chemical  composition  is  sufficiently 
changed  to  provide  higher  production  of  oil  insoluble  organic-iron  compounds. 

The  amount  of  iron  particles  and  iron  oxides  data  are  treated  similarly. 

These  results  are  shown  in  Figures  40  and  41.  The  production  of  iron  particles 
depends  on  the  frequency  of  contact  and  the  effectiveness  of  the  lubrication. 

The  production  of  iron  oxides  will  depend  on  the  amount  of  oxygen  available  to 
the  surface  and  the  temperature.  As  explained  before,  the  collection  of  these 
data  was  difficult  and  the  resulting  data  scatter  make  any  interpretation 
difficult.  Nevertheless,  the  general  trends  are  there.  In  Figure  40,  a steady 
rise  in  the  amount  of  iron  particles  produced  per  unit  wear  area  can  be  seen, 
suggesting  increasing  amount  of  contact  between  the  surfaces.  In  Figure  41, 
there  is  enough  scatter  in  the  data  that  virtually  several  lines  can  be  drawn 
with  equal  validity.  The  lines  drawn  in  Figure  40  are  the  lines  that  best 
correlate  the  MLO  7789  data. 
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3.  Dynamic  Wear  Tests.  Two  series  of  experiments  were  made  at 
moderate  speed,  600  r.p.m.  (linear  velocity  22.8  cm/sec  or  0.75  ft/sec)  on  the 
four-ball  wear  testers.  The  bulk  fluid  temperature  was  maintained  at  167°F. 

The  other  operating  parameters  were  maintained  the  same  as  the  slow-sliding 
experiments  namely:  load  = AO  kg.,  duration  = 100  minutes,  atmosphere  = dry  air. 
At  the  end  of  the  run,  the  wear  debris  were  analyzed.  In  Table  59,  the  results 
for  MLO  7789A  (percolated  through  activated  alumina)  are  shown.  There  is  very 
little  oil  soluble  organic-iron  generated  but  the  oil  insoluble  organic-iron 
generation  is  quite  significant.  This  demonstrates  again  that  the  production  of 
these  oil  insoluble  organometallic^jompounds  are  strongly  surface  related.  The 
amount  of  iron  in  oxide  form  is  greater  than  the  amount  of  iron  particles.  This 
suggests  iron  oxidation  plays  quite  important  role  in  the  lubrication.  It  is 
not  clear  whether  the  surface  oxidizes  first  and  then  the  oxides  are  removed  by 
shear  stresses  or  whether  the  surfaces  produce  small  iron  particles  from  contact 
fatique  and  these  particles  subsequently  oxidize  in  the  conjunction  at  elevated 
temperatures . 

Table  60  shows  the  data  obtained  from  the  untreated  oil,  MLO  7789  in  the 
Shell  four-ball  wear  tester.  Higher  amounts  of  both  oil  soluble  and  oil  insoluble 
organic-iron  are  produced.  The  mean  wear  scar  diameters  appear  to  fluctuate 
from  0.60  mm  to  0.87  mm.  It  was  later  discovered  that  the  particular  wear 
tester,  due  to  its  design  had  different  alignment  when  a different  ball-pot  was 
used.  Runs  K1A  and  K25  used  the  same  ball-pot  while  Runs  K1B  and  K26  had  the 
other  ball-pot.  With  this  in  mind,  the  precision  of  the  wear  scar  diameters  is 
actually  quite  good.  Apparently,  the  other  results  were  not  affected  by  this 
problem. 


d-  Junction  Hot-Spot  Temperature  Estimation  with  Chemical  Conversion 
Data.  When  one  considers  the  overall  chemical  reaction  systems  involved  in 
boundary  lubrication,  one  usually  concludes  that  some  kind  of  simplification  is 
necessary  in  order  that  meaningful  data  can  be  obtained  and  interpreted.  In 
this  work,  without  going  to  an  exhaustive  and  probably  fruitless  analytical 
route  of  separations  and  purifications,  the  numerous  reaction  systems  are  lumped 
into  one  gross  chemical  reaction  system.  The  reactants  in  this  system  are,  on 
the  one  hand,  iron  on  the  metal  surface  and  numerous  polar  hydrocarbon  species 
on  the  other.  Because  of  the  little  known  nature  of  the  chemical  products  and 
the  reaction  mechanisms  involved,  the  usual  reaction  kinetic  analysis  is 
difficult  to  apply  here.  The  difficulties  encountered  in  obtaining  reproducible 
data  add  another  restraint  to  the  task.  Ideally,  if  the  concentration  of  the 
iron  as  or ganometallics  can  be  measured  at  the  initial  stages  of  the  rubbing 
process  in  the  slow-sliding  and  the  dynamic  systems,  a plot  of  lineatized  initial 
rates  of  such  reaction  at  various  temperatures  can  be  used  to  correlate  the  rate 
data  between  the  dynamic  system  and  the  static/slow-sliding  system.  This  proves 
to  be  impractical.  In  the  dynamic  experiment,  at  the  start  of  the  rubbing  process, 
due  to  the  tremendously  high  pressure  (very  small  load  supporting  area),  a lot 
of  wear  is  generated  in  a relatively  short  time  and  the  concentration  of  organo- 
metallics  is  very  low  if  at  all  detectable  in  the  present  experimental  set-up. 

There  are  some  questions  as  to  whether  or  not  this  "wear-in"  period  is  in  the 
chemical  lubrication  regime,  i.e.,  the  lubricant  is  not  functioning  effectively, 
hence  it  may  be  more  of  a physical  rubbing  and  tearing  rather  than  chemical 
reactions  controlling  the  lubricating  mode.  Next,  the  precision  of  data  gathered 
in  the  early  stages  due  to  the  low  values,  is  of  such  quality  that  correlation 
built  on  such  rate  date  is  of  dubious  value.  In  the  static  or  slow-sliding 
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experiments,  at  high  temperatures,  this  initial  rate  concept  also  encounters 
some  problems.  It  takes  about  ten  minutes  for  the  ball-pot  assembly  to  be 
heated  up  to  the  required  temperature  and  as  temperature  rises,  reactions  occur. 
The  same  analytical  problems  present  in  the  dynamic  experiment  are  present  here 
also  except  that  the  theoretical  analysis  of  a transient  temperature  region 
complicates  the  interpretation  unnecessarily. 

1.  Steady-State  Concentration  Correlation.  One  reasonable 
alternative  to  the  initial  rate  method  is  the  correlation  based  on  the  steady 
state  conversion  data.  If  the  rates  of  formation  of  the  products  are  constant 
(zero  order  reaction)  or  if  the  product  concentration  reaches  a steady  state  at 
a fixed  temperature,  then  at  a fixed  time  or  after  enough  time  is  allowed  to 
reach  the  steady  state  concentration,  the  product  concentration  will  be  determined 
by  the  temperature  alone.  To  investigate  the  validity  of  this  assumption,  a 
series  of  experiments  was  performed.  In  Table  61,  data  from  dynamic  experiments 
as  a function  of  time  at  167°F  and  600  r.p.m.  are  presented.  In  Figure  42,  the 
change  of  the  mean  wear  scar  diameter  with  time  is  plotted.  As  can  be  seen  from 
the  graph,  after  30  minutes  the  rate  of  wear  is  essentially  constant.  During 

the  first  half  hour,  the  wear  rate  can  be  approximated  reasonably  well  with  a 
straight  line.  In  Figure  43,  a similar  conclusion  can  be  drawn  with  respect  to 
the  oil  insoluble  iron  concentration.  It  rises  very  quickly  and  levels  off  at  a 
steady  state  concentration.  Therefore,  at  least  for  the  dynamic  experiments, 
the  steady  state  correlation  is  justifiable.  In  the  slow-sliding  or  static 
experiments,  since  the  temperature  at  the  junction  is  constant  throughout  the 
run,  the’ concentration  of  oil  insoluble  iron  has  to  be  the  function  of  the 
particular  temperature  of  that  experiment. 

2.  Surface  Area  Considerations.  In  the  high  temperature  region, 
the  ball-pot  assembly  in  the  static  and  slow-sliding  experiments  was  heated  and 
maintained  at  high  temperature  throughout  the  run.  The  bulk  oil  reacted  with  the 
metal  surfaces  in  the  surrounding.  In  the  dynamic  run,  only  a small  fraction 

of  the  oil  was  subjected  to  high  temperatures  at  the  junction  and  in  terms  of 
surface  area,  only  the  wearing  areas  and  the  wear  metal  particles  existing  at 
the  conjunction  were  involved.  To  account  for  these  differences,  the  oil 
insoluble  iron  has  to  be  expressed  in  Mg  per  unit  area.  For  the  static  experi- 
ments, the  total  active  area  throughout  the  run  is  easily  obtained  by  measuring 
and  summing  the  total  surface  area  of  the  ball-pot  assembly.  In  the  slow- 
sliding  and  dynamic  systems,  the  surface  area  of  the  wear  particles  presents  some 
problems.  In  the  slow-sliding  case,  all  of  the  wear  particLes  generated  are 
subjected  tc  the  same  temperature  at  all  times.  Therefore,  the  total  surface 
area  will  be  the  sum  of  all  the  particles  and  the  ball-pot  assembly.  In  AFML-TR- 
74-201,  Part  I,  typical  particle  size  distribution,  measured  under  a microscope, 
is  presentea.  Surface  area  can  then  be  estimated  from  these  data.  The  detailed 
calculations- of  the  surface  area  are  listed  in  Appendix  I.  In  the  dynamic  case, 
onry  those  wear  particles  which  remained  at  the  conjunction  were  subjected  to 
high  temperature.  The  particles  that  were  removed  from  the  junction  would  have 
a temperature  of  167  F,  the  bulk  oil  temperature.  And  during  a dynamic  wear  test, 
the  wearing  area  increased  with  time,  so  the  amount  of  particles  at  the  junction 
would  also  increase  with  time.  To  simplify  the  situation,  all  wear  particles 
are  assumed  to  be  active  at  all  times.  This  assumption  would  reduce  the  average 
temperature  somewhat.  In  view  of  the  approximate  nature  of  the  correlation,  and 
the  precision  of  other  data,  a more  accurate  analysis  is  not  warranted. 
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3.  Estimation  of  the  Temperature  at  the  Conjunction.  As  discussed 
in  the  previous  sections,  the  amount;  of  oil  insoluble  organometallic  iron  per 
unit  surface  area  is  used  to  correlate  the  dynamic  runs  at  600  r.p.m.  with  the 
static  and  the  slow-sliding  runs.  Table  62  tabulates  all  the  pertinent  data. 

The  correlation  between  the  static  and  the  dynamic  runs  with  MLO  7789A  is 
first  considered.  In  Figure  44,  the  data  from  the  static  runs  is  plotted  against 
temperature.  The  datum  point  at  500°F  is  weighted  less  significantly  based  on 
the  general  correlation  presented  in  Figures  38  and  39.  The  dynamic  test  at 
167°F,  40  kg.  load  and  600  r.p.m.  for  100  minutes  generated  7.76  x 10-3  pg/mm2 
of  oil  insoluble  organometallic  iron.  The  static  test  results  showed  that  for 
the  same  amount  of  reaction  products  to  be  formed  in  100  minutes,  the  temperature 
would  have  to  be  715°F  ± 50°F.  For  the  same  fluid,  but  instead  of  the  static 
runs,  the  slow-sliding  with  2.9  r.p.m.,  40  kg.  load  results  showed  that  the 
temperature  would  have  to  be  660°F  ± 15°F  as  indicated  in  Figure  45.  The 
difference  between  the  two  sets  of  data  could  be  attributed  to  several  factors. 

In' the  slow-sliding  tests,  fresh  iron  surfaces  were  continuously  generated.  The 
fresh  iron  surface  would  have  a much  higher  surface  energy  hence  higher  reactivity. 
Another  factor  would  be  that  in  the  slow-sliding,  although  according  to  theoretical 
calculations  the  temperatures  at  the  junction  would  be  the  same  as  the  bulk  fluid 
temperature,  in  reality,  it  might  be  slightly  higher  than  the  bulk  to  generate 
more  reaction  products.  Lastly,  the  surface  area  of  the  wear  particles  of  the 
slow-sliding  runs  were  estimated  from  approximate  measured  profiles  in  which  all 
particles  were  assumed  spherical.  It  is  conceivable  that  the  assumptions  built- 
in  the  surface  area  calculations  could  produce  some  errors. 

MLO  7789  is  the  untreated  paraffinic  mineral  oil,  which  exhibits  different 
behavior  compared  to  the  activated-alumina-percolated  oil  (MLO  7789A).  In  the 
same  dynamic  run  at  167°F,  600  r.p.m.  and  100  minutes  duration,  it  has  generated 
12.34  x 10” 3 pg/mm2  of  oil  insoluble  organic  iron  and  at  the  same  time  generated 
8.64  x 10" 3 pg/mm2  oil  soluble  organometallic  iron.  Therefore  in  this  case,  both 
the  oil-soluble  and  the  oil  insoluble  organometallic  iron  compounds  can  be  used 
to  correlate  the  temperature.  The  data  on  Figure  46  shows  the  oil-insoluble 
organic  iron  versus  temperature.  For  the  same  amount  of  reaction  products,  the 
temperature  at  the  junction  would  have  to  be  640  ± 25CF  in  the  dynamic  test. 

For  the  oil  soluble  iron.  Figures  36  and  37  can  be  used.  The  result  is  640°F  ± 

50^.  Thus,  for  different  purity  levels  (different  high  temperature  behavior) 
different  correlating  parameters,  the  answers  all  agree  very  well.  The  results 
are  summarized  in  Table  63. 

4.  Consecutive  Dynamic  Test  to  Eliminate  "Run-in"  Effects.  In  the 
dynamic  runs,  for  the  first  30  minutes  the  wear  rate  was  high  as  shown  in 
Figure  42,  then  a constant  wear  rate  set  in.  At  the  initial  wearing  stages,  the 
temperature  at  the  junction  could  get  higher  than  the  steady  state  temperature 
associated  with  the  constant  wear  rate.  To  check  this  effect,  a consecutive  run 
was  performed.  In  the  first  100  minutes  run,  the  wear  scar  reached  certain 
dimensions.  The  ball-pot  assembly  was  then  thoroughly  cleaned  with  naphtha, 
acetone  and  hot  pyridine  (180°F)  without  disturbing  the  locking  positions  of  the 
ball-bearings.  Then  fresh  oil  was  used  and  the  second  100  minutes  run  commenced. 
The  wear  debris  were  analyzed  and  the  results  are  shown  in  Table  64.  It  can  be 
seen  that  the  result  falls  right  in  the  range  of  the  four  first  100  minutes  runs, 
indicating  that  as  far  as  the  temperature  correlation  is  concerned,  the  "wear-in" 
has  very  little  effect. 
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e.  Theoretical  Analysis  of  Junction  Temperatures.  The  frictional  process 
whether  between  lubricated  or  unlubricated  surfaces,  is  always  accompanied  by 
dissipation  of  heat.  This  heating  process  raises  the  temperatures  of  the  sliding 
surfaces  and  the  surrounding  environment.  In  the  last  decade,  a great  deal  of 
research  have  been  done  both  experimentally  and  theoretically  studying  the 
temperatures  at  the  frictional  junction.  From  the  works  of  Earles  (79,  80,  81), 
Poweil  (156,  78),  Furey  (77),  Cameron  (157)  and  others,  it  became  increasingly 
clear  that  there  are  two  temperatures  at  least,  at  the  junction;  0s,  the  sub- 
surface temperatures,  and  0p,  the  flash  temperatures  or  hot-spot  temperatures. 

In  the  lubricated  junction,  there  possibly  exists  another  temperature  ©l,  the 
lubricant  film  temperatures.  Since  these  temperatures  play  an  important  role  in 
boundary  lubrication,  they  are  discussed  below. 

The  flash  temperature  is  the  localized  temperature  of  the  flashes  caused  by 
the  disslp  tion  of  frictional  (and  viscous  if  lubricated)  energy  at  the 
asperities  on  one  or  both  surfaces.  The  measurement  of  these  temperatures,  0f  is 
extremely  difficult  as  they  are  of  only  a few  milli-seconds  duration  and  occur 
over  an  extremely  small  area  relative  to  the  contact  area.  Most  of  the  experi- 
mental investigations  and  theories  are  centered  around  dry  wear  systems.  The 
addition  of  lubricant  into  an  already  complex  system,  at  this  point,  has  not 
been  very  fruitful.  The  sub-surface  temperature  is  the  steady  state  temperature 
existing  just  below  the  asperity  (0.1  mm)  surface.  Usually  the  bulk  of  the  metal 
in  the  system  is  large  enough  that  the  heat  is  treated  as  conducting  away  into  a 
semi-infinite  slab.  The  temperature  at  or  near  the  surface  interface  is  a 
dynamic  steady  state  temperature  established  by  che  heat  sources  (flash 
temperatures)  and  heat  transfer  away  from  the  surface  by  conduction  into  the  bulk 
material.  In  the  case  of  a lubricated  junction,  the  hydrocarbon  lubricant  due 
to  its  relative  high  heat  capacity  and  low  heat  conductivity  will  act  as  an 
insulating  material.  From  the  continuous  generation  of  temperature  flashes,  the 
oil  film  temperature  builds  up  near  the  surfaces  very  fast.  This  temperature 
gradient  in  the  oil  film  has  never  been  experimentally  measured  (the  film  may  be 
only  100  to  500  A thick) , and  it  has  only  been  qualitively  analyzed  from  a 
mathematical  point  of  view  by  Archard  (70) . 

Furey  (77)  in  his  work  studying  the  dry  wear  of  metals,  used  both  the 
dynamic  thermocouple  and  imbedded  thermocouple  methods  to  measure  0s  and  0p. 

The  dynamic  thermocouple  consisted  of  rubbing  the  two  dissimilar  metals,  e.g., 
iron  and  constantan,  against  one  another  and  the  junction  temperatures  measured 
by  the  e.m.f.  generated.  A ball-on-cylinder  device  was  used  in  his  study.  An 
oscilloscope  was  also  employed  to  monitor  the  temperature  flashes.  Before  the 
temperature  data  obtained  in  Furey' s work  are  examined,  some  discussion  of  the 
two  temperature  measuring  devices  (imbedded  thermocouples  and  dynamic  thermo- 
couples) are  in  order.  The  use  of  imbedded  thermocouples  close  to  the  contact 
region  interferes  with  the  heat  flow  and  it  is  doubtful  that  they  indicate  any 
real  temperatures- reached  at  the  tips  of  asperities  in  contact.  They  will  give 
an  estimation  of  the  sub-surface  temperatures.  The  dynamic  thermocouples,  by 
using  dissimilar  metals  and  measuring  the  thermoelectric  voltages  generated  by 
contact,  theoretically  would  give  the  true  temperatures  at  contact.  In  reality, 
only  a time  average  temperature  is  obtained.  The  measured  e.m.f.  is  some  average 
value  pertaining  to  all  the  points  in  contact  at  any  given  instant.  If  all 
points  in  contact  are  at  the  same  temperature,  the  e.m.f.  generated  will  give  a 
true  response.  However,  if  the  points  are  at  different  temperatures,  the  measured 
e.m.f  will  be  a kind  of  average  value  which  is  not  necessarily  even  an 
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arithmetical  average,  it  will  depend  on  the  relative  ratio  of  the  contact 
resistance.  The  temperature  reading  probably  will  approach  the  temperatures  of 
the  contacts  that  have  the  least  electrical  resistance.  Since  the  electrical 
resistance  increases  in  contacts  at  elevated  temperatures,  one  may  expect  to 
read  a temperature  considerably  below  the  mean  temperature.  The  use  of  an 
oscilloscope  only  increases  the  response  time  of  the  equipment  but  the  basic 
problem  remains.  Whether  the  asperities  have  uniform  temperature  across  the 
contact  area  or  not  can  be  resolved  from  the  following  considerations.  The  area 
of  contact  is  often  formed  by  elastic  deformation  of  curved  bodies  in  nominal 
point  contacts  (as  in  4-ball  machine,  pin-on-disc,  ball-on- cylinder  devices). 

In  these  cases,  the  contact  pressure  is  exactly  or  approximately  elliptically 


distributed  over  the  contact  area  according  to  the  classic  equations  of  Hertz. 
Thus,  the  power  density  (rate  of  frictional  heat  production  per  unit  area)  in  a 
sliding  Hertzian  contact  will  be  likewise  elliptically  distributed.  Therefore, 
different  temperatures  will  result  at  areas  under  varying  power  density  distri- 
bution. The  hottest  spot  will  be  at  the  center  of  contact  area.  Also  in  a time 
frame,  at  any  instant,  different  stages  of  contact  will  have  different 
temperatures.  With  these  considerations  in  mind,  Furey's  data  will  be  examined 
below. 


When  the  ball  and  rotating  cylinder  first  come  into  contact  at  the  start 
of  a test,  the  surface  temperature  increases  instantaneously  to  a high  level 
and  remains  essentially  at  this  level  during  a test  at  a given  load  and  speed. 
Although  there  are  variations  in  surface  temperature  with  time,  these  are 
generally  small  and  can  be  explained  on  the  basis  of  frictional  force  change 
during  the  test.  When  the  rotating  cylinder  is  stopped,  the  surface  temperature 
drops  instantaneously  to  the  original  value.  These  temperatures  are  recorded  as 
a recorder  trace  in  the  experiments  using  the  dynamic  thermocouple.  The  same 
temperatures  are  recorded  with  a high-speed  camera  taking  a picture  of  the 
oscilloscope  trace.  It  was  found  that  although  the  average  surface  temperature 
is  about  149"“?,  the  maximum  surface  temperature  is  as  high  as  300  F.  The 
system  temperature  is  77°F.  These  temperatures  are  obtained  from  a load  of 
250  gm  and  56  cm/ sec  linear  velocity.  The  duration  of  these  temperature  flashes 
or  spikes  revealed  in  the  oscilloscope  picture  is  about  a 10  H second.  One 
important  result  of  these  data  is  the  relatively  constant  temperature  with  time. 
Since  wear  is  occurring  throughout  the  test,  this  suggests  that  the  temperature 
rise  is  affected  very  little  by  the  gross  wear  area.  Furey  rechecked  this  point 
by  running  more  experiments.  For  a total  58-fold  increase  in  wear  area,  the 
surface  temperatures  remain  relatively  constant,  124. 5CF  versus  123. 5°F.  This  is 
further  evidence  that  the  surface  area  does  not  appear  to  affect  the  temperature 
at  the  junction  significantly. 


Furey's  work  is  with  dry  sliding.  From  his  results,  one  gains  the  basic 
concept  of  the  nature  of  surface  temperature  flashes.  The  temperatures  measured 
or  theoretically  estimated  are  generally  average  temperatures  over  the  contact 
area.  In  terms  of  chemical  reactions  localized  at  the  asperity  tips,  the 
average  temperature  at  ihe  junction  is  irrelevant.  It  is  the  flash  temperatures 
under  which  the  chemical  reactions  (oxidation/degration/polymerization)  and 
adsorpt lon/desorption  processes  are  initiated  and/or  terminated  that  are  the 
key  parameters  in  the  system. 


In  the  lubricated  contact,  one  additional  heat  transfer  consideration  enters 
the  picture.  The  oil  film  between  the  surfaces  acts  as  an  insulation  When 


heat-energy  is  continuously  fed  into  the  oil  film  inside  which  there  is  a net 
accumulation  of  energy  until  a steady  state  temperature  is  reached  (at  this 
point,  the  net  gain  of  heat  will  be  zero),  the  temperature  gradient  across  the 
film  could  be  several  times  the  surface  temperatures.  Therefore,  in  lubricated 
contacts,  the  chemical  reactions  could  take  place  at  the  surface  where  the 
temperature  flashes  are  occurring  or  inside  the  oil  film  and  the  resultant  polar 
species  move  toward  the  metal  surfaces. 


1.  Theoretical  Hot-Spot  Temperature  Calculations.  Several  theoretical 
studies  have  been  made  on  surface  temperatures  generated  by  friction.  In  1937, 
Blok  (68)  investigated  the  problem  of  stationary  and  moving  source  of  heat  in 
plane  surfaces  with  infinite  heat  sinks.  Jaeger  (69)  in  1942  similarly  studied 
the  problem  and  included  considerable  numerical  detail.  Holm  (158)  in  1948  and 
Bowden  & Tabor  (64)  in  1950  presented  a simplified  version  of  Blok's  theory. 
Archard  (70)  in  1959  simplified  and  reorganized  Blok's  temperature  theory  and 
included  such  considerations  as  elastic  and  plastic  deformation.  This  is  often 
called  the  Blok-Archard  flash  temperature  theory.  The  basic  assumptions 
involved  are  a)  the  heat  is  generated  at  the  area  of  true  contact,  considering 
a single  area  of  contact  as  a plane  heat  source,  and  that  b)  this  heat  is 
conducted  away  into  the  bulk  of  the  rubbing  members.  The  theory  required  the 
determination  of  the  heat  flow  into  each  body  which  must  result  in  the  same 
average  surface  temperature  over  the  contact  area.  The  surface  temperatures  are 
expressed  in  terms  of  the  rate  of  supply  of  heat,  the  size  and  speed  of  the  heat 
source,  and  the  heat  capacities  and  thermal  conductivities  of  the  surface 
materials.  Barber  (159)  in  1967,  and  Francis  (160)  in  1970  both  derived 
equations  to  calculate  the  surface  temperature  with  further  refinements  of  the 
Blok-Archard  theory . However,  the  existence  of  a layer  of  low  thermoconductivity 
material,  i.e.,  the  lubricant  is  not  allowed  for.  Therefore  the  theory  will  only 
predict’ the  temperature  in  a dry  sliding  case. 


There  are  several  equations  which  can  be  used  in  calculating  the  flash 
temperatures  based  on  different  or  varying  degree  of  assumptions.  Fein  (161) 
derived  an  equation  for  the  Shell  four-ball  wear  tester  based  on  the  Blok-Archard 
theory,  Furey  (77)  derived  simplified  equations  using  orchard's  results.  Most 
recently,  in  1975,  Bos  (162)  compared  several  simplified  mathematical  models 
with  his  experimental  work  in  determining  the  temperature  at  the  wear  scars  of 
the  four-ball  apparatus.  In  the  following  sections,  various  equations  for 
calculating  the  surface  temperature  rise  are  discussed. 


a.  Blok-Archard  Equations.  The  model  is  based  on  a protuberance  on  the 
surface  of  body  B which  forms  a circular  area  of  contact  A = TTy2 , which  moves  with 
a velocity  V over  the  flat  surface  of  body  C.  Body  B therefore  receives  heat 
from  a stationary  heat  source  and  body  C from  a moving  heat  source.  Let  Qg  be 
the  quantity  of  heat  supplied  per  second  to  body  B and  be  the  quantity  of  heat 
supplied  per  second  to  body  C.  kg  and  k^  are  the  thermal  conductivities,  Cg  and 
Cc  the  specific  heats  and  pg  and  the  densities  of  the  two  bodies.  The  mean 
rise  in  temperature  over  the  contact  area,  On,  for  a stationary  heat  source  can 
be  calculated  from  the  following  equation: 


It  can  be  seen  that  the  mean  rise  in  temperature  is  directly  proportional  to  the 
rate  of  heat  supplied  and  inversely  proportional  to  the  radius  of  the  contact 
area  and  thermal  conductivity. 


For  a slow  moving  heat  source,  the  contact  time  is  long  enough  for  temperature 
distribution  of  a stationary  contact  to  be  established  in  C,  and  thus 


At  high  speeds,  the  contact  time  is  insufficient  to  establish  a temperature 
distribution  and  the  temperature  rise  is  given  by: 


The  speed  criterion  is  established  by  the  dimensionless  parameter  developed  by 
Jaeger  (69),  L = (Vr)/(2a).  Equation  [85]  applies  when  L < 0.1  and  Equation  [86] 
for  L > 5.  For  intermediate  values  of  L,  Archard  presents  a graphical  procedure. 
In  these  equations,  the  rate  of  heat  supply  is  related  to  the  total  energy  input 
to  the  system  by: 

Q = fWgV/J  = Q + Q [87] 


The  distribution  of  the  energy  between  the  contact  members  in  dry  sliding  is 
assumed  by  Archard  tc  be  divided  as  follows.  For  L < 0.1  (low  speeds),  and  for 
bodies  made- of  the  same  material,  the  heat  supply  is  divided  equally,  i.e., 

Qb  = Qc  = %Q»  At  L > 100  (high  speeds)  essentially  all  the  heat  is  supplied  to 
body  C and  Q = Q^.  For  intermediate  speeds,  the  following  criterion  applies: 


Where  the  temperature  rise  for  each  body  B and  C is  calculated  on  the  assumption 
that  all  of  the  heat  is  supplied  to  it  and  the  mean  temperature  rise  0m 
calculated  from  Equf  tior.  [88], 


For  contacts  whose  size  is  determined  by  the  applied  load,  the  mean 
temperature  rise  is  given  by  the  following  equations: 


With  plastic  deformation  at  low  speeds  (L  < 0.1) 


where 


2ci(ttP  ) 


With  plastic  deformation  at  high  speeds  (L  > 100) 

6m  = 0.435  N [92] 

At  moderately  low  speeds  (0.1  < L < 5), 

0m  = 0.25  N L [93] 

Where  9m  ranges  from  a value  of  about  0.95  at  L = 0.1  to  about  0.5  at 
L = 5.  At  moderate  speeds  (5  < L < 100) 


9m  = 0.435  y N 

[94] 

Y = 1/(1  + 0.87  Lh 

[95] 

Generally,  y ranges  from  a value  of  0.72  at  L = 5 to  0.92  at  L = 100.  Since 
the  maximum  flash  temperature  occurs  when  the  maximum  possible  load  is 
concentrated  at  the  smallest  possible  area,  i.e.,  when  the  total  load  is  borne 
by  plastic  deformation  at  a single  area.  Equations  [89],  [92],  [93],  and  [94] 
can  be  used  to  calculate  the  maximum  flash  temperature  in  dry  sliding.  Using 
these  formulas,  the  average  junction  temperature  without  taking  into  account 
the  existence  of  the  oil  film  for  runs  42  through  45  (40  kg,  600  rpm,  167JF)  is 
found  to  be  260°F  and  the  maximum  temperature,  353'F.  Detailed  calculations  are 
shown  in  Appendix  II. 

Fein  (161)  in  using  Blok's  formula  to  calculate  the  flash  temperature  in 
lubricated  contacts  in  the  four-ball  wear  tester,  found  the  temperatures  thus 
determined  were  too  low  to  account  for  the  observed  transition  temperature  data. 
When  he  arbitrarily  increased  his  calculated  temperatures  by  33  percent,  a 
reasonable  correlation  of  the  data  was  obtained.  Following  his  procedure,  ©m  = 
290° F and  9max  = 4153F. 

b.  Francis'  Equations.  Francis  (160)  derived  an  analytical  expression 
for  the  steady-state  interfacial  temperature  field  in  a sliding  circular 
Hertzian  contact.  He  accounted  for  the  ellipsoidal  distribution  of  the 
frictional  power,  hence  heat  generation  in  the  contact  area.  Blok-Archard ' s 
equations  assume  uniform  heat  generation  and  equal  temperatures  of  the  two 
contacting  bodies.  In  an  apparatus  such  as  the  four-ball  machine,  this  is  not 
the  case.  The  area  of  contact  is  formed  by  elastic/plastic  deformation  of  curved 
bodies  in  nominal  contact  and  the  pressure  is  distributed  el liptically  according 
to  the  classical  equations  of  Hertz.  Thus,  given  that  the  friction  force 
contributed  by  a differential  element  of  area  is  proportional  to  the  normal  load 
on  that  element  the  power  density  (i.e.,  the  rate  of  frictional  heat  production 
per  unit  area)  in  a sliding  Hertzian  contact  will  be  elliptlcaliy  distributed. 

For  the  stationary  surface,  the  average  temperature  rise  is  given  by 
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for  the  moving  surface 


where : 


1.016  Q . 1.047  0.774, 

°V  " ~ T"*  U + B + —TfT] 

TrrkB7  D B ' 


Q = rate  of  heat  generation,  cal/sec 

k = thermal  conductivity  of  the  metal,  cal/cm  cal/cm-°C-sec . 
B = dimensionless  parameter  = rV/a 


For  B = 10,  Equation  [97]  simplifies  to 


[97] 


A0V  = 1.148  (-^-r) 


irrkB‘ 


[98] 


The  average  interfacial  temperature  is  given  by 


A0„ 


1/(A0S  + A0V} 


[99] 


The  maximum  temperature  can  also  be  calculated  by 


A6max  = ~vlk  ^ /(f-996  - 1.091B  Q318  - 0.537B  0,271+  b^)  [100] 


The  detail  derivation  for  the  above  equation  is  given  by  Francis  in 
Reference  (160)  The  calculated  temperatures  from  Francis'  equations  are  higher 
than  those  calculated  using  the  Blok-Archard  equation  by  33  to  38  percent.  Thus 
the  empirical  adjustment  by  Fein  to  increase  the  temperatures  by  33  percent  are 
seemingly  justified  by  Francis'  theory.  Using  the  above  equations,  the  average 
temperatures  at  the  junction  in  the  four-ball  machine  is  determined  to  be  338°F 
and  the  maximum  temperature,  442’F.  The  detail  calculations  are  shown  in 
Appendix  II. 


c.  Bos'  Work.  Bos  (162)  continued  Francis'  work  and  using  a computer, 
obtained  the  solution  for  the  whole  region  of  B.  From  regression  analysis,  he 
obtained  the  following  equation: 


where: 
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A6m  = 


Q 0.281 
kr (B  + 4)^ 


rate  of  heat  generation,  Joules/sec. 
thermal  conductivity,  Joules/m2-sec-rK 
radius  of  the  heat  source,  m 
dimensionless  parameter  = Vr/a 


1 


[101] 
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V = sliding  velocity,  m /sec 

2 

a = thermal  dlffusivlty,  m /sec. 

Which  is  good  for  0 < B < 100.  Bos  further  carried  out  a series  of  experimental 
studies  measuring  the  mean  temperature  at  the  wear  scars  of  the  four-ball  machine. 

The  technique  used  was  an  imbedded  thermocouple  0.5  mm  below  the  wear  scar  and 
0.3  mm  in  diameter.  The  experiments  were  made  with  HVI  oil  containing  one  percent 
sulfur.  In  comparing  the  experimentally  measured  mean  temperatures  with  Francis' 
theoretical  model,  the  measured  temperatures  are  always  much  higher  than  the 
calculated  temperatures.  Figure  47  presents  Bos'  data  for  30  cm/sec  linear  velocity 
in  a four-ball  apparatus  (FBA) . It  can  be  seen  that  at  r = 0.15  mm,  the 
difference  between  the  two  temperatures  is  about  83  percent,  and  at  r = 0.4  mm, 
the  difference  is  300  percent.  Bos  theorizes  that  after  the  start  of  an  FBA 
experiment,  the  distribution  of  frictional  heat  over  the  rotating  and  fixed  balls 
changes  rapidly,  the  heat  flow  into  the  lower  balls  increases  because  the  bulk 
temperature  of  the  upper  ball,  which  is  heated  more  intensively,  rises  more 
rapidly  than  that  of  the  lower  balls.  To  prove  this  hypothesis.  Bos  used  an 
aluminum  oxide  upper  ball  which  is  non-heat  conducting.  The  result  gave  a very 
good  agreement  between  theory  and  measured  temperature. 

A different  explanation  of  the  observed  phenomena  is  also  possible.  There 
is  an  oil  film  existing  between  the  contact  area  and  the  conductivity  of  the  oil 
is  low  enough  that  in  the  physical  system,  the  oil  film  is  equivalent  to  the 
non-conducting  AI2O3  upper  ball.  Furthermore,  the  oil  film  will  absorb  enough 
heat  in  a very  short  while  that  the  resulting  temperature  gradient  across  the 
oil  film  reaches  such  proportions  that  there  will  be  heat  flow  towards  both 
metal  surfaces  and  this,  in  essence,  would  also  create  a situation  similar  to 
the  non-conducting  AI2O3  ball. 

From  Bos'  work,  the  temperatures  at  the  junction  can  be  similarly  adjusted 
and  the  resulting  temperatures  are  tabulated  in  Table  65.  Calculational  details 
are  compiled  in  Appendix  II. 

d.  Discussion  of  the  Juncticn  Temperatures.  From  the  previous  discussions, 
it  is  evident  that  the  temperatures  at  the  lubricated  junction  are  ill-defined 
both  experimentally  and  theoretically.  The  experimental  difficulties  in 
physically  measuring  the  actual  temperatures  are  enormous  and  their  interpretation 
equally  arduous.  The  theories  generally  have  too  many  assumptions  and  none  of 
them  account  for  the  presence  of  the  oil  film  and  the  heat  transfer  character- 
istics of  the  testing  machines.  Furthermore,  the  inability  to  predict  the  values 
of  the  physical  properties  of  the  oil  and  the  physical  condition  and  frequency 
of  contact  make  it  all  but  impossible  to  predict  the  temperatures. 

This  work,  so  far,  has  studied  the  temperature  problem  at  great  length. 

The  nature,  and  the  magnitude  of  the  temperatures  have  been  discussed  in  concert 
with  other  researcher's  work.  In  the  following  sections,  the  validity  of  the 
results  obtained  in  this  study  and  the  effects  of  several  parameters  on  the 
results  will  be  discussed  at  length. 

1.  Validity  of  Present  Result  In  this  work,  chemical  reactions 
which  by  themselves  are  not  clearly  defined,  are  utilized  to  back  calculate  the 
average  hot-spot  temperatures  over  a fixed  period  of  time  in  the  junction. 
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This  method  avoids  the  many  physical  difficulties  in  actually  measuring  the 
temperature.  On  the  other  hand,  this  method  relies  on  several  factors  to  be 
successful:  the  ability  of  the  solvent  to  extract  a fixed  percentage  (100% 

would  be  the  best)  of  the  organometallic  compounds  that  are  formed  during 
lubrication;  the  variation  of  temperatures  in  the  junction  throughout  the 
dynamic  run  remain  small;  and  same  types  of  chemical  reaction  products  are  formed 
in  the  static  experiments  so  that  these  products  have  the  same  solubility  in  the 
extracting  solvent  as  the  dynamic  tests. 

Pyridine  is  the  solvent  used  in  this  study.  From  the  repeat  test  data  shown 
in  Table  53,  the  solution  rate  of  pyridine  seems  quite  constant.  Also  from  the 
temperature  organometallic  formation  relationship,  it  is  quite  obvious  that  at 
high  temperatures  much  more  organic-iron  compound  is  formed.  The  exponential 
nature  of  the  formation  rate  of  organometallic  compounds  with  temperature  makes 
the  constancy  of  pyridine  solution  rate  less  critical.  The  existence  of  these 
oil  insoluble  organic-iron  materials  is  supported  in  part  in  that  correspondent 
amounts  are  found  in  the  oil  phase  as  oil  soluble  compounds. 

During  a dynamic  run,  the  pressure  at  the  junction  varies  considerably. 
Initially,  the  load  (16200  gm  in  40  kg.  load)  is  supported  entirely  on  a single 
point  and  the  mean  specific  pressure  is  about  300,000  psi.  At  the  conclusion  of 
the  test,  for  a wear  scar  diameter  of  0.9  mm,  the  same  load  produces  about 

30.000  psi  average  pressure  at  the  junction.  This  is  a tenfold  decrease. 
Traditionally,  people  when  considering  the  temperatures  at  the  junction, 
logically  assume  that  the  temperatures  will  be  proportional  to  load,  i.e.,  large 
change  in  junction  temperatures  during  wear  tests.  Two  considerations  point  to 
the  fact  that  there  may  be  some  temperature  change  during  wear  experiments,  but 
the  degree  of  change  is  quite  small.  Theoretically,  it  must  be  pointed  out 
that  the  tenfold  decrease  in  pressure  relates  only  to  the  average  or  mean 
specific  pressure  at  the  junction.  In  reality,  the  pressure  or  the  load  is  not 
uniformly  distributed  but  ellipsoidal  in  nature  as  shown  in  Francis'  work  (160). 
This  implies  that  even  though  the  average  pressure  is  changing,  the  maximum 
pressure  at  the  center  of  the  junction  may  change  very  little.  Also,  the  wear-in 
process  is  very  quick  in  this  study,  within  a few  minutes  the  wear  scar  is  already 
in  the  0. 5-0.6  mm  range.  Experimentally,  the  consecutive  run  data  presented  in 
Table  64  certainly  implies  the  existence  of  constant  temperatures  at  the  junction. 
The  run  was  started  after  the  normal  100  minutes  of  a wear  experiment  in  which 
the  wear  scar  diameter  had  reached  0.9  mm.  The  whole  ball-pot  assembly  was 
carefully  washed  and  run  for  another  100  minutes.  The  change  in  wear  scar 
diameter,  this  way,  is  drastically  reduced  (from  0.9  mm  to  0.995),  and  the 
corresponding  average  pressure  change  is  from  30,000  psi  to  28,000  psi,  a change 
of  seven  percent.  The  amount  of  organometallic  iron  determined  from  this  run, 
however,  is  about  the  same  as  from  the  first  100  minute  run.  Another  supporting 
evidence  is  offered  by  Furey  (77)  in  his  work  measuring  surface  temperatures  in 

a ball-on-cylinder  device  (the  load  is  also  ellipsoidal  distributed).  He  found 
that  for  a 58  fold  variation  of  wear  area,  the  temperatures  (average  temperatures) 
changed  only  lO’F.  This  Is  shown  in  Figure  47  It  can  be  seen  that  the  result 
is  quite  impressive.  Another  consideration  is  that  the  load  is  supported  by  the 
contacting  surface.  The  real  contact  area  may  be  only  a fraction  of  the  gross 
wear  area.  For  an  asperity  which  is  supporting  the  load,  the  local  pressure  at 
that  instant  could  be  just  as  high  as  the  very  initial  point  of  contact.  One 
interesting  point  is  that  the  flow  pressure  of  the  bearing  is  approximately 

100.000  psi.  The  asperity  which  is  load  catrying  when  subjected  to  any  pressure 


7 4 A % m * M. 


1 


greater  or  equal  to  100,000  psl  will  deform  plastically.  Below  this  pressure, 
the  deformation  Is  elastic.  Theoretically,  at  least,  elastohydrodynamic 
lubrication  should  produce  little  wear  if  any  at  all.  From  the  average  junction 
pressure  consideration,  when  the  wear  scar  diameter  reaches  0.53  mm,  the  average 
pressure  is  about  100,000  psi  and  the  wear  rate  after  this  should  drop  quite 
significantly.  In  Figure  43,  where  the  wear  scar  diameter  versus  time  is  plotted, 
it  can  be  seen  that  there  is  no  sudden  break  in  the  wear  curve  and  the  wear 
increases  at  a constant  rate.  This  indicates  that  from  the  previous  considerations, 
such  as  the  ellipsoidal  pressure  distribution,  small  actual  real  contact  area  and 
the  experimental  evidence  of  this  work  and  that  of  Furey's,  the  temperatures  at 
the  junction  of  the  4-ball  apparatus  in  the  dynamic  runs  in  this  work  change  little. 

To  ensure  that  the  same  type  of  chemical  reaction  products  are  formed  during 
the  static  experiments,  the  same  apparatus  and  the  same  conditions  are  used  so 
that  the  heat  and  mass  transfer  characteristics  between  the  two  sets  of  experiments 
are  identical.  Since  impurities  are  excluded  from  the  oil,  the  same  reactants 
under  similar  - condit ions  would  lead  to  similar  type  of  products 

The  temperature  determined  in  this  study  is  the  average  mean  temperature  of 
the  aggregate  of  hot  flashes  and  mean  surface  temperatures.  It  is  the  temperature 
under  which  chemical  reactions  take  place.  If  the  temperature  varies  during  a 
run,  the  temperature  determined  will  be  the  Integral  average  temperature  under 
which  the  integral- of  the  rate  of  reaction  gives  the  total  amount  of  reaction 
products.  The  theoretical  calculated  temperatures,  on  the  other  hand,  are  all 
equilibrium  or  steady  state  average  temperatures  at  the  junction.  These  two 
temperatures- can  be  different  and  the  fact  that  they  agree  (though  only  after 
the  correction  factor  offered  by  Bos's  measurements)  lends  support  to  one  another's 
validity. 


2.  Effect  of  Surface  Activity  on  Chemical  Reactions.  It  on  has  been 
known  as  a dehydrogenation  catalyst  under  relative  high  temperature  foi  a long 
time.  The  recent  studies  of  Morecroft  (150)  and  Buckley  (163,  164)  indicate 
that  fresh  Iron  surface  produced  by  filament  evaporation  and  condensation  show 
catalytic  activity  when  in  contact  with  otganic  vapors.  Whether  the  catalytic 
effect- is  important  or  nor  in  boundary  lubrication  has  long  been  a subject  of 
conjecture  and  hypothesis  In  this  work,  sets  of  experiments  are  designed  to 
give  some  indications  in  this  regard.  The  difference  between  the  static  experi- 
ments and  the  slow-sliding  runs  is  the  creation  of  wear  or  ftesh  iron  surface 
throughout  the  run.  The  difference  obtained  in  terms  of  the  temperature  deter- 
mination byconversion  data  is  small,  about  55'F  which  lies  within  the  precision 
range  of  the  data.  This  indicates  that  catalytic  effect  may  be  present  but  is 
not  significant.  Further,  in  run  K34  and  108,  excellent  contrast  can  be  observed 
in  this  respect.  The  data  are  presented  in  Table  6b.  For  an  increase  in  surface 
area  of  1800  percent,  the  total  chemical  reaction  product  changes  58  percent  and 
the  oil  insoluble  portion  increases  by  45  percent  only  If  catalytic  effects 
are  important  in  chemical  reactions  in  boundary  lubrication  regime,  the  change 
in  the  integral  conversion  data  should  be  much  more  pronounced.  It  appears  that 
the  iron  surface,  whether  fresh  or  oxide  covered,  reacts  directly  with  polar 
components  rather  than  serving  as  a catalyst.  In  this  respect,  the  reactivity 
between  the  two  surfaces  is  about  the  same,  while  ftesh  iron  surface  may  exhibit 
higher  reactivity  at  high  tempera rur es . Data  from  stop-and-go  intermittent  wear 
experiments  included  in  AFML-TR-74-201,  Part  I,  also  suggest  the  above  conclusion. 


3.  Effect  of  Bulk  Fluid  Concentration  of  Polar  Species  on  Chemical 
Reactions . In  the  dynamic  experiment,  167°F,  600  r.p.m.,  the  bulk  fluid  temperature 
is  low  enough  that  very  little  oil  soluble  iron  compounds  are  produced  indicating 
that  bulk  oil  oxidation  is  insignificant.  Preliminary  studies  indicate  that,  in 
the  same  experiment,  if  a five  percent  octadecyl  alcohol  is  added,  the  amount  of 
oil  soluble  iron  compounds  and  to  a less  extent,  oil  insoluble  iron  are  increased 
drastically.  In  both  the  static  and  slow-sliding  experiments,  the  bulk  oil  is 
subjected  to  constant  high  temperature  at  run  temperatures  of  600°  and  700°F. 

The  bulk  oil  oxidation  is  very  significant  producing  polar  species  such  as 
alcohols,  ketones,  etc.  Whether  this  increased  polar  species  concentration  in 
the  bulk  oil  will  affect  the  resulting  chemical  reactions  and  what  extent  its 
influence  has  on  the  results  is  an  important  question. 

First  let  the  differences  between  MLO  7789  and  MLO  7789A  be  examined  in  this 
regard.  Two  fluids,  having  quite  different  polar  material  concentration  as 
evidenced  by  different  amount  of  organometallic  products  at  different  temperatures, 
upon  correlating,  give  the  same  junction  temperature.  This  would  indicate  that 
within  a single  chemical  system,  this  concentration  effect  is  insignificant. 

When  foreign  materials  which  by  themselves  may  initiate  different  reactions,  are 
introduced,  the  effects  are  dramatic.  At  this,  if  the  same  chemical  system  is 
used  for  correlation,  the  effect  may  again  be  small.  It  can  alsc  be  argued  that 
the  concentration  of  these  polar  species  in  the  "valleys"  of  the  lubricated 
junction  may  be  just  as  high  as  the  high  temperature  static  runs.  Furthermore, 
even  if  the  concentration  effect  of  oxidation  products  at  high  temperatures  is 
significant,  it  would  mean  that  for  the  dynamic  test,  where  low  concentration 
is  observed,  the  temperature  at  the  junction  must  be  much  higher  than  the  present 
results  in  order  to  compensate  for  the  concentration  difference. 


Table  44 


RANGE  OF  SURFACE  ROUGHNESS  OBTAINED  WITH 
CERTAIN  FINISHING  OPERATIONS 


Finish 

- - - - pin. 
h 

max 

h ** 

rms 

Very  rough  surfaces  (brick,  castings,  etc,) 

>40,000 

Rough  machining 

4 5 

10  to  10 

3 

>10 

Medium  machining  (furrows  visible-rough  grind) 

2,000-10,000 

>400 

Fine  machining  (furrows  invisible-fine  grind) 

200-400 

30-100 

Very  fine  grind  (turned) 

20-200 

5-30 

Honed 

40-80 

2-10 

Lapped,  polished 

2-40 

<4 

Fine  metallographic  polish 

<10 

<2 

Cleavage  plane 

<1 

**  hmg  - Root  Mean  Square  of  the  height  of  the  asperity 
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Table  4 


PROPERTIES  OF  THE  OXIDES  OF  IRON* 


Properties 

FeO 

Fe2°3 

Fe3°4 

Formula  Weight 

71.9 

159.7 

231.5 

Heat  of  Formation,  kcal/mole 

-63.7 

-196.5 

-267.0 

Gibb's  Free  Energy,  kcal/mole 

-58.4 

-177.1 

-242.4 

Color 

black 

red/black 

black 

Min.  Temperature  at  which  formed. 

°C 

570°C** 

200°C 

200*C 

Min.  Temperature  at  which  formed. 

°F 

1058' F** 

392°F 

392°F 

**  Bernard,  J.,  L'oxydation  de  metaux, 

Gauthier 

Villars,  Paris,  1964 

* Lange's  Handbook  of  Chemistry, 

10th 

Edition, 

McGraw  Hill,  1966 
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Fi-are  24.  TYPICAL  STRUCTURES  IN  LUBE  OIL 


Table  50 

HEAT  OF  FORMATION  DATA  ON  BONDS  INVOLVING  C,  02,  H,  Fe,  S 


ds  AH^.  kcal/mole* 


>C0-0 

-50.5 

Fe-0 

-39.0** 

O-H 

-27.0 

Fe-S 

-22.8** 

>C0-H 

-13.9 

+ 1.4** 


* Data  taken  from  Benson  (98) 

**  Estimated  from  AH^  of  F^O^  = “ 196.3;  FeS  = -22.8;  Fe^C  + 5.4 
from  Kubaschewski  (100)  as  suggested  by  Mortimer  (100; 

> Indicates  double  bond  C=C 


8BE3t «.  mSm ‘Mi  3ha jjS ^ ■>- 


Figure  26.  A TYPICAL  ATOMIC  ABSORPTION  TRACE. 


figure  27.  A TYPICAL  ATOMIC  ABSORPTION  CALIBRATION  CHART 


Table  51 


PERCENTAGE  IRON  AND  ITS  OXIDES  DISSOLVED  IN  3.92N  HN03 


Substance 

Sample  Weight,  gm 

Iron  in  Filtrate,  gm 

Dissolved,  % 

Fe 

0.02863 

0.02360 

83 

FeO 

0,68017 

0.12780 

19 

Fe2°3 

0.58606 

0.00003 

0.005 

Fe3°4 

0.71947 

0.00300 

0.038 

DISSOLUTION  RATES  FOR  IRON  AND  IRON  OXIDE  IN  HYDROCHLORIC  ACID  INHIBITED 
WITH  ACRIDINE , 


standard  deviation  0 4-39  24.0  34,1 


83%  of  Data  Included 
in  the  Band 


Test  Conditions 


Speed  = 2.9  rpm 

Load  = 40  kg 


I1L0  7789A  w/  He  Atm 


□ MLO  7789A 


IRON  PARTICLES  GENERATION  VS.  THE  RECIPROCAL 
OF  ABSOLUTE  TEMPERATURE. 


Test  Conditions: 

Speed  = 2,9  rpm 
Load  = 40  kg 

□ 

A MLO  7789 
□ MLO  7789A 

O MLO  7789A  w/  Additive/or  He  Atm. 


Table  53 


PROPERTIES  OF  SUPER-REFINED  PARAFFINIC  MINERAL  OIL 


Fluid  Type,  MLO  Designation 


Viscosity  at  210=F.,  Cs. 

Viscosity  at  100°F.,  Cs. 

Viscosity  Index  (ASTM) 

ASTM  Slope 

API  Gravity,  Degrees 
20 

Density,  D^  g./ml. 

20 

Refractive  Index,  n^ 

Calculated  Molecular  Weight 

Hydrocarbon  Analysis: 

Wt.  % C as  Aromatic  Rings 

Wt.  % C as  Naphthenic  Rings 

Wt.  % C as  Paraffinic  Chains 


MLO  7789 


0.763 


0.838 


1.4637 


Boiling  Point  Range  at 


600-770 


MICROCOPY  RESOLUTION  TEST  CHART 

NATIONAL  BUREAU  OF  STANDARDS -196  3-^ 


ORGANOMETALLIC  FORMATION  IN  STATIC  TESTS  AS  A 
FUNCTION  OF  TEMPERATURES 


Test  Fluid:  MLO  7789A 

Test  Machine:  GE/Brown  Modified  Four-ball  Wear  Tester 

Test  Conditions:  Load  = 0,  Speed  = 2.9  r.p.m.,  Duration 

Atmosphere  = dry  air  at  0.25  H/min. 


Run  No. 

Temperature,  °F 

- - - Organometallic,  pg  - 
Oil  Solu.  Pyridine 

Solu. 

105 

400 

1 

7 

106 

500 

8 

20 

107 

63 

20 

108 

244 

58 

ORGANOMETALLIC  FORMATION  IN  STATIC  TESTS  AS  FUNCTION  OF  TEMPERATURE 


Table  56 

EFFECT  OF  OXYGEN  ON  THE  FORMATION  OF  ORGANOMETALLIC  COMPOUNDS 
Test  Fluid:  MLO  7789A 

Test  Machine:  GE/Brown  Modified  Four-ball  Wear  Tester 
Test  Conditions:  Speed  = 2.9  r.p-m.,  Load  = 40  kg.,  Duration  = 100  min. 
Gas  flowrate  = 0.25  l/min 


- - - - -Test  Parameters  - 
Run  No.  Temp.  °F 


- - -Organometallic , yg  - - - 
Oil  Solu.  Pyridine  Solu. 


Atm. 


Oil  Soluble 


Temperature,  °F 


Figure  34.  IRON  AND  IRON  OXIDES  FORMATION  PER  UNIT  AREA 
VS.  TEMPERATURE  IN  SLOW-SLIDING  WEAR  TESTS 
(MLO  7789). 


187 


▲ Iron  in  Oxides 


Temperature 


Figure  35 


IRON  AND  IRON  OXIDES  FORMATION  PER  UNIT  AREA  VS. 
TEMPERATURE  IN  SLOW- SLIDING  WEAR  TESTS  (MLO  7789A) 


Oil  Soluble  Iron,  yg 


1000 


□ MLO  7789A  Slow-Sliding 


A MLO  7789  Slow-Sliding 


O MLO  7789A  Static 


Figure  36.  OIL  SOLUBLE  IRON  VS.  THE  RECIPROCAL  OF  THE 
ABSOLUTE  TEMPERATURE. 


O MLO  7789A  Static 


A MLO  7789A  Slow- 

Sliding 

□ MLO  7789  Slow- 

Sliding 


OIL  SOLUBLE  IRON  PER  UNIT  AREA  VS.  THE 
RECIPROCAL  OF  THE  ABSOLUTE  TEMPERATURE 


ABSOLUTE  TEMPERATURE. 


**  BDC  = Below  Detectable 


Figure  42.  MEAN  WEAR  SCAR  DIAMETER  AS  A FUNCTION  OF  TIME 
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OIL  INSOLUBLE  ORGANIC  IRON  AS  A FUNCTION  OF 
TEMPERATURE  IN  THE  STATIC  TESTS  (MLO  7789A) 


Insoluble  Organic 


660  F ± 15  F 


Temperature,  F 

Figure  45, 

OIL  INSOLUBLE  ORGANIC  IRON  AS  A 
TEMPERATURE  IN  THE  SLOW-SLIDING 

FUNCTION  OF 
TESTS  (MLO  7789A) 

mhi t'Wl/  .*5 


Figure  47 


COMPARISON  BETWEEN  BOS'  EXPERIMENTAL  AND 
FRANCIS'  THEORETICAL  JUNCTION  TEMPERATURES 
IN  A FOUR-BALL  APPARATUS, 


D,  Conclusions . An  analysis  of  the  viscosity-pressure  measurement  In 
the  PRL  high  pressure  viscometer  shows  that  gas  solubility  In  the  fluid,  as  a 
result  of  diffusion  at  the  gas-liquid  interface,  represents  the  principal  experi- 
mental problem.  Studies  show  that  diffusion  rate  of  the  gas  into  the  liquid 
varies  over  a significant  range  and  appears  to  be  related  to  the  fluid  class. 
Significant  errors  can  be  shewn  after  only  one  hour  soak  time  at  pressures  as 
low  as  3000  psig.  The  gas  solubility  results  in  too  low  a measured  viscosity  at 
a given  pressure  Polyolefins  which  show  a significant  diffusion  rate  are  among 
the  hydrocarbons  that  show  a poor  predicted  pressure-viscosity  value.  In  order 
to  reduce  the  error  due  to  gas  solubility,  the  limiting  slope  of  the  viscosity- 
pressure  curve  at  0 psig  should  be  used  to  predict  the  viscosities  at  high 
pressures . 

Two  empirical  correlations  were  developed  from  the  100'F  viscosity-pressure 
data  from  the  PRL  pressure  viscometer  and  the  Penn  State  API  Project  42.  The 
final  form  of  the  equations  was  obtained  by  regression  analysis. 


For  fluids  with  average  of  about  4 5 formula  weights  per  atom.  Equation  [4  ] 
or  Figures  16  through  18  adequately  describe  the  pressure  coefficient  of  mineral 
oils,  pure  hydrocarbons,  nonhydrocarbcns  and  blends  of  mineral  oils,  polymer,  and 
resin  at  temperatures  ranging  from  32v  to  275' F. 

a = 0.8382  + 2.8563  V3,0627 

■+  0 8338  ASTM5,1903  V1-5976 
- 2.7574  V3  0975  D°-U62 


[4  ) 


where: 


a 


pressure  coefficient,  psig  1 x 10  \ 


V = kinematic  atmospheric  viacosity  at  temperature  of 
int  crest , cs, 

ASTM  « ASTM  slope,  and 

D - atmospheric  density  at  temperature  of  interest,  gm./cc. 


When  atmospheri:  density  is  nuc  available,  Equation  [4  ] or  Figure  22  can  be  used 
ter  a close  approximation.  For  fluids  with  average  formula  weights  far  from  4.5 


per  atom,  Equation 
coef  f ic lent . 


is  a better  equation  t_  u=e  for  approximating  the  pressure 

.1.0459 


0. 7236  + 0.2919  V 


+ 0.5860  ASTM 


3 5371  „i.670i 


l ) 


where: 


pressure  coefficient,  psig  1 x 10 


-.-4 
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kinematic  atmospheric  viscosity  at  temperature  of 
interest,  cs,  and 

ASTM  =*  ASTM  slope. 


The  parameters  used  for  these  correlations  appear  to  be  related  to  such 
physical  characteristics  as  molecular  interlocking,  molecular  packing,  molecular 
rigidity,  and  size  of  the  flow  unit.  The  behavior  of  the  correlations  was  shown 
to  be  in  agreement  with  the  void  volume  theory. 

A new  experimental  technique  has  been  developed  to  study  the  chemical 
reactions  in  boundary  lubrication.  The  method  combines  the  use  of  the  four-ball 
wear  tester  and  atomic  absorption  spectroscopy  with  successive  solvent  extractions. 
It  offers  new  insight  into  the  mechanism  of  lubrication  in  the  boundary  regime. 

The  junction  temperature  rise  has  been  determined  using  a chemical  conversion 
correlation  between  static  and  dynamic  runs.  It  is  substantially  higher  than  the 
generally  used  Blok-Archard  flash  temperatures. 

The  catalytic  aspect  of  metal  surfaces  in  chemical  reactions  under  boundary 
lubrication,  at  least  in  paraffinic  mineral  oils,  appears  to  be  insignificant. 

Chemical  reactions  do  indeed  play  a very  important  role  in  boundary 
lubrication. 

Oxygen  has  been  demonstrated  to  have  a large  effect  on  the  nature  of 
chemical  reaction  products. 
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Surface  Area  Calculation  of  the  3all-Pot  Assembl 


1.  Balls:  Surface  area/ball 


Total  exposed  surface  area 


2.  Bottom  piece:  d = 0.85",  h = 0.18"  circular 

Top  surface  = TTd2/4  = 0.567  in^ 
Lateral  sides  = 2(iTdh)  = 0.961  in' 


3.  Lock  ring:  di  = 0.9",  do  = 1.73",  h = 0.3"  donut  shape 

Lateral  sides  = 2iTdh  = 3.26  in2 
Surface  area  = 2ff(di/2  + do/2) (do/2  - di/2) 

= 3.428  in2 


4.  Ball-Pot 


Top  half  = 0.68"/2  = 0.34 
Surface  area  = irdh  = tt  (1 . 


Incline  surface 


Net  surface  area  of  incline 


From  trigonometry,  in  triangle  ABC 


I 

I 


in  ACFD 


< C = < A = 38.5' 


0.425 


= 0,425  tan  C = 0.338  in. 


Surface  area  = it  (r^/r^+h^  - r^r^+h,,2) 

= ti(0.865  / [0. 865) 2 + (0.688)2 

- 0.425  /(0.4252  + (0.338)2 


= 2.278  in^ 


. Total  surface  area  = 2.75  + 0.567  + 0.961  + 3.26  + 3.428 

+ 1.6  + 2.278 

= 14.844  in^ 

= 9576.75  mm2 


Surface  Area  Calculation  of  the  Wear  Particles 


A method  for  the  measurement  of  the  particle  diameters  of  the  wear 
debris  is  discussed  in  AFML-TR- 74-201,  Part  I.  From  this  information,  a 
particle  size  distribution  can  be  calculated  and  assuming  all  particles 
are  spheres,  a surface  area  distribution  can  be  estimated.  From  density 
data,  by  assuming  the  smallest  particles  are  iron  oxides  and  the  layer 
particles  iron,  the  total  weight  of  the  distribution  can  be  calculated. 
The  end  result  is  surface  area  per  unit  weight  ratio.  This  constant 
equals  4565.52  mm  /mg  for  this  particular  particle  distribution.  The 
calculations  are  shown  in  "able  67. 


*PIL\.  f 9 -T • .y 


Sample  Calculation 


Run  38  (MLO  7789  slow-sliding  at  600°F) 

Data:  Oil  soluble  iron  = 128  Mg 

Pyridine  soluble  iron  = 80  yg 
Iron  = 94  yg 

Iron  as  iron  oxides  = 149  yg 

Assuming  the  amount  of  iron  oxides  are  a 50:50  mixture  of 

Fe  0 and  Fe„0. 

2 3 3 4 

Atomic  weight  of  Fe  * 55.85 
Atomic  weight  of  0 = 16 

Molecular  weight  of  = 159.7 

Molecular  weight  of  Fe^^  = 231.55 

.*.  Amount  of  iron  in  = 111.7/159.7 

Amount  of  iron  in  Fe  0.  = 167.55/231.55 

3 4 

.’.  To  get  the  total  weight  of  iron  oxides  from  the  weight  of 
measured  iron,  one  needs  to  multiply  by  the  reciprocal  of  the  per- 
centage 

For  Fe^,  it  is  159.7/11.7  = 1.42 
Fe304,  it  is  231.5/167.5  = 1.38 
.'.  For  50:50  mixture,  the  multiplier  = 1.40 

Amount  of  iron  oxides  = 149  x 1.4  = 208.6  yg 
Amount  of  iron  particles  = 94.0  yg 
Total  iron  = 302.6  yg 

From  the  previous  section,  the  surface  area  of  the  particles  can  be 
estimated  from  the  measured  particle  size  distribution.  From  Table  67 
Surface  area/mg  = 4565.52 

.'.  Surface  area  of  the  iron  particles  = (317.5  yg) (4565. 52) (10  ^) 

= 1381.5  mm2 

Surface  area  of  the  ball-pot  assembly  = 9569.7  mm2 
.'.  Total  surface  area  = 1381.5  + 9569.7  = 10951.2  mm2 
.'.  Oil  soluble  Iron  = 128  yg/10951.2  = 11.688  x 10-2  yg/mm2 
Oil  insoluble  iron  = 80  yg/10951.2  = 7.306  x 10“  ^ yg/mm2 
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FLASH  TEMPERATURE  CALCULATIONS 


Appendix  II  of  Annual  Report 
AFML-TR-74-20I , Part  II,  Prepared  by 
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Calculation  of  Flash  Temperatures 
Data 


Runs  42-45 
Load  = 40  kg. 

Speed  = 600  r.p.m.  * 22.7  cm/sec. 

Temperature,  bulk  fluid  = 167°F 

Strain  gauge  force  = 215  gm 


Coefficient  of  Friction  Determination 


From  geometric  considerations,  the  force  normal  to  the  wear 
scar  Fjj  of  a stationary  ball  is  given  by: 


Fn  = 0.408  W 

= (0.408) (40kg)  = 16320  gm. 


(102) 


The  coefficient  of  friction,  f can  be  calculated  by 
f = 2.23  F £/W 


(103) 


where 


F 

SL 

W 


force  measured  by  strain  gauge,  gm. 

lever  arm  distance  from  the  center  of  rotation,  cm. 

machine  loading,  gm. 


(2. 23) (215) (10. 48)/(40,000)  = 0.117 


Blok-Archard  Theory 


N 


fg  11  pm 


PC 


(10«0 


where 


f = 

g = 
J = 


Pm  = 


P = 
C = 


r , gm-cm 
erg[=]  £ r 


coefficient  of  friction  = 0.117  (from  Run  42) 
acceleration  due  to  gravity  = 980  cm/sec^ 
mechanical  equivalent  of  heat  = 4.18  x 10?  ergs/ cal 
flow  pressure  = 7.04  x 10^  gm/cm^ 

density  = 7.78  gm/cm^ 
specific  heat  = 0.109  cal/gm°C 
2 


sec' 


dimensionally 


N[  = ] 


(cm/ sec^) (gm/cm^) 


2 2 3 

(gm-cm  )/sec  -cal) (gm/cm  )(cal/gm-°C 


[ = ] 


2 1 7 


Pjj  = normal  load  = 16320  gm. 


sliding  velocity  = 22.7  cm/ sec. 
thermal  dif fusitivity  = 0.12  cm^/sec. 
flow  pressure  of  steel  = 7.04  x 10^  gm/cm^ 


[=]  dimensionless 


dimensionally 


(16320) (22.7 


numerically 


(2)  (0.12)  [ (3. 1416X7. 04  x 10°)] 


Equation  (41)  applies 


6 is  estimated  from  Archard's  paper  (26)  to  be  about  0.75 


the  average  temperature 


0-117)(16200)(980)(22.7 


(3.1416) (7.04  x 10 


4(2-7  x 10-  ; (0-0546) 


Oc  is  given  by  Fig.  2 in  Archard's  paper  [Ref.  70]. 


Lc  = 2.57  from  previous  calculation 

0 

t ' °-85 

.".  0C  = 0.85  N = 0.85  (71.54) 

' " 0m  0B  0c  171  60.81 

0m  = 44.86’C  = 112. 7°C 

.*.  0m  = 167  + 112.7  = 279. 7°F 


= 60.81°C 
= 0.022 


f 


Blok's  formula  used  by  R,  S.  Fein: 

A spot  of  radius  r rubbing  against  an  extensive  surface  with 
uneven  heat  distribution  between  surfaces. 


0 % 

max 


2fP1 

pc 


<1  + A • ♦ 

1 + it/2  • 4> 


(108; 


where : 


W 

k 


c 

P 

♦ 

&o 


coefficient  of  friction  = 0.117 

junction  pressure,  kg/cm2  = 6.85  x 103  W3-/3 

= 2.343  x 104  kg/cm2 

machine  load,  kg  = 40  kg. 

thermal  conductivity  = 0.108  cal/cm-°C-sec 

= 4.609  kg-cm/cm=,C-sec 
specific  heat  = 0.109  cal/gn-’C 

= 4.652  x lQl  kg-cm/kg-°C 
density  of  steel  = 7.78  x 10~3  kg/cm3 
thermal  dif fusitivity  = 0.127  cm V sec 
ambient  temperature  = 167°F 


substituting  these  values  into  (B— 7 ) 


q n.  2(0. 117)(2.343  x 104)  (1 . 707)  (0 . 67)  ,, 

max  ° ^ ^ < 

(7.78  x 10  )(4.652  x 10  )(1  + 1.0259) 
0max~0°  ^ 85'5’c  = 186°F 


. 0n,ax  = 167  + 186  = 353°F 
or  as  in  Fein's  work,  0max  is  adjusted  by  1.33 

. *.  Qmax-Qo  ^ 186°  x 1-33  = 248°F 
0max  = 248  + 167°  = 415°F 


Francis'  Equation: 


For 

the  stationary  surface. 

A0S 

_ 9"  JL 
32  irrk 

(109) 

where 

Q = 

IT  = 

r = 

k = 

rate  of  heat  supply  = 1.0087  cal/sec 
3.1416 

wear  scar  radius  = 0.015  cm 
0.108  cal/cm-°C-sec 

•’  • 

A0S= 

175°C 

For  the  moving  surface 

A0v 

= 1.016  Vt  [1  + + ^372] 

iTrkB1’2  B B3/2 

(110) 

where 

B = 

Vr/Q  = (22. 7) (0.15)/ (0.127)  = 2.68 

A0V 

= 192. 8°C 

Therefore 

the 

mean  surface  temperature  0m  is  given  by 

AQm 

= i/(|-  + i-) 

Us  Uv 

( HI  ) 

A°m 

= 91. 7°C  = 197°F 

, ’ . 

®m 

= 167°  + 197°  = 364°F 

The  maximum  temperature  is  given  by 


A 0, 


0.852)^) 


max 


(112  ) 


1.996-1. 09 1B_0-818-0.537B_0- 2 L+B1/2 
Equation  (B-ll)  is  suitable  to  use  when  B >_  10.  In  this  work,  B = 2.68, 
so  Equation  (B-ll),  straightly  speaking,  is  not  suitable.  However  the 
error  introduced  this  way  is  small  compared  with  other  assumptions . 


A0. 


max 


Q, 


max 


137°C  = 275°F 
275  + 167  = 442#F 
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Bos'  Equation: 


I 


where 


AG, 


'm 


kr 


0.281 


(B+A) 


1/2 


Q = rate  of  heat  supply  = 4.22  Joules/sec 
k = thermal  conductivity  = 45.2  Joules/m3-sec-#K 
B = dimensionless  parameter  = 2.68 
r = radius  of  wear  scar  = 1.5  x 10-<!l  m 


(113)) 


A0m  = 68°K  = 154°F 

.“.  0m  = 167  + 154  = 321°F 

m 

This  result  differs  from  that  calculated  from  Francis'  Equations  (109), 
(110),  and  (111).  Francis  Equations  are  really  limiting  case  studies 
and  apply  at  B > 10.  Bos'  Equation  (113)  is  obtained  by  regression 
analysis  from  complete  numerical  solutions.  The  two  answers  are  close 
enough  to  justify  both  sets  of  equations. 

From  Figure  47  (Fig.  7 in  Bos'  work),  the  theoretical  and 
experimental  temperatures  for  r = 0.015  cm  are  estimated. 


Experimental:  EN  31  Steel 


k = 26.6  J/m3-sec-°K 

p = 7.81  x 103  kg/m3 

c = 460  J/kg-°K 

Q = 4.22  J/sec 


Theoretical:  A0  ^ „„ 

Q 

A0  ^ (23) (4.22)  = 97°K 

A0  experimental  _ 177^  1 
A0  theoretical  97 


Therefore  the  temperatures  calculated  from  previous  equations  are 
similarly  adjusted  to  determine  the  kind  of  flash  temperatures  that 
would  be  obtained  reiying  on  Bos'  findings 
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